Bayesian optimization of solution-processed thermoelectric polymers using a database of ab initio electronic structure data

Bayesian optimization of solution-processed thermoelectric polymers using a database of ab initio electronic structure data | 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 Bayesian optimization of solution-processed thermoelectric polymers using a database of ab initio electronic structure data Connor Ganley, Howard E. Katz, Paulette Clancy This is a preprint; it has not been peer reviewed by a journal. https://doi.org/ 10.21203/rs.3.rs-5850990/v1 This work is licensed under a CC BY 4.0 License Status: Posted Version 1 posted You are reading this latest preprint version Abstract Historically, organic materials design takes place wherein the components in the system are chosen based on expert knowledge of organic synthetic chemistry. Experimental synthesis and processing, often followed by computational validation, is a standard route to design materials such as polymers. The advent of data-driven machine learning approach can now shift that paradigm to one in which computational information \textit{leads} experimental validation to select a much smaller, more promising and perhaps unexpected set of new materials. This work describes one such investigation to characterize novel, hypothetical, and promising polymer candidates by their calculated molecular-scale electronic and chemical parameters. Subsequently, we optimize them for a key chemical property using a chemically-informed Bayesian surrogate model representing those properties. As a test case, over 7300 combinations of novel, hypothetical semiconducting diketopyrrolopyrrole-based (DPP) polymers and commonly used solvents, generated by density functional theory, were screened based on their free energies of solvation, G$ {solv}$. From this synthetic data set, we trained a physics-informed Gaussian process model that linked molecular-scale electronic structure properties to G$ {solv}$, and then used Bayesian optimization via the PAL 2.0 software to identify key descriptors for predicting optimal polymer solvation. Implicit and explicit solvation models of these synthetic polymers, combined with Bayesian optimization, exhibited a minimum in $\Delta$G$ {solv}$ as a function of the solvent dielectric. As a result, we predict an ''optimal'' solvent dielectric (around 10) for the entire DPP-based polymer class. The chemically-informed nature of our approach allowed us to identify a molecular electronic property that demonstrates near-linear additivity when constructing an aggregate synthetic polymer from functional group sub-units. This linearly related property was the isotropic quadrupole moment, which we found to be highly correlated to $\Delta$G$ {solv}$. As validation of this chemical insight, we showed that this minimum $\Delta$G$ {solv}$ was associated with the highest experimentally measured conductivity among a set of synthesized doped polymers. The importance of these observations is that PAL 2.0 has provided the chemical insight that can now allow us to quickly screen solvents for potential DPP-based polymers based on the polymer repeat unit's quadrupole moments and to identify preferred compatible solvents corresponding to a minimized $\Delta$G$ {solv}$. This study also highlights the effectiveness of the PAL 2.0 algorithm to find the optimal solution using just 6% of the parameter space and in a very short execution time (on the order of minutes), a feat which cannot be duplicated experimentally. Physical sciences/Engineering/Chemical engineering Physical sciences/Physics/Condensed matter physics/Electronic properties and materials Physical sciences/Chemistry/Materials chemistry/Electronic materials Physical sciences/Materials science/Theory and computation/Electronic structure Full Text Additional Declarations No competing interests reported. Supplementary Files Polymerpapersi.pdf Cite Share Download PDF Status: Posted Version 1 posted 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-5850990","acceptedTermsAndConditions":true,"allowDirectSubmit":true,"archivedVersions":[],"articleType":"Article","associatedPublications":[],"authors":[{"id":411893140,"identity":"53cf674f-d9af-49ba-a9a3-491b4acc4021","order_by":0,"name":"Connor Ganley","email":"","orcid":"","institution":"The Johns Hopkins University","correspondingAuthor":false,"prefix":"","firstName":"Connor","middleName":"","lastName":"Ganley","suffix":""},{"id":411893142,"identity":"69c93ccd-fee9-4ad9-9750-79affecff7f9","order_by":1,"name":"Howard E. Katz","email":"","orcid":"","institution":"The Johns Hopkins University","correspondingAuthor":false,"prefix":"","firstName":"Howard","middleName":"E.","lastName":"Katz","suffix":""},{"id":411893147,"identity":"c40321d9-0c2d-4a7c-a378-182b7a0f4f6f","order_by":2,"name":"Paulette Clancy","email":"data:image/png;base64,iVBORw0KGgoAAAANSUhEUgAAAZAAAAAyAQMAAABI0h/eAAAABlBMVEX///8AAABVwtN+AAAACXBIWXMAAA7EAAAOxAGVKw4bAAAAx0lEQVRIiWNgGAWjYFAC5gYQacfPcAbMI0YLI1hLsmQDqVoYNxzgIVKLOfvBxs8FNYeZjQ+ePSbBUGGd2EBIi2VPYrP0jGOH+cwOnEuTYDiTTliLwQ3GBmnehsPMZgfOmEkwth0mSkvzb6AWxs0NIC3/iNPSBrKFcQMDSEsDEVqAfmmz5jmWnixx4IyxRcKxdGOCWszZDx++zVNjbcc/44zhjQ811rKEHQahmhkYJA4wMCQQUo6kpY6BgZ+g8aNgFIyCUTBSAQDBqEDM9BZVCgAAAABJRU5ErkJggg==","orcid":"","institution":"The Johns Hopkins University","correspondingAuthor":true,"prefix":"","firstName":"Paulette","middleName":"","lastName":"Clancy","suffix":""}],"badges":[],"createdAt":"2025-01-17 17:23:13","currentVersionCode":1,"declarations":"","doi":"10.21203/rs.3.rs-5850990/v1","doiUrl":"https://doi.org/10.21203/rs.3.rs-5850990/v1","draftVersion":[],"editorialEvents":[],"editorialNote":"","failedWorkflow":false,"files":[{"id":83092864,"identity":"ca1e0338-eaed-4843-ad60-50394b0b66e7","added_by":"auto","created_at":"2025-05-20 02:31:51","extension":"pdf","order_by":1,"title":"","display":"","copyAsset":false,"role":"manuscript-pdf","size":1448864,"visible":true,"origin":"","legend":"","description":"","filename":"Polymerpapermanuscript.pdf","url":"https://assets-eu.researchsquare.com/files/rs-5850990/v1_covered_ecaf3818-6b47-4fbd-8325-2eb0efa7a73a.pdf"},{"id":75701624,"identity":"2d5e1051-a951-48d0-9432-7c4d9a766566","added_by":"auto","created_at":"2025-02-07 09:28:00","extension":"pdf","order_by":1,"title":"","display":"","copyAsset":false,"role":"supplement","size":552164,"visible":true,"origin":"","legend":"","description":"","filename":"Polymerpapersi.pdf","url":"https://assets-eu.researchsquare.com/files/rs-5850990/v1/94a581510df7553353749f0c.pdf"}],"financialInterests":"No competing interests reported.","formattedTitle":"Bayesian optimization of solution-processed thermoelectric polymers using a database of ab initio electronic structure data","fulltext":[],"fulltextSource":"","fullText":"","funders":[],"hasAdminPriorityOnWorkflow":false,"hasManuscriptDocX":false,"hasOptedInToPreprint":true,"hasPassedJournalQc":"","hasAnyPriority":false,"hideJournal":true,"highlight":"","institution":"","isAcceptedByJournal":false,"isAuthorSuppliedPdf":true,"isDeskRejected":"","isHiddenFromSearch":false,"isInQc":false,"isInWorkflow":false,"isPdf":true,"isPdfUpToDate":true,"isWithdrawnOrRetracted":false,"journal":{"display":true,"email":"[email protected]","identity":"researchsquare","isNatureJournal":false,"hasQc":true,"allowDirectSubmit":true,"externalIdentity":"","sideBox":"","snPcode":"","submissionUrl":"/submission","title":"Research Square","twitterHandle":"researchsquare","acdcEnabled":true,"dfaEnabled":false,"editorialSystem":"","reportingPortfolio":"","inReviewEnabled":false,"inReviewRevisionsEnabled":true},"keywords":"","lastPublishedDoi":"10.21203/rs.3.rs-5850990/v1","lastPublishedDoiUrl":"https://doi.org/10.21203/rs.3.rs-5850990/v1","license":{"name":"CC BY 4.0","url":"https://creativecommons.org/licenses/by/4.0/"},"manuscriptAbstract":"Historically, organic materials design takes place wherein the components in the system are chosen based on expert knowledge of organic synthetic chemistry. Experimental synthesis and processing, often followed by computational validation, is a standard route to design materials such as polymers. The advent of data-driven machine learning approach can now shift that paradigm to one in which computational information \\textit{leads} experimental validation to select a much smaller, more promising and perhaps unexpected set of new materials. This work describes one such investigation to characterize novel, hypothetical, and promising polymer candidates by their calculated molecular-scale electronic and chemical parameters. Subsequently, we optimize them for a key chemical property using a chemically-informed Bayesian surrogate model representing those properties. As a test case, over 7300 combinations of novel, hypothetical semiconducting diketopyrrolopyrrole-based (DPP) polymers and commonly used solvents, generated by density functional theory, were screened based on their free energies of solvation, G$_{solv}$. From this synthetic data set, we trained a physics-informed Gaussian process model that linked molecular-scale electronic structure properties to G$_{solv}$, and then used Bayesian optimization via the PAL 2.0 software to identify key descriptors for predicting optimal polymer solvation. Implicit and explicit solvation models of these synthetic polymers, combined with Bayesian optimization, exhibited a minimum in $\\Delta$G$_{solv}$ as a function of the solvent dielectric. As a result, we predict an ''optimal'' solvent dielectric (around 10) for the entire DPP-based polymer class. The chemically-informed nature of our approach allowed us to identify a molecular electronic property that demonstrates near-linear additivity when constructing an aggregate synthetic polymer from functional group sub-units. This linearly related property was the isotropic quadrupole moment, which we found to be highly correlated to $\\Delta$G$_{solv}$. As validation of this chemical insight, we showed that this minimum $\\Delta$G$_{solv}$ was associated with the highest experimentally measured conductivity among a set of synthesized doped polymers. The importance of these observations is that PAL 2.0 has provided the chemical insight that can now allow us to quickly screen solvents for potential DPP-based polymers based on the polymer repeat unit's quadrupole moments and to identify preferred compatible solvents corresponding to a minimized $\\Delta$G$_{solv}$. This study also highlights the effectiveness of the PAL 2.0 algorithm to find the optimal solution using just 6\\% of the parameter space and in a very short execution time (on the order of minutes), a feat which cannot be duplicated experimentally.","manuscriptTitle":"Bayesian optimization of solution-processed thermoelectric polymers using a database of ab initio electronic structure data","msid":"","msnumber":"","nonDraftVersions":[{"code":1,"date":"2025-02-07 09:27:55","doi":"10.21203/rs.3.rs-5850990/v1","editorialEvents":[{"type":"communityComments","content":0}],"status":"published","journal":{"display":true,"email":"[email protected]","identity":"researchsquare","isNatureJournal":false,"hasQc":true,"allowDirectSubmit":true,"externalIdentity":"","sideBox":"","snPcode":"","submissionUrl":"/submission","title":"Research Square","twitterHandle":"researchsquare","acdcEnabled":true,"dfaEnabled":false,"editorialSystem":"","reportingPortfolio":"","inReviewEnabled":false,"inReviewRevisionsEnabled":true}}],"origin":"","ownerIdentity":"d2203126-fe15-4bf0-bb1d-2ce6e9de2775","owner":[],"postedDate":"February 7th, 2025","published":true,"recentEditorialEvents":[],"rejectedJournal":[],"revision":"","amendment":"","status":"posted","subjectAreas":[{"id":43920675,"name":"Physical sciences/Engineering/Chemical engineering"},{"id":43920676,"name":"Physical sciences/Physics/Condensed matter physics/Electronic properties and materials"},{"id":43920677,"name":"Physical sciences/Chemistry/Materials chemistry/Electronic materials"},{"id":43920678,"name":"Physical sciences/Materials science/Theory and computation/Electronic structure"}],"tags":[],"updatedAt":"2025-05-20T02:23:41+00:00","versionOfRecord":[],"versionCreatedAt":"2025-02-07 09:27:55","video":"","vorDoi":"","vorDoiUrl":"","workflowStages":[]},"version":"v1","identity":"rs-5850990","journalConfig":"researchsquare"},"__N_SSP":true},"page":"/article/[identity]/[[...version]]","query":{"redirect":"/article/rs-5850990","identity":"rs-5850990","version":["v1"]},"buildId":"8U1c8b4HqxoKbykW_rLl7","isFallback":false,"isExperimentalCompile":false,"dynamicIds":[84888],"gssp":true,"scriptLoader":[]}