A New Approach for Simulating Inhomogeneous Chemical Kinetics | Research Square window.SnipcartSettings = { analytics: { enabled: false } }; (function() { var accessVector = localStorage.getItem('access_vector') || ''; window.dataLayer = window.dataLayer || []; if (accessVector) { window.dataLayer.push({ user: { profile: { profileInfo: { snid: accessVector } } } }); } })(); (function(w,d,s,l,i){w[l]=w[l]||[];w[l].push({'gtm.start':new Date().getTime(),event:'gtm.js'});var f=d.getElementsByTagName(s)[0],j=d.createElement(s),dl=l!='dataLayer'?'&l='+l:'';j.async=true;j.src='https://www.googletagmanager.com/gtm.js?id='+i+dl;f.parentNode.insertBefore(j,f);})(window,document,'script','dataLayer','GTM-K279D39R'); Browse Preprints In Review Journals COVID-19 Preprints AJE Video Bytes Research Tools Research Promotion AJE Professional Editing AJE Rubriq About Preprint Platform In Review Editorial Policies Our Team Advisory Board Help Center Sign In Submit a Preprint Cite Share Download PDF Article A New Approach for Simulating Inhomogeneous Chemical Kinetics Georgia Bradshaw, Mel O'Leary, Arthur Purser, Balder Villagomez-Bernabe, and 3 more This is a preprint; it has not been peer reviewed by a journal. https://doi.org/ 10.21203/rs.3.rs-2752229/v1 This work is licensed under a CC BY 4.0 License Status: Published Journal Publication published 28 Aug, 2023 Read the published version in Scientific Reports → Version 1 posted 10 You are reading this latest preprint version Abstract In this paper, inhomogeneous chemical kinetics are simulated by describing the concentrations of interacting chemical species by a linear expansion of basis functions in such a manner that the coupled reaction and diffusion processes are propagated through time efficiently by tailor-made numerical methods. The approach is illustrated through modelling α− and γ−radiolysis in thin layers of water and at their solid interfaces from the start of the chemical phase until equilibrium was established. The method’s efficiency is such that hundreds of such systems can be modelled in a few hours using a single core of a typical laptop, allowing the investigation of the effects of the underlying parameter space. Illustrative calculations showing the effects of changing dose-rate and water-layer thickness are presented. Other simulations are presented which show the approach’s capability to solve problems with spherical symmetry (an approximation to an isolated radiolytic spur), where the hollowing out of an initial Gaussian distribution is observed, in line with previous calculations. These illustrative simulations show the generality and the computational efficiency of this approach to solving reaction-diffusion problems. Furthermore, these example simulations illustrate the method’s suitability for simulating solid-fluid interfaces, which have received a lot of experimental attention in contrast to the lack of computational studies. Physical sciences/Chemistry Physical sciences/Mathematics and computing Full Text Additional Declarations No competing interests reported. Cite Share Download PDF Status: Published Journal Publication published 28 Aug, 2023 Read the published version in Scientific Reports → Version 1 posted Editorial decision: Major revision 29 Jun, 2023 Reviews received at journal 20 Jun, 2023 Reviews received at journal 23 May, 2023 Reviewers agreed at journal 17 May, 2023 Reviewers agreed at journal 16 May, 2023 Reviewers invited by journal 16 May, 2023 Editor assigned by journal 16 May, 2023 Editor invited by journal 18 Apr, 2023 Submission checks completed at journal 18 Apr, 2023 First submitted to journal 29 Mar, 2023 You are reading this latest preprint version Research Square lets you share your work early, gain feedback from the community, and start making changes to your manuscript prior to peer review in a journal. As a division of Research Square Company, we’re committed to making research communication faster, fairer, and more useful. We do this by developing innovative software and high quality services for the global research community. Our growing team is made up of researchers and industry professionals working together to solve the most critical problems facing scientific publishing. Also discoverable on Platform About Our Team In Review Editorial Policies Advisory Board Help Center Resources Author Services Accessibility API Access RSS feed Manage Cookie Preferences © Research Square 2026 | ISSN 2693-5015 (online) Privacy Policy Terms of Service Do Not Sell My Personal Information {"props":{"pageProps":{"initialData":{"identity":"rs-2752229","acceptedTermsAndConditions":true,"allowDirectSubmit":false,"archivedVersions":[],"articleType":"Article","associatedPublications":[],"authors":[{"id":192995660,"identity":"3d5c52ea-388b-4387-abfc-ab796c11aa0a","order_by":0,"name":"Georgia Bradshaw","email":"data:image/png;base64,iVBORw0KGgoAAAANSUhEUgAAAZAAAAAyAQMAAABI0h/eAAAABlBMVEX///8AAABVwtN+AAAACXBIWXMAAA7EAAAOxAGVKw4bAAABWElEQVRIie2RMWvCQBSA3xFolkuznmjbv3AhYAfF/BWPQFwMFApdFAwU4hJ3S6C/oUUQx8BBXDIXxaG6ODm0UMRAS5tgKCa1nTvcN9y9O97He48HIBD8Q1QnPSmAkVwEJS81/UDOQQ7JKyTIFM3JlJKTKfgXhX4HQaakwd+KfBu+vl01QPP5ah5PGmz0xNfLeFI3DHmwhF0X2J2TV3Bo+hVqgnFq6TUlMtl4YV1qg8hiHp5S5IXA/HwVbdbWJUIl0DyolpEr2eNFMw14ExMLQHGA3Rcae96kSg9oJG9LsduzR35rmyifBr5YA/o4osywjl4oTxRcJYrL7YdyO60SII+cgJRWKTRGIsuUgE6x5uHrmuJOe8NF+6Y0cE3mRRbwSkj0wvhqn3O0e++cG1h+nMduR1f91pjEbsOQ+yFabbr1s2EARSR8sIQ8wY+t7EG74/kCgUAg2PMFHdlzIDjNMr8AAAAASUVORK5CYII=","orcid":"","institution":"University of Manchester","correspondingAuthor":true,"prefix":"","firstName":"Georgia","middleName":"","lastName":"Bradshaw","suffix":""},{"id":192995661,"identity":"5c284ece-5777-42a2-baf6-49bbfb44cb3b","order_by":1,"name":"Mel O'Leary","email":"","orcid":"","institution":"Dalton Cumbrian Facility","correspondingAuthor":false,"prefix":"","firstName":"Mel","middleName":"","lastName":"O'Leary","suffix":""},{"id":192995662,"identity":"5eef1817-55ab-41f5-8e61-15030b48166e","order_by":2,"name":"Arthur Purser","email":"","orcid":"","institution":"University of Manchester","correspondingAuthor":false,"prefix":"","firstName":"Arthur","middleName":"","lastName":"Purser","suffix":""},{"id":192995663,"identity":"421fb9d7-e3c2-45f2-8dda-328b81a35e3d","order_by":3,"name":"Balder Villagomez-Bernabe","email":"","orcid":"","institution":"Royal Surrey County Hospital","correspondingAuthor":false,"prefix":"","firstName":"Balder","middleName":"","lastName":"Villagomez-Bernabe","suffix":""},{"id":192995664,"identity":"d2737263-3e53-4b62-8601-8731cf5d5909","order_by":4,"name":"Cyrus Wyett","email":"","orcid":"","institution":"University of Manchester","correspondingAuthor":false,"prefix":"","firstName":"Cyrus","middleName":"","lastName":"Wyett","suffix":""},{"id":192995665,"identity":"99af9f44-4c2b-416e-80a6-a5ca17a08169","order_by":5,"name":"Frederick Currell","email":"","orcid":"","institution":"University of Manchester","correspondingAuthor":false,"prefix":"","firstName":"Frederick","middleName":"","lastName":"Currell","suffix":""},{"id":192995666,"identity":"7b3e654c-ca76-47ae-a464-4400e24f8afc","order_by":6,"name":"Marcus Webb","email":"","orcid":"","institution":"University of Manchester","correspondingAuthor":false,"prefix":"","firstName":"Marcus","middleName":"","lastName":"Webb","suffix":""}],"badges":[],"createdAt":"2023-03-29 13:44:38","currentVersionCode":1,"declarations":"","doi":"10.21203/rs.3.rs-2752229/v1","doiUrl":"https://doi.org/10.21203/rs.3.rs-2752229/v1","draftVersion":[],"editorialEvents":[{"content":"https://doi.org/10.1038/s41598-023-39741-y","type":"published","date":"2023-08-28T15:10:28+00:00"}],"editorialNote":"","failedWorkflow":false,"files":[{"id":42782373,"identity":"6d85af70-d96f-416b-82a0-7bf2b1fb6ab9","added_by":"auto","created_at":"2023-09-07 15:16:56","extension":"pdf","order_by":1,"title":"","display":"","copyAsset":false,"role":"manuscript-pdf","size":474281,"visible":true,"origin":"","legend":"","description":"","filename":"ANewApproachforSimulatingInhomogeneousChemicalKinetics.pdf","url":"https://assets-eu.researchsquare.com/files/rs-2752229/v1_covered_fb26fd80-b50d-4738-81fc-eb2c3f6fb525.pdf"}],"financialInterests":"No competing interests reported.","formattedTitle":"A New Approach for Simulating Inhomogeneous Chemical Kinetics","fulltext":[],"fulltextSource":"","fullText":"","funders":[],"hasAdminPriorityOnWorkflow":false,"hasManuscriptDocX":false,"hasOptedInToPreprint":true,"hasPassedJournalQc":"","hasAnyPriority":false,"hideJournal":false,"highlight":"","institution":"","isAcceptedByJournal":true,"isAuthorSuppliedPdf":true,"isDeskRejected":"","isHiddenFromSearch":false,"isInQc":false,"isInWorkflow":false,"isPdf":true,"isPdfUpToDate":true,"isWithdrawnOrRetracted":false,"journal":{"display":true,"email":"
[email protected]","identity":"scientific-reports","isNatureJournal":false,"hasQc":true,"allowDirectSubmit":false,"externalIdentity":"scirep","sideBox":"Learn more about [Scientific Reports](http://www.nature.com/srep/)","snPcode":"","submissionUrl":"","title":"Scientific Reports","twitterHandle":"","acdcEnabled":true,"dfaEnabled":true,"editorialSystem":"stoa","reportingPortfolio":"Scientific Reports","inReviewEnabled":true,"inReviewRevisionsEnabled":true},"keywords":"","lastPublishedDoi":"10.21203/rs.3.rs-2752229/v1","lastPublishedDoiUrl":"https://doi.org/10.21203/rs.3.rs-2752229/v1","license":{"name":"CC BY 4.0","url":"https://creativecommons.org/licenses/by/4.0/"},"manuscriptAbstract":"In this paper, inhomogeneous chemical kinetics are simulated by describing the concentrations of interacting chemical species by a linear expansion of basis functions in such a manner that the coupled reaction and diffusion processes are propagated through time efficiently by tailor-made numerical methods. The approach is illustrated through modelling α− and γ−radiolysis in thin layers of water and at their solid interfaces from the start of the chemical phase until equilibrium was established. The method’s efficiency is such that hundreds of such systems can be modelled in a few hours using a single core of a typical laptop, allowing the investigation of the effects of the underlying parameter space. Illustrative calculations showing the effects of changing dose-rate and water-layer thickness are presented. Other simulations are presented which show the approach’s capability to solve problems with spherical symmetry (an approximation to an isolated radiolytic spur), where the hollowing out of an initial Gaussian distribution is observed, in line with previous calculations. These illustrative simulations show the generality and the computational efficiency of this approach to solving reaction-diffusion problems. Furthermore, these example simulations illustrate the method’s suitability for simulating solid-fluid interfaces, which have received a lot of experimental attention in contrast to the lack of computational studies.","manuscriptTitle":"A New Approach for Simulating Inhomogeneous Chemical Kinetics","msid":"","msnumber":"","nonDraftVersions":[{"code":1,"date":"2023-04-20 04:05:23","doi":"10.21203/rs.3.rs-2752229/v1","editorialEvents":[{"type":"communityComments","content":0},{"type":"decision","content":"Major revision","date":"2023-06-29T05:01:08+00:00","index":"","fulltext":""},{"type":"editorInvitedReview","content":"","date":"2023-06-20T08:11:04+00:00","index":"hide","fulltext":""},{"type":"editorInvitedReview","content":"","date":"2023-05-23T18:52:21+00:00","index":"hide","fulltext":""},{"type":"reviewerAgreed","content":"4600c742-fb2b-482b-b575-fb123c84808e","date":"2023-05-17T17:31:49+00:00","index":"hide","fulltext":""},{"type":"reviewerAgreed","content":"b6ebb8d4-fadb-4f48-9dd5-a1a4fd07df50","date":"2023-05-16T13:07:22+00:00","index":"hide","fulltext":""},{"type":"reviewersInvited","content":"","date":"2023-05-16T13:03:36+00:00","index":"","fulltext":""},{"type":"editorAssigned","content":"","date":"2023-05-16T06:18:21+00:00","index":"","fulltext":""},{"type":"editorInvited","content":"","date":"2023-04-18T15:20:19+00:00","index":"","fulltext":""},{"type":"checksComplete","content":"","date":"2023-04-18T14:50:00+00:00","index":"","fulltext":""},{"type":"submitted","content":"Scientific Reports","date":"2023-03-29T13:33:32+00:00","index":"","fulltext":""}],"status":"published","journal":{"display":true,"email":"
[email protected]","identity":"scientific-reports","isNatureJournal":false,"hasQc":true,"allowDirectSubmit":false,"externalIdentity":"scirep","sideBox":"Learn more about [Scientific Reports](http://www.nature.com/srep/)","snPcode":"","submissionUrl":"","title":"Scientific Reports","twitterHandle":"","acdcEnabled":true,"dfaEnabled":true,"editorialSystem":"stoa","reportingPortfolio":"Scientific Reports","inReviewEnabled":true,"inReviewRevisionsEnabled":true}}],"origin":"","ownerIdentity":"ad925275-f784-459b-af19-4eee8ab1da11","owner":[],"postedDate":"April 20th, 2023","published":true,"recentEditorialEvents":[],"rejectedJournal":[],"revision":"","amendment":"","status":"published-in-journal","subjectAreas":[{"id":20809577,"name":"Physical sciences/Chemistry"},{"id":20809578,"name":"Physical sciences/Mathematics and computing"}],"tags":[],"updatedAt":"2023-09-07T15:14:53+00:00","versionOfRecord":{"articleIdentity":"rs-2752229","link":"https://doi.org/10.1038/s41598-023-39741-y","journal":{"identity":"scientific-reports","isVorOnly":false,"title":"Scientific Reports"},"publishedOn":"2023-08-28 15:10:28","publishedOnDateReadable":"August 28th, 2023"},"versionCreatedAt":"2023-04-20 04:05:23","video":"","vorDoi":"10.1038/s41598-023-39741-y","vorDoiUrl":"https://doi.org/10.1038/s41598-023-39741-y","workflowStages":[]},"version":"v1","identity":"rs-2752229","journalConfig":"researchsquare"},"__N_SSP":true},"page":"/article/[identity]/[[...version]]","query":{"redirect":"/article/rs-2752229","identity":"rs-2752229","version":["v1"]},"buildId":"_2-kVJe1T_tPrBINL-cwx","isFallback":false,"isExperimentalCompile":false,"dynamicIds":[84888],"gssp":true,"scriptLoader":[]}
Text is read by the "Ask this paper" AI Q&A widget below.
Extraction quality varies by source — PMC NXML preserves structure
cleanly, OA-HTML may include some navigation residue, and OA-PDF can
have broken hyphenation. The publisher copy
(via DOI)
is the canonical version.