Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “DFT Calculations”

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.

At least 847 records · Page 47
All sourcesNASA NTRSDOE OSTI
DOE OSTI2026

First-Principles Insights into Proton-Coupled Electron Transfer versus Hydrogen Evolution Reaction Selectivity from a Base-Appended Cobaltocene Mediator

Performing selective proton-coupled electron transfer (PCET) to substrates such as N 2 , CO 2 , and unsaturated organic molecules under electrochemical conditions requires the suppression of the competing hydrogen evolution reaction (HER). To address this challenge, our laboratory previously demonstrated a PCET mediator strategy using a dimethylaniline-appended cobaltocene complex, [(CpCoCp NMe2 )H] + , which performs selective reductive chemistry while suppressing the HER. However, the origin of the suppressed, yet still observable, HER has not been thoroughly established. In this work, we perform density functional theory (DFT) calculations to elucidate the HER mechanism involving this redox mediator and to provide atomistic insights into the bifurcation between the PCET and HER pathways. We find that protonation of the aniline moiety to form [CpCoCp NMe2H ] + is more favorable, both kinetically and thermodynamically, than formation of the ring-protonated species [(CpCo(Cp-H) NMe2 )] + . Furthermore, PCET to acetophenone is energetically more favorable via [CpCoCp NMe2H ] + than via [(CpCo(Cp-H) NMe2 )] +1/0 . In contrast, the most favorable HER pathway involves the ring-protonated Co(I) species. These results offer mechanistic insights into HER versus PCET bifurcation and establish guiding principles for designing PCET mediators for selective electroreductive transformations.

evolution reactions↗
DOE OSTI2020

Dehydration of UO 2 Cl 2 ·3H 2 O and Nd(NO 3 ) 3 ·6H 2 O with a Soft Donor Ligand and Comparison of Their Interactions through X-ray Diffraction and Theoretical Investigation [plus Supplemental Information]

Herein we investigated whether the relatively Lewis basic imidazole-2-thiones could be used to substitute water ligands bound to f-element cations and generate f-element soft donor complexes. Reactions of 1,3-diethylimidazole-2-thione (C 2 C 2 ImT) with Nd(NO 3 ) 3 ·6H 2 O and UO 2 Cl 2 ·3H 2 O led to the isolation of the anhydrous thione complexes Nd(NO 3 ) 3 (C 2 C 2 ImT) 3 and UO 2 Cl 2 (C 2 C 2 ImT) 2 , characterized by single crystal X-ray diffraction. Differences in the strength of metal–thione interactions have been examined by means of the crystal structure analysis and density functional theory (DFT) calculations. The C 2 C 2 ImT ligands were found to be affected by both coordination and noncovalent interactions, making it impossible to deconvolute the effects of one from the other. Calculated partial atomic charges indicated greater ligand-to-metal charge transfer in the [UO 2 ] 2+ complex, indicative of a stronger interaction. The reactivity of C 2 C 2 ImT demonstrates its usefulness in the preparation of f-element soft donor complexes from readily available hydrates that could be useful intermediates for promoting the coordination and studying the effects of soft donor anions.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
DOE OSTI2025

A Comprehensive Machine Learning Model for Metal–Ligand Binding Prediction: Applications in Chemistry and Biology

A machine-learning (ML) model that predicts metal–ligand binding constants was developed using the open-source Chemprop software. The model was trained on over 30,000 experimental log K 1 values, which include both protonation and metal–ligand stability constants, comprising over 3500 ligands and 10 2 metal ions from 73 total elements, thus generalizing beyond existing limited approaches, which focus only on specific metals or ligand families. The best-performing model included a combination of SMILES-based molecular representations along with descriptors for the metal ion and experimental conditions. It had an external test R 2 value of 0.942, and MAE value of 0.834. A “SMILES-only” simpler version also produced accurate predictions and preserved the binding trends, serving as a quick and easily accessible alternative for users without computational expertise. The SMILES-only model performed comparably to density functional theory (DFT) calculations but utilized a fraction of the computational resources. The model was successfully applied across diverse domains, including bioinorganic chemistry, heavy metal remediation, and sensor development and demonstrated its effectiveness as a rapid and reliable screening tool for both academic and industrial uses.

Ligands↗
DOE OSTI2020

Reconciling Work Functions and Adsorption Enthalpies for Implicit Solvent Models: A Pt (111)/Water Interface Case Study

Implicit solvent models are a computationally efficient method of representing solid/liquid interfaces prevalent in electrocatalysis, energy storage, and materials science. However, electronic structure changes induced at the metallic surface by the dielectric continuum are not fully understood. To address this, we perform DFT calculations for the Pt(111)/water interface, in order to compare Poisson–Boltzmann continuum solvation methods with ab initio molecular dynamics (AIMD) simulations of explicit solvent. We show that the implicit solvent cavity can be parametrized in terms of the electric dipole moment change at the equilibrated explicit Pt/water interface to obtain the potential of zero charge (PZC). We also compare the accuracy of aqueous enthalpies of adsorption of phenol on Pt(111) using geometry and charge density based dielectric cavitation methods. The ability to parametrize the cavity according to individual atoms, as afforded in the geometry based approach, is key to obtaining accurate enthalpy changes of adsorption under aqueous conditions. Additionally, we show that the electronic structure changes induced by explicit solvent and our proposed implicit solvent parametrization scheme yield comparable density difference profiles and d-band projected density of states. We therefore demonstrate the capability of implicit solvent approaches to capture both the energetics of adsorption processes and the main electronic effects of aqueous solvent on the metallic surface. Therefore, this work provides a scheme for computationally efficient simulations of interfacial processes for applications in areas such as heterogeneous catalysis and electrochemistry.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
DOE OSTI2021

Semi-Automated Creation of Density Functional Tight Binding Models through Leveraging Chebyshev Polynomial-Based Force Fields

Density functional tight binding (DFTB) is an attractive method for accelerated quantum simulations of condensed matter due to its enhanced computational efficiency over standard density functional theory (DFT) approaches. However, DFTB models can be challenging to determine for individual systems of interest, especially for metallic and interfacial systems where different bonding arrangements can lead to significant changes in electronic states. In this regard, we have created a rapid-screening approach for determining systematically improvable DFTB interaction potentials that can yield transferable models for a variety of conditions. Our method leverages a recent reactive molecular dynamics force field where many-body interactions are represented by linear combinations of Chebyshev polynomials. This allows for the efficient creation of multi-center representations with relative ease, requiring only a small investment in initial DFT calculations. Here, we have focused our workflow on TiH 2 as a model system and show that a relatively small training set based on unit-cell-sized calculations yields a model accurate for both bulk and surface properties. Our approach is easy to implement and can yield reliable DFTB models over a broad range of thermodynamic conditions, where physical and chemical properties can be difficult to interrogate directly and there is historically a significant reliance on theoretical approaches for interpretation and validation of experimental results.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
DOE OSTI2023

Even Faster Exact Exchange for Solids via Tensor Hypercontraction

Hybrid density functional theory (DFT) remains intractable for large periodic systems due to the demanding computational cost of exact exchange. We apply the tensor hypercontraction (THC) (or interpolative separable density fitting) approximation to periodic hybrid DFT calculations with Gaussian-type orbitals using the Gaussian plane wave approach. This is done to lower the computational scaling with respect to the number of basis functions (N) and k-points (N k ) at a fixed system size. Additionally, we propose an algorithm to fit only occupied orbital products via THC (i.e., a set of points, N ISDF ) to further reduce computation time and memory usage. This algorithm has linear scaling cost with k-points, no explicit dependence of N ISDF on basis set size, and overall cubic scaling with unit cell size. Significant speedups and reduced memory usage may be obtained for moderately sized k-point meshes, with additional gains for large k-point meshes. Adequate accuracy can be obtained using THC-oo-K for self-consistent calculations. In conclusion, we perform illustrative hybrid density function theory calculations on the benzene crystal in the basis set and thermodynamic limits to highlight the utility of this algorithm.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
DOE OSTI2024

Efficient Parameterization of Density Functional Tight-Binding for 5 f -Elements: A Th–O Case Study

Density functional tight binding (DFTB) models for f-element species are challenging to parametrize owing to the large number of adjustable parameters. The explicit optimization of the terms entering the semiempirical DFTB Hamiltonian related to f orbitals is crucial to generating a reliable parametrization for f-block elements, because they play import roles in bonding interactions. However, since the number of parameters grows quadratically with the number of orbitals, the computational cost for parameter optimization is much more expensive for the f-elements than for the main group elements. In this work we present a set of efficient approaches for mitigating the hurdle imposed by the large size of the parameter space. A novel group-by-orbital correction functions for two-center bond integrals was developed. With this approach the number of parameters is reduced, and it grows linearly with the number of elements, maintaining the accuracy and the number of parameters, in the case of f elements, by more than 40%. The parameter optimization step was accelerated by means of the mini-batch BFGS method. This method allows parameter optimizations with much larger training sets than other single batch methods. A stochastic optimizer was employed that helped overcome shallow local minima in the objective function. The proposed algorithm was used to parametrize the DFTB Hamiltonian for the Th–O system, which was subsequently applied to the study of ThO 2 nanoparticles. The training set consisted of 6322 unique structures, which is barely feasible with conventional optimization methods. The optimized parameter set, LANL-ThO, displays good agreement with DFT-calculated properties such as energies, forces, and structures for both clusters and bulk ThO 2 . Benefiting from the fewer number of parameters and lower computational costs for objective function evaluations, this new approach shows its potential applications in DFTB parametrization for elements with high angular momentum, which present a challenge to conventional methods.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
DOE OSTI2025

Wavelength-Selective Reactivity of Iron(III) Halide Salts in Photocatalytic C–H Functionalization

The utility of halogen radicals in hydrocarbon functionalization extends from early examples of photochemical halogenation to recent reports using photoredox catalysis with iridium complexes and simple transition metal salts such as FeCl 3 . The majority of these methods (uncatalyzed and iron-catalyzed) require UV light (λ ≤ 390 nm), and systematic efforts to enable the use of visible light remain valuable. We report the use of a simple Fe(III) salt that enables a C–H to C–C and C–N functionalization under visible light. The reactivity and selectivity profile using different light sources demonstrates wavelength-selective behavior, which was further investigated with deuterium kinetic isotope effect experiments and DFT calculations. These results show that control over the reactive intermediates in this iron-catalyzed reaction can be achieved through proper choice of the wavelength of irradiation.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
DOE OSTI2020

Phase Change-Driven Negative Activation Energies in Pd/Carbon-Based/Organic Getter Hydrogenation Reactions

In this study, the hydrogenation of 1,4-diphenylbutadiyne (DPB) blended with carbon-supported Pd (DPB–Pd/C) in the form of pellets was investigated by isothermal–isobaric experiments at 1333 Pa of H 2 and in the temperature range of 291–315 K. The extracted kinetics were then used in conjunction with a complementary constant rate of H 2 input experimentation to model the performance of a DPB-catalysis/support system as a function of temperature and H 2 partial pressure. First-principles density functional theory (DFT) calculations were also performed to shed light on the molecular level energetics of DPB and its intermediate states. A seemingly puzzling formation of alternate positive activation energy barrier (higher reaction rate with higher temperature) and negative activation energy barrier (higher reaction rate with lower temperature) zones during the hydrogenation process was discovered. However, this observed phenomenon can be logically explained in terms of the associated phase changes and H 2 transport in the material. This work provides a good illustration of a rarely encountered chemical process with a negative activation energy barrier.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
DOE OSTI2022

Proton-Coupled Electron Transfer in a Ruthenium(II) Bipyrimidine Complex in Its Ground and Excited Electronic States

In this work, proton-coupled electron transfer (PCET) was studied for the ground and excited electronic states of a [Ru(terpy)(bpm)(OH 2 )(PF 6 ) 2 ] complex, Ru-bpm. Cyclic voltammetry measurements show that the Ru(II)-aqua moiety undergoes PCET to form a Ru(IV)-oxo moiety in the anodic region, while the bpm ligand undergoes PCET to form bpmH 2 in the cathodic region. The photophysical behavior of Ru-bpm was studied using steady-state and femtosecond transient UV–vis absorption spectroscopy, coupled with density functional theory (DFT) calculations. The lowest-lying excited state of Ru-bpm is described as a (Ru → bpm) metal-to-ligand charge-transfer (MLCT) state, while the metal-centered (MC) excited state was found computationally to be close in energy to the lowest-energy bright MLCT state (MC state was 0.16 eV above the MLCT state). The excited-state kinetics of Ru-bpm were found via transient absorption spectroscopy to be short-lived and were fit well to a biexponential function with lifetimes $τ_1$ = 4 ps and $τ_2$ = 65 ps in aqueous solution. Kinetic isotope effects of 1.75 ($τ_1$) and 1.61 ($τ_2$) were observed for both decay components, indicating that the solvent plays an important role in the excited-state dynamics of Ru-bpm. Based on the pH-dependent studies and the results from prior studies of similar Ru-complexes, we hypothesize that the 3 MLCT state forms an excited-state hydrogen-bond adduct with the solvent molecules and that this process occurs with a 4 ps lifetime. The formation of such a hydrogen-bond complex is consistent with the electronic density accumulation at the peripheral N atoms of the bpm moiety in the 3 MLCT state. The hydrogen-bonded state 3 MLCT decays to the ground state with a 65 ps lifetime. Such a short lifetime is likely associated with the efficient vibrational energy transfer from the 3 MLCT state to the solvent.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
DOE OSTI2022

Photophysical Properties and Electronic Structure of Hydroporphyrin Dyads Exhibiting Strong Through-Space and Through-Bond Electronic Interactions

Electronic interactions between tetrapyrroles are utilized in natural photosynthetic systems to tune the light-harvesting and energy/charge-transfer processes in these assemblies. Such interactions also can be employed to tailor the electronic properties of tetrapyrrolic dyads and larger arrays for use in materials science and biomedical research. Here we have utilized static and time-resolved optical spectroscopy to characterize the optical absorption and emission properties of a set of chlorin and bacteriochlorin dyads with varied degrees of through-bond (TB) and through-space (TS) interactions between the constituent macrocycles. The dyads consist of two chlorins or two bacteriochlorins joined by a linker that utilizes a triple-double-triple bond (enediyne) motif in which the double bond portion is a ester-substituted ethylene or o-phenylene unit. The photophysical studies are coupled with density functional theory (DFT) calculations to probe the ground-state molecular orbital (MO) characteristics of the dyads and time-dependent calculations (TDDFT) to elucidate excited-state properties. The latter include electronic characteristics of the singlet excited state manifold and the absorption transitions to these states from the electronic ground state. Comparison of the MO and calculated spectral properties of each dyad with the linker present versus disrupted (by eliminating the double-bond portion) gives insight into the relative contributions of TB versus TS interactions to the electronic properties of the dyads. The results show that the TB and TS contributions are additive (constructively interfere), which is not always the case for molecular dyads. Most of the dyads have shorter lifetimes of the lowest singlet excited state compared to the parent monomer, which derives from increased S 1 → S 0 internal conversion. The enhancement is greater for the dyads in benzonitrile than in toluene. Furthermore, the studies provide insights into the nature of the electronic interactions between the constituents in the tetrapyrrole arrays and how these interactions dictate the spectral properties and excited-state decay characteristics.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
DOE OSTI2022

Computational Predictions of the Hydrolysis of 2,4,6-Trinitrotoluene (TNT) and 2,4-Dinitroanisole (DNAN)

Hydrolysis is a common transformation reaction that can affect the environmental fate of many organic compounds. In this study, three proposed mechanisms of alkaline hydrolysis of 2,4,6-trinitrotoluene (TNT) and 2,4-dinitroaniline (DNAN) were investigated with plane-wave density functional theory (DFT) combined with ab initio and classical molecular dynamics (AIMD/MM) free energy simulations, Gaussian basis set DFT calculations, and correlated molecular orbital theory calculations. Most of the computations in this study were carried out using the Arrows web based tools. For each mechanism: Meisenheimer complex formation, nucleophilic aromatic substitution, and proton abstraction reaction energies and activation barriers were calculated for the reaction at each relevant site. For TNT, it was found that the most kinetically favorable first hydrolysis steps involve Meisenheimer complex formation by attachment of OH$^{-}$ at the C1 and C3 arene carbons and proton abstraction from the methyl group. The nucleophilic aromatic substitution reactions at the C2 and C4 arene carbons were found to be thermodynamically favorable. However, the calculated activation barriers were slightly lower than previous studies, but still found to be $\Delta$G$^{\ddagger} \approx 18$ kcal/mol using PBEo AIMD/MM free energy simulations, suggesting the reactions are not kinetically significant. For DNAN, the barriers of nucleophilic aromatic substitution were even greater ($\Delta$G$^{\ddagger} > 29$ kcal/mol PBEo AIMD/MM). Further, the most favorable hydrolysis reaction for DNAN was found to be a two-step process in which the hydroxyl first attacks the C1 carbon to form a Meisenheimer complex at the C1 arene carbon C1-(OCH$_3$)OH$^{-}$, subsequently, the methoxy anion (-OCH$_3$) at the C1 arene carbon dissociates and the proton shuttles from the C1-OH to the dissociated methoxy group, resulting in methanol and an aryloxy anion.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
DOE OSTI2025

Structural Changes in Metal Chalcogenide Nanoclusters Associated with Single Heteroatom Incorporation

Atomically precise nanoclusters (NCs) are promising building blocks for designing materials and interfaces with unique properties. By incorporating heteroatoms into the core, the electronic and magnetic properties of NCs can be precisely tuned. To accurately predict these properties, density functional theory (DFT) is often employed, making the rigorous benchmarking of DFT results particularly important. In this study, we present a benchmarking approach based on metal chalcogenide NCs as a model system. We synthesized a series of bimetallic, iron-cobalt chalcogenide NCs [Co 6-x Fe x S 8 (PEt 3 ) 6 ] + (x = 0-6) (PEt = triethyl phosphine) and investigated the effect of heteroatoms in the octahedral metal chalcogenide core on their size and electronic properties. Using ion mobility-mass spectrometry (IM-MS), we observed a gradual increase in the collision cross section (CCS) with an increase in the number of Fe atoms in the core. DFT calculations combined with trajectory method CCS simulations successfully reproduced this trend, revealing that the increase in cluster size is primarily due to changes in metal-ligand bond lengths, while the electronic properties of the core remain largely unchanged. Moreover, this method allowed us to exclude certain multiplicity states of the NCs, as their CCS values were significantly different from those predicted for the lowest-energy structures. Here, this study demonstrates that gas-phase IM-MS is a powerful technique for detecting subtle size differences in atomically precise NCs, which are often challenging to observe using conventional NC characterization methods. Accurate CCS measurements are established as a benchmark for comparison with theoretical calculations. The excellent correspondence between experimental data and theoretical predictions establishes a robust foundation for investigating structural changes of transition metal NCs of interest to a broad range of applications.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
DOE OSTI2025

How Thermodynamic, Electronic, and Steric Factors Influence Mesitylcopper Oligomers

Mesitylcopper (CuMes) is a highly versatile organocopper reagent used in both organic and inorganic syntheses. It has previously been shown that CuMes exists as a tetrameric or pentameric cyclic oligomer [CuMes] n (n = 4, 5), both in solution and in the solid state. The bonding arrangement between the [CuMes] units has qualitatively been described as localized three-center twoelectron (3c-2e) bonds. However, the electronic, structural, and thermodynamic forces driving this aggregation are still not well understood. For this reason, we employed density functional theory (DFT) calculations to study mesitylcopper as a monomeric [CuMes] unit and [CuMes]n oligomers with n = 2 to n = 7. We found that there is a strong electronic driving force for aggregation caused by strong mixing between the Cu’s d orbitals and Mes’s π orbitals in oligomers larger than the dimer. This mixing is only optimized in oligomers with n ≥ 3, where the mesityl group is no longer bonded to a single copper center but instead becomes a bridging ligand. Beyond the trimer, steric and entropic factors become relevant for determining the relative stabilities of the different aggregates, with midsized oligomers (n = 4−5) having the optimal balance between the electronic Cu−C bonding character, Cu···Cu attractive forces, entropy, reduced internal ring strain, and reduced steric interactions between the mesityl groups.

Copper↗
DOE OSTI2026

Exploring Excited State Proton Transfer Dynamics upon Ultraviolet Excitation

Excited state intramolecular proton transfer (ESIPT) has been investigated in two prototypical systems─salicylaldehyde azine (SAA) and 1,5-dihydroxyanthraquinone (DHAQ)─using transient absorption spectroscopy upon ultraviolet excitation into the less studied higher excited (S n ) manifold. Excitation with sub-30 fs pulses and broadband visible probing has allowed for direct measurement of the ESIPT rate. In conjunction with steady-state measurements and TD-DFT calculations, a complete delineation of the ultrafast photophysics has been carried out. In SAA, ESIPT remains ultrafast (∼30 fs), consistent with previous S 1 excitation studies. Coherent vibrational beats maps reveal significant wavelength dependence, however. Theoretical analysis suggests that the observed modes and their intensities in coherent vibrational spectra are modulated by the nature of the electronically excited state. In DHAQ, the first direct observation of ESIPT presents a time-constant of ∼85 fs, and a slower component of 9 ps, akin to previous reports on double-proton transfer systems. Collectively, the results suggest that while the ultrafast ESIPT rate remains largely invariant vis-à-vis the excitation energy, the product yield, as well as accompanying coherent oscillations, may be substantively altered, owing to the existence of alternative decay pathways.

excited states↗
DOE OSTI2020

Deciphering the Solvent Effect for the Solvation Structure of Ca 2+ in Polar Molecular Liquids

Although the crystal structures for many inorganic compounds are readily available, researchers are still working hard to understand the relations between the structures and chemical properties of solutions because most of the chemical reactions take place in solutions. A huge amount of effort has been put toward modeling the ion solvation structure from the perspectives of both experiments and theories. In this study, the solvation structures of Ca 2+ ions in aqueous and alcoholic solutions at different concentrations were carefully evaluated by Ca K-edge X-ray absorption near-edge structure (XANES) and extended X-ray absorption fine structure (EXAFS) analyses. Further, density functional theory (DFT) calculations were also performed to correlate the results with the experimental data and then further extended to other similar systems. It was found that the number of coordinating solvent molecules decreases with increasing Ca 2+ concentration and increasing solvent molecule sizes. From the EXAFS data, it was observed that the first solvation shell of Ca 2+ splits into two Ca-O distances in a methanol solution and the counter ion Cl - might also be within the first shell at high concentrations. For the first time, the effects of solvents with different polarities and sizes on the ion solvation environment were systematically evaluated.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
DOE OSTI2020

Extracting the Frequency-Dependent Dynamic Stokes Shift from Two-Dimensional Electronic Spectra with Prominent Vibrational Coherences

The dynamic Stokes shift is a common means for characterizing ultrafast solvation dynamics of electronically excited states. Here we extract the excitation frequency-dependent dynamic Stokes shift from two-dimensional electronic spectra (2DES) of cresyl violet, a molecule with a well-defined vibronic progression. The extracted dynamic Stokes shift function, S(t), exhibits oscillatory behavior, and the oscillatory components are assigned to intramolecular vibrational modes through DFT and TD-DFT calculations. Here, the well-characterized oscillations are incorporated into the fitting procedure of S(t). The excitation frequency dependence of the ultrafast response is examined through the analysis of S(t) obtained from slices taken at different excitation frequencies of the 2DES spectra. The extracted ultrafast timescales range from 36 to 98 fs, and we interpret the frequency dependence of the timescales in the context of other dynamic processes that also lead to lineshape changes in the 2DES spectrum, such as vibrational energy relaxation and spectral diffusion. Through comparison of the extracted timescales, we find that the fastest timescales are extracted over a range of excitation frequencies, where contributions from vibrational relaxation and spectral diffusion can be minimized.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
DOE OSTI2023

Development of a Reactive Force Field for Simulating Photoinitiated Acrylate Polymerization

Light-driven and photo-curable polymer based additive manufacturing (AM) has enormous potential due to its excellent resolution and precision. Acrylated radical chain-growth polymerized resins are widely used in photopolymer AM due to their fast kinetics, and often serve as a departure point for developing other resin materials for photopolymer-based AM technologies. For successful control of the photopolymer resins, the molecular basis of the acrylate free-radical polymerization has to be understood in detail. We present an optimized reactive force field (ReaxFF) for molecular dynamics (MD) simulations of acrylate polymer resins that captures radical polymerization thermodynamics and kinetics. The force field is trained against an extensive training set including density functional theory (DFT) calculations of reaction pathways along the radical polymerization from methyl acrylate to methyl butyrate, bond dissociation energies, and structures and partial charges of several molecules and radicals. We also found that it was critical to train the force field against an incorrect, nonphysical reaction pathway observed in simulations that used parameters not optimized for acrylate polymerization. As a result, the parameterization process utilizes a parallelized search algorithm, and the resulting model can describe polymer resin formation, crosslinking density, conversion rate, and residual monomers of the complex acrylate mixtures.

36 MATERIALS SCIENCE↗