Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “intermolecular potentials”

Search indexed NASA NTRS and DOE OSTI research on propulsion, heat transfer, battery materials and energy systems. Follow report and document links to the original sources.

Quote a phrase for an exact phrase match. Source license links do not imply unrestricted reuse.

123 records · Page 7

The Effective Fragment Molecular Orbital Method: Achieving High Scalability and Accuracy for Large Systems

The effective fragment molecular orbital (EFMO) method has been developed to predict the total energy of a very large molecular system accurately (with respect to the underlying quantum mechanical method) and efficiently by taking advantage of the locality of strong chemical interactions and employing a two-level hierarchical parallelism. The accuracy of the EFMO method is partly attributed to the accurate and robust intermolecular interaction prediction between distant fragments, in particular, the many-body polarization and dispersion effects, which require the generation of static and dynamic polarizability tensors by solving the coupled perturbed Hartree–Fock (CPHF) and time-dependent HF (TDHF) equations, respectively. Solving the CPHF and TDHF equations is the main EFMO computational bottleneck due to the inefficient (serial) and I/O-intensive implementation of the CPHF and TDHF solvers. In this work, the efficiency and scalability of the EFMO method are significantly improved with a new CPU memory-based implementation for solving the CPHF and TDHF equations that are parallelized by either message passing interface (MPI) or hybrid MPI/OpenMP. Here, the accuracy of the EFMO method is demonstrated for both covalently bonded systems and noncovalently bound molecular clusters by systematically examining the effects of basis sets and a key distance-related cutoff parameter, R cut . R cut determines whether a fragment pair (dimer) is treated by the chosen ab initio method or calculated using the effective fragment potential (EFP) method (separated dimers). Decreasing the value of Rcut increases the number of separated (EFP) dimers, thereby decreasing the computational effort. It is demonstrated that excellent accuracy (<1 kcal/mol error per fragment) can be achieved when using a sufficiently large basis set with diffuse functions coupled with a small R cut value. With the new parallel implementation, the total EFMO wall time is substantially reduced, especially with a high number of MPI ranks. Given a sufficient workload, nearly ideal strong scaling is achieved for the CPHF and TDHF parts of the calculation. For the first time, EFMO calculations with the inclusion of long-range polarization and dispersion interactions on a hydrated mesoporous silica nanoparticle with explicit water solvent molecules (more than 15k atoms) are achieved on a massively parallel supercomputer using nearly 1000 physical nodes. In addition, EFMO calculations on the carbinolamine formation step of an amine-catalyzed aldol reaction at the nanoscale with explicit solvent effects are presented.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Long-Range Fit: A Software Package for the Representation and Study of Long-Range Molecular Interactions

Describing intermolecular forces is fundamental to modeling and predicting the behavior of molecular systems. In particular, long-range molecular interactions─with electrostatic, induction, and dispersion as the main components─play a critical role, especially for low-temperature and low-density regimes. Long-range interactions are often described through perturbation theory, representing the electronic charge distribution via a multipolar series of the moments and polarizability tensors corresponding to each molecule. However, while the theory is well established, obtaining the resulting analytical expressions (and their practical implementation) constitutes a highly complex and system-dependent task. To address this challenge, we developed long-range-fit (LRF), an interactive and user-friendly software package designed to automate the generation and fitting of long-range interaction terms for arbitrary molecules in nondegenerate (ground or excited) electronic states. We have derived and implemented all terms up to 15th order, without approximations, via a spherical tensor representation, with symmetry adaptation to all molecular point-group symmetries. The resulting potential energy surface is compatible with most representations of the close interaction region.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Unveiling the transport properties of protic ionic liquids: Lithium ion dynamics modulated by the anion fluorine reservoir

Protic ionic liquids (PILs) show great potential as electrolyte components for energy storage devices. A comprehensive understanding of their transport properties must be achieved to optimize the design of safer and efficient electrolytes. This study focuses on a series of PILs based on the DBUH + cation (protonated 1,8-diazabicyclo[5,4,0]–undec-7-ene superbase) and three anions derived from strong acids: TFO – (triflate), IM14 – (perfluorobutyl-trifluoromethylsulfonylimide) and TFSI – (bis(trifluoromethylsulfonyl)imide). Neat PILs and PILs doped with LiTFO, LiIM14, and LiTFSI were studied using temperature-dependent NMR diffusion and relaxation techniques. The ionicity of these systems was also evaluated. Results revealed that the dynamic behaviour of lithium ions, as well as ionicity, strongly depend on the structural features of the anions, particularly in the case of IM14 – , whose main feature is the uneven distribution of the fluorinated sidegroups. The 19 F relaxation rates in IM14 – provide insights into the rotational reorientation of that anion. DBUH-IM14 exhibited diffusion coefficients lower than the expected ones on the basis of its viscosity, likely due to fluorophilic intermolecular interactions involving the fluorinated terminal groups. The presence of Li + in the DBUH-IM14 electrolyte led to unexpected and relatively faster translational mobility of Li + ions, resulting in a higher lithium apparent transference number. However, the trends observed in ionicity indicate a more complex interplay between intermolecular interactions and ion correlations. While DBUH-TFSI showed minimal effect of Li + addition, DBUH-TFO and DBUH-IM14 exhibited a significant decrease in ionicity, possibly attributed to strong interactions between ions.

25 ENERGY STORAGE↗

Internally Catalyzed Hydrogen Atom Transfer (I-CHAT)—A New Class of Reactions in Combustion Chemistry

The current paradigm of low-T combustion and autoignition of hydrocarbons is based on the sequential two-step oxygenation of fuel radicals. The key chain-branching occurs when the second oxygenation adduct (OOQOOH) is isomerized releasing an OH radical and a key ketohydroperoxide (KHP) intermediate. The subsequent homolytic dissociation of relatively weak O–O bonds in KHP generates two more radicals in the oxidation chain leading to ignition. Based on the recently introduced intramolecular “catalytic hydrogen atom transfer” mechanism (J. Phys. Chem. 2024, 128, 2169), abbreviated here as I-CHAT, we have identified a novel unimolecular decomposition channel for KHPs to form their classical isomers—enol hydroperoxides (EHP). The uncertainty in the contribution of enols is typically due to the high computed barriers for conventional (“direct”) keto–enol tautomerization. Remarkably, the I-CHAT dramatically reduces such barriers. The novel mechanism can be regarded as an intramolecular version of the intermolecular relay transfer of H-atoms mediated by an external molecule following the general classification of such processes (Catal. Rev.-Sci. Eng. 2014, 56, 403). Here, we present a detailed mechanistic and kinetic analysis of the I-CHAT-facilitated pathways applied to n-hexane, n-heptane, and n-pentane models as prototype molecules for gasoline, diesel, and hybrid rocket fuels. We particularly examined the formation kinetics and subsequent dissociation of the γ-enol-hydroperoxide isomer of the most abundant pentane-derived isomer γ-C5-KHP observed experimentally. To gain molecular-level insight into the I-CHAT catalysis, we have also explored the role of the internal catalyst moieties using truncated models. All applied models demonstrated a significant reduction in the isomerization barriers, primarily due to the decreased ring strain in transition states. In addition, the longer-range and sequential H-migration processes were also identified and illustrated via a combined double keto–enol conversion of heptane-2,6-diketo-4-hydroperoxide as a potential chain-branching model. To assess the possible impact of the I-CHAT channels on global fuel combustion characteristics, we performed a detailed kinetic analysis of the isomerization and decomposition of γ-C5-KHP comparing I-CHAT with key alternative reactions—direct dissociation and Korcek channels. Calculated rate parameters were implemented into a modified version of the n-pentane kinetic model developed earlier using RMG automated model generation tools (ACS Omega, 2023, 8, 4908). Simulations of ignition delay times revealed the significant effect of the new pathways, suggesting an important role of the I-CHAT pathways in the low-T combustion of large alkanes.

Biochemistry & Molecular Biology↗

Many-Body Effects in Aqueous Systems: Synergies Between Interaction Analysis Techniques and Force Field Development

Interaction analysis techniques, including the many-body expansion (MBE), symmetry-adapted perturbation theory, and energy decomposition analysis, allow for an intuitive understanding of complex molecular interactions. We review these methods by first providing a historical context for the study of many-body interactions and discussing how nonadditivities emerge from Hamiltonians containing strictly pairwise-additive interactions. We then elaborate on the synergy between these interaction analysis techniques and the development of advanced force fields aimed at accurately reproducing the Born–Oppenheimer potential energy surface. In particular, we focus on ab initio–based force fields that aim to explicitly reproduce many-body terms and are fitted to high-level electronic structure results. These force fields generally incorporate many-body effects through ( a) parameterization of distributed multipoles, ( b) explicit fitting of the MBE, ( c) inclusion of many-atom features in a neural network, and ( d) coarse-graining of many-body terms into an effective two-body term. We also discuss the emerging use of the MBE to improve the accuracy and speed of ab initio molecular dynamics.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

X-ray crystallographic and hydrogen deuterium exchange studies confirm alternate kinetic models for homolog insulin monomers

Despite the crucial role of various insulin analogs in achieving satisfactory glycemic control, a comprehensive understanding of their in-solution dynamic mechanisms still holds the potential to further optimize rapid insulin analogs, thus significantly improving the well-being of individuals with Type 1 Diabetes. Here, we employed hydrogen-deuterium exchange mass spectrometry to decipher the molecular dynamics of newly modified and functional insulin analog. A comparative analysis of H/D dynamics demonstrated that the modified insulin exchanges deuterium atoms faster and more extensively than the intact insulin aspart. Additionally, we present new insights derived from our 2.5 Å resolution X-ray crystal structure of modified hexamer insulin analog at ambient temperature. Furthermore, we obtained a distinctive side-chain conformation of the Asn3 residue on the B chain (AsnB3) by operating a comparative analysis with a previously available cryogenic rapid-acting insulin structure (PDB_ID: 4GBN). The experimental conclusions have demonstrated compatibility with modified insulin’s distinct cellular activity, comparably to aspart. Additionally, the hybrid structural approach combined with computational analysis employed in this study provides novel insight into the structural dynamics of newly modified and functional insulin vs insulin aspart monomeric entities. It allows further molecular understanding of intermolecular interrelations driving dissociation kinetics and, therefore, a fast action mechanism.

59 BASIC BIOLOGICAL SCIENCES↗

Molecular To Mesoscale Targeting of Oxoanions with Multi-Tasking Hosts

Achieving a better understanding of anion interactions both in solution and crystalline state was the overarching goal of this project. Anions are everywhere throughout Nature and play important roles in biological and environmental processes. They can be beneficial or deleterious or both in different situations and concentrations. For either reason it is important to have molecules that can bind anions for key needs that benefit society. However, recognition of specific anions is challenging due to the diffuse nature of their negative charge(s) as well as their various shapes and sizes. Understanding the basic properties of anions and how they interact with other molecules and ions in surrounding environments is key to selective recognition. In this project multi-tasking molecules for selective binding of targeted anions were designed to achieve cooperativity and synergism in one rather than multiple host molecules, including (1) cation:anion pair hosts for anions with charges of -2 or greater; (2) pH and redox activated hosts for on-off binding and release; and (3) multiple anion capture in extended host networks. Our design strategy was to combine the use of simple inexpensive building blocks and high yield synthetic pathways to provide economically feasible scale-up for applications. Oxoanions representing multiple shapes and charges were chosen based on having the potential for significant impact on DOE separations needs. Amide/amine-based macrocycles and urea/amine-based chelates and macrocycles with multiple hydrogen bonding sites provided the basic anion-binding frameworks. Successful multi-tasking outcomes were forthcoming in all three tasks. In Task 1, successful ion pair binding for anions with multiple charges was achieved. Furthermore, the ion pair molecules were capable of extended interactions through supramolecular intertwining, like fishing nets for capturing pools of fish (also fitting with Task 3). In Task 2, molecules were synthesized possessing on-off switches. These included a pH sensitive sensor for on-off binding of anions in general, as well as an electrochemical sensor selective for sulfate capture. Three new classes of extended anion host networks capable of binding multiple ions was a major outcome of Task 3. These systems included: anion sensitive, fluorescent organogels; channel-forming macrocycles for studying anion-water including larger macrocyclic cluster sandwiches; and, the offshoot of Task 1, fishing net ion-pair networks for higher valent anions. These strategies can be expanded in the future to other ions and molecules for a better understanding of intermolecular and interionic interactions.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Exploring the Ion Solvation Environments in Solid-State Polymer Electrolytes through Free-Energy Sampling

The success of poly(ethylene oxide) (PEO) in solid-state polymer electrolytes for lithium-ion batteries is well established. Recently, in order to understand this success and to explore possible alternatives, we studied polyacetal electrolytes to deepen the understanding of the effect of the local chemical structure on ion transport. Advanced molecular dynamics techniques using newly developed, tailored interaction potentials have helped elucidate the various coordination environments of ions in these systems. In particular, the competition between cation-Anion pairing and coordination by the polymer has been explored using free-energy sampling (metadynamics). At equivalent reduced temperatures, with respect to the polymer-specific glass-Transition temperature, two-dimensional free-energy plots reveal the existence of multiple coordination environments for the lithium (Li) ions in these systems and their relative stabilities. Furthermore, we observe that the Li-ion movement in PEO follows a serial, stepwise pathway when moving from one coordination state to another, whereas this happens in a more continuous and concerted fashion in a polyacetal such as poly(1,3-dioxalane) [P(EO-MO)]. The implication is that interconversion between coordination states of the Li ions may be easier in P(EO-MO). However, the overarching observation from our free-energy analysis is that Li-ion coordination is dominated by the polymer (in either case) and contact-ion pairs are rare. We rationalize the observed higher increase in glass-Transition temperature (Tg) with salt loading in polyacetals as due to intermolecular Li-ion coordination involving multiple polymer chains, rather than just one chain for PEO-based electrolytes. This interchain coupling in the polyacetals, resulting in the higher Tg, works against any gains due to variations in Li-ion coordination that might enhance transport processes over PEO. Further research is required to overcome the interdependence between local coordination and macroscopic properties to compete with PEO electrolytes at the same absolute working temperature.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Impacts of hydrogen bonding interactions with Np(v/vi)O 2 Cl 4 complexes: vibrational spectroscopy, redox behavior, and computational analysis

The neptunyl (Np(V)O 2 + /Np(VI)O 2 2+ ) cation is the dominant form of 237 Np in acidic aqueous solutions and the stability of the Np(V) and Np(VI) species is driven by the specific chemical constituents present in the system. Hydrogen bonding with the oxo group may impact the stability of these species, but there is limited understanding of how these intermolecular interactions influence the behavior of both solution and solid-state species. In the current study, we systematically evaluate the interactions between the neptunyl tetrachloride species and hydrogen donors in coordination complexes and in the related aqueous solutions. Both Np(V) compounds (N 2 C 4 H 12 ) 2 [Np(V)O 2 Cl 4 ]Cl (Np(V)pipz) and (NOC 4 H 10 ) 3 [Np(V)O 2 Cl 4 ] (Np(V)morph) exhibit directional hydrogen bonding to the neptunyl oxo group while Np(VI) compounds (NC 5 H 6 ) 2 [Np(VI)O 2 Cl 4 ] (Np(VI)pyr) and (NOC 4 H 10 ) 4 [Np(VI)O 2 Cl 4 ]·2Cl (Np(VI)morph) assemble via halogen interactions. Here, the Raman spectra of the solid-state phases indicate the activation of vibrational bands when there is asymmetry of the neptunyl bond, while these spectral features are not observed within the related solution phase spectra. Density functional theory calculations of the Np(V)pipz system suggest that activation of the ν 3 asymmetric stretch and other combination modes lead to additional complexity within the solid-state spectra. Electrochemical analyses of complexes in the solution phases are consistent with the results of the crystallization experiments as the voltammetric potentials of Np(V)/Np(VI) complexes in the presence of protonated heterocycles differ from the potentials of pure Np(V) and may correlate with the hydrogen bonding interactions.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Probing the influence of ion-pairing on ligand-field excited-state dynamics

Exploration of the photophysical and photochemical properties of transition metal complexes has driven ground-breaking advancements in solar energy conversion technologies, including photoredox catalysis. While significant research has been devoted to understanding excited state properties of second- and third-row transition metal complexes, earth-abundant first-row metal complexes have received comparatively little attention in this context until very recently. In particular, the role of ion-pairing – which has been identified as a potentially significant factor for Ir(III)-based photosensitizers – has not been examined with regard to its influence on the ligand-field excited states that dominate much of first-row photophysics. A key challenge in studying ion-pair interactions lies in quantifying the extent and nature of ion-pairing, particularly in non-aqueous media where the vast majority of photophysical studies are performed. Cobalt(III) polypyridyl complexes provide an attractive platform to address such questions due to their demonstrated potential for applications in photoredox catalysis involving ligand-field excited states. In the present study, we prepared a cobalt(III) polypyridyl complex, [Co(4,4′-OMebpy) 3 ](BAr F 4 ) 3 (where 4,4′-OMebpy is 4,4′-dimethoxy-2,2′-bipyridine and BAr F 4 is tetrakis(3,5-bis(trifluoromethyl)-phenyl)borate) to probe ion-pairing in non-aqueous solutions. Specifically, analysis of data acquired from both variable-temperature diffusion ordered spectroscopy (DOSY) NMR and 1-D rotating-frame nuclear Overhauser effect (ROE) experiments allowed us to identify and differentiate between solvent-separated ion pairs in high-dielectric media and contact ion pairs in a low-dielectric solvent. Time-resolved absorption spectroscopy was then used to measure ground-state recovery dynamics under these varying conditions of ion-pairing, the results of which revealed an increase in excited-state lifetime for contact ion-pairs that we suggest arises from a reduction in outer-sphere reorganization energy relative to conditions which favored solvent-separated ion pairs. We believe this study demonstrates that one can leverage broadly available NMR-based methods to understand ion-pairing in non-aqueous solutions, which in turn can provide a microscopic picture of intermolecular interactions that can impact the photophysical properties of transition metal-based chromophores.

Ghosh, Atanu [Michigan State University, East Lans↗

Predictive coarse-grained (CG) modeling of morphologies in polymer nanocomposites with specific and directional intermolecular interactions (Final Report)

The overarching goal of the proposed work was to develop predictive models for investigating structure and dynamics in soft materials with chemistries that have specific and directional molecular interactions. The motivation behind studying materials with specific and directional interactions lies in the many desirable features these interactions provide when designing novel soft materials. Soft materials with specific and direction interactions (such as hydrogen bonds or H-bonds) can have a) thermally reversible phase behavior with different functions with varying temperature, b) precisely tuned nanostructure with desirable geometries that afford unique physical properties (e.g., color response, mechanical properties) and c) well-mixed/blended morphologies that are useful for variety of applications in energy field (e.g., materials for batteries require use of blended polymers where one domain gives superior mechanical properties and one domain promotes electrical conduction). Notably, biology makes extensive use of specific and directional interactions, in many cases based on H-bonds, to construct materials with precisely defined architectures and properties. Engineered soft materials with precisely tuned nanostructures and improved processiblity through thermoresponsive phase behavior are useful in numerous applications that are relevant to the Department of Energy (DOE) including high efficiency electronic devices, light-weight high-strength composite materials for batteries and fuel cells, and polymer membranes for separations, etc. While past computational studies have been tremendously useful in understanding molecular phenomena and guiding synthesis of new macromolecular soft materials for a wide variety of applications, the inability to capture small scale specific and directional interactions alongside macromolecular length and time scales represented a key limitation of most studies to date. Our work in this project addressed this grand challenge in computational materials chemistry, i.e., ability to model the anisotropic, directional, and specific interactions that govern the behavior of many macromolecular soft matter systems of interest, thus, has greatly expanded the predictive potential of simulations. Specifically the key outcomes were: successful development of new coarse-grained (CG) polymer models to study generic and specific polymer chemistries in which hydrogen-bonding interactions are dominant. These CG models were then used in molecular simulations to study structure and thermodynamics in polymer nanocomposites and blends; some studies were conducted in collaboration with experimentalists. We also published a perspective and a viewpoint which included some of the work we completed in this DOE project; we believe these perspective and viewpoint articles guide other researchers in the soft materials community on how to extend the computational approaches and models we have developed for the purposes of their studies.

36 MATERIALS SCIENCE↗

Stockmayer fluid simulations for viscosity and glass transition temperature of ionic liquids

We develop a Stockmayer fluid model for molecular dynamics simulations of ionic liquids that captures molecular polarization, ionic conductivity, viscosity, and glass transition temperature, using ethylammonium nitrate (EAN) as an example. The ions in EAN are treated as spheres interacting via the Lennard-Jones potential with an embedded point charge and a permanent dipole moment. We show that our simulation results for EAN are consistent with experimental data and then explore the effects of the molecular parameters on the viscosity of ionic liquids. Our results indicate that viscosity monotonically increases with ionic charge and dipole moment but non-monotonically changes with ionic diameter (or molar volume). This non-monotonic trend arises from the competition among the electrostatic interactions, molecular packing, and size asymmetry between the cation and anion. In conclusion, our model also shows that long-lived ion pairs result in higher viscosities.

Coarse-grained simulations↗

Influence of counterion substitution on the properties of imidazolium-based ionic liquid clusters

Due to their unique physiochemical properties that may be tailored for specific purposes, ionic liquids (ILs) have been investigated for various applications, including chemical separations, catalysis, energy storage, and space propulsion. The different cations and anions comprising ILs may be selected to optimize a range of desired properties, such as thermal stability, ionic conductivity, and volatility, leading to the designation of certain ILs as designer “green” solvents. The effect of counterions on the properties of ILs is of both fundamental scientific interest and technological importance. Herein, we report a systematic experimental and theoretical investigation of the size, charge, stability toward dissociation, and geometric/electronic structure of 1-ethyl-3-methyl imidazolium (EMIM)-based IL clusters containing two different atomic counterions (i.e., bromide [Br − ] and iodide [I − ]). This work extends our studies of EMIM + cations with atomic chloride (Cl − ) and molecular tetrafluoroborate (BF 4 − ) anions reported previously by Baxter et al. [Chem. Mater. 34, 2612 (2022)] and Zhang et al . [J. Phys. Chem. Lett. 11, 6844 (2020)], respectively. Distributions of anionic IL clusters were generated in the gas phase using electrospray ionization and characterized by high mass resolution mass spectrometry, energy-resolved collision-induced dissociation, and negative ion photoelectron spectroscopy experiments. The experimental results reveal anion-dependent trends in the size distribution, relative abundance, ionic charge state, stability toward dissociation, and electron binding energies of the IL clusters. Complementary global optimization theory provides molecular-level insights into the bonding and electronic structure of a selected subset of clusters, including their low energy structures and electrostatic potential maps, and how these fundamental characteristics are influenced by anion substitution. Collectively, our findings demonstrate how the fundamental properties of ILs, which determine their suitability for many applications, may be tuned by substituting counterions. These observations are critical in the sub-nanometer cluster size regime where phenomena do not scale predictably to the bulk phase, and each atom counts toward determining behavior.

cluster↗

Crystal Engineering of Hydrogen Bonding for Direct Air Capture of CO2: A Quantum Crystallography Perspective

Rising atmospheric CO2 levels demand efficient and sustainable carbon capture solutions. Direct air capture (DAC) via crystallizing hydrogen-bonded frameworks such as carbonate salts has emerged as a promising approach. This review explores the potential of crystal engineering, in tandem with advanced quantum crystallography techniques and computational modeling, to unlock the full potential of DAC materials. We examine the critical role of hydrogen bonding and other noncovalent interactions within a family of bis-guanidines that governs the formation of carbonate salts with high CO2 capture capacity and low regeneration energies for utilization. Quantum crystallography and charge density analysis prove instrumental in elucidating these interactions. A case study of a highly insoluble carbonate salt of a 2,6-pyridine-bis-(iminoguanidine) exemplifies the effectiveness of these approaches. However, challenges remain in the systematic and precise determination of hydrogen atom positions and atomic displacement parameters within DAC materials using quantum crystallography, and limitations persist in the accuracy of current energy estimation models for hydrogen bonding interactions. Future directions lie in exploring diverse functional groups, designing advanced hydrogen-bonded frameworks, and seamlessly integrating experimental and computational modeling with machine learning. This synergistic approach promises to propel the design and optimization of DAC materials, paving the way for a more sustainable future.

36 MATERIALS SCIENCE↗

Tuning the Intermolecular Electrostatic Interaction toward High-Efficiency and Low-Cost Organic Solar Cells

Organic solar cells (OSCs) have achieved much progress with rapidly increasing power conversion efficiencies (PCEs). It should be noted that the top-performance OSCs are generally consisted of active materials with complex chemical structures, resulting in high costs. Here, combining the material design and morphology control, high-efficiency OSCs are fabricated by a low-cost donor: acceptor blend. A completely non-fused electron acceptor named Tz is designed and synthesized via introducing thiazole units on both sides of a bithiophene core, which shows an outstanding PCE of 13.3% with a typical polythiophene donor. More importantly, optimization guidelines are presented to get excellent morphology for low-cost donor:acceptor systems. Three polythiophenes are selected, poly(3-hexylthiophene) and its two derivatives with electron-withdrawing substitutions (PDCBT and PDCBT-2F), as donors to fabricate the cell devices. The computational and experimental data reveal that decreasing the electrostatic interaction between polythiophene and Tz is the key to getting a suppressed miscibility and thus a high phase purity. Finally, this study provides insight into the molecular design and donor:acceptor matching requirements for high-efficiency and low-cost OSCs.

14 SOLAR ENERGY↗