Search NASA⌕ Search

SEARCH · Search NASA

Results for “molecular geometry”

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 19 records
All sourcesNASA NTRSDOE OSTI
DOE OSTI2024

On the role of asymmetric molecular geometry in high-performance organic solar cells

Although asymmetric molecular design has been widely demonstrated effective for organic photovoltaics (OPVs), the correlation between asymmetric molecular geometry and their optoelectronic properties is still unclear. To access this issue, we have designed and synthesized several symmetric-asymmetric non-fullerene acceptors (NFAs) pairs with identical physical and optoelectronic properties. Interestingly, we found that the asymmetric NFAs universally exhibited increased open-circuit voltage compared to their symmetric counterparts, due to the reduced non-radiative charge recombination. From our molecular-dynamic simulations, the asymmetric NFA naturally exhibits more diverse molecular interaction patterns at the donor (D):acceptor (A) interface as compared to the symmetric ones, as well as higher D:A interfacial charge-transfer state energy. Moreover, it is observed that the asymmetric structure can effectively suppress triplet state formation. These advantages enable a best efficiency of 18.80%, which is one of the champion results among binary OPVs. Therefore, this work unambiguously demonstrates the unique advantage of asymmetric molecular geometry, unveils the underlying mechanism, and highlights the manipulation of D:A interface as an important consideration for future molecular design.

14 SOLAR ENERGY↗
DOE OSTI2023

Importance of Dispersion in the Molecular Geometries of Mn(III) Spin-Crossover Complexes

The computational investigation of the molecular geometries of a pair of manganese(III) spin-crossover complexes is reported. For the geometry of the quintet high-spin state, density functionals significantly overestimate Mn–Namine bond distances, although the geometry for the triplet intermediate-spin state is well described. Here, comparisons with several wave function-based methods demonstrate that this error is due to the limited ability of commonly used density functionals to recover dispersion beyond a certain extent. Among the methods employed for geometry optimization, restricted open-shell Møller–Plesset perturbation theory (MP2) appropriately describes the high-spin geometry but results in a slightly shorter Mn–O distance in both spin states. On the other hand, extended multistate complete active space second-order perturbation theory (XMS-CASPT2) provides a good description of the geometry for the intermediate-spin state but also sufficiently recovers dispersion, performing well for the high-spin state. Despite the fact that the electronic structure of both spin states is dominated by one-electron configuration, XMS-CASPT2 offers a balanced approach, leading to molecular geometries with much better agreement with experiment than MP2 and DFT. A scan along the Mn–N amine bond demonstrates that for these complexes coupled cluster methods (i.e., DLPNO-CCSD(T)) also yield bond distances in agreement with experiment while multiconfiguration pair density functional theory (MC-PDFT) is unable to recover dispersion well enough, analogous to single-reference DFT.

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

CASPT2 molecular geometries of Fe( II ) spin-crossover complexes

Using fully internally contracted (FIC)-CASPT2 analytical gradients, geometry optimizations of spin-crossover complexes are reported. This approach is tested on a series of Fe(II) complexes with different sizes, ranging from 13 to 61 atoms. A combination of active space and basis set choices are employed to investigate their role in determining reliable molecular geometries. The reported strategy demonstrates that a wave function-based level of theory can be used to optimize the geometries of metal complexes in reasonable times and enables one to treat the molecular geometry and electronic structure of the complexes using the same level of theory. For a series of smaller Fe(II) SCO complexes, strong field ligands in the LS state result in geometries with the largest differences between DFT and CASPT2; however, good agreement overall is observed between DFT and CASPT2. For the larger complexes, moderate sized basis sets yield geometries that compare well with DFT and available experimental data. Here we recommend using the (10e,12o) active space since convergence to a minimum structure was more efficient than with truncated active spaces despite having similar Fe–ligand bond distances.

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

Molecular Geometry and Electronic Structure of Copper Corroles

Copper corroles are known for their unique multiconfigurational electronic structures in the ground state, which arise from the transfer of electrons from the π orbitals of the corrole to the d-orbital of copper. While density functional theory (DFT) provides reasonably good molecular geometries, the determination of the ground spin state and the associated energetics is heavily influenced by functional choice, particularly the percentage of the Hartree–Fock exchange. Using extended multireference perturbation theory methods (XMS-CASPT2), the functional choice can be assessed. The molecular geometries and electronic structures of both the unsubstituted and the meso-triphenyl copper corroles were investigated. A minimal active space was employed for structural characterization, while larger active spaces are required to examine the electronic structure. The XMS-CASPT2 investigations conclusively identify the ground electronic state as a multiconfigurational singlet (S o ) with three dominant electronic configurations in its lowest energy and characteristic saddled structure. In contrast, the planar geometry corresponds to the triplet state (T o ), which is approximately 5 kcal/mol higher in energy compared to the S o state for both the bare and substituted copper corroles. Notably, the planarity of the T o geometry is reduced in the substituted corrole compared with that in the unsubstituted one. By analyzing the potential energy surface (PES) between the S o and T o geometries using XMS-CASPT2, the multiconfigurational electronic structure is shown to transition toward a single electron configuration as the saddling angle decreases (i.e., as one approaches the planar geometry). Despite the ability of the functionals to reproduce the minimum energy structures, only the TPSSh-D 3 PES is reasonably close to the XMS-CASPT2 surface. Significant deviations along the PES are observed with other functionals.

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

Supporting Data: Inverse Trans Influence and Uranium-Arene σ-Bonding Drive Molecular Geometry: Ligand Modification from Hard to Soft Flips the Oxide from Axial to Equatorial

Supporting Computational Data from the manuscript titled "ITI and Uranium-Arene s-Bonding Drive Molecular Geometry: Ligand Modification from Hard to Soft Flips the Oxide from Axial to Equatorial" by Benedikt Kestel, Daniel Pividori, Sabyasachi Roy Chowdhury, Frank W. Heinemann, Bess Vlaisavljevich, and Karsten Meyer published in J. Am. Chem. Soc. (DOI: 10.1021/jacs.5c05029).

Chowdhury, Sabyasachi Roy [Department of Chemistry↗
DOE OSTI2021

Modulating the Molecular Geometry and Solution Self-Assembly of Amphiphilic Polypeptoid Block Copolymers by Side Chain Branching Pattern

Solution self-assembly of coil-crystalline diblock copolypeptoids has attracted increasing attention due to its capability to form hierarchical nanostructures with tailorable morphologies and functionalities. While the N-substituent (or side chain) structures are known to affect the crystallization of polypeptoids, their roles in dictating the hierarchical solution self-assembly of diblock copolypeptoids are not fully understood. Herein, we designed and synthesized two types of diblock copolypeptoids, i.e., poly(N-methylglycine)-b-poly(N-octylglycine) (PNMG-b-PNOG) and poly(N-methylglycine)-b-poly(N-2-ethyl-1-hexylglycine) (PNMG-b-PNEHG), to investigate the influence of N-substituent structure on the crystalline packing and hierarchical self-assembly of diblock copolypeptoids in methanol. With a linear aliphatic N-substituent, the PNOG blocks pack into a highly ordered crystalline structure with a board-like molecular geometry, resulting in the self-assembly of PNMG-b-PNOG molecules into a hierarchical microflower morphology composed of radially arranged nanoribbon subunits. By contrast, the PNEHG blocks bearing bulky branched aliphatic N-substituents are rod-like and prefer to stack into a columnar hexagonal liquid crystalline mesophase, which drives PNMG-b-PNEHG molecules to self-assemble into symmetrical hexagonal nanosheets in solution. A combination of time-dependent small/wide-angle X-ray scattering and microscopic imaging analysis further revealed the self-assembly mechanisms for the formation of these microflowers and hexagonal nanosheets. These results highlight the significant impact of the N-substituent architecture (i.e., linear versus branched) on the supramolecular self-assembly of diblock copolypeptoids in solution, which can serve as an effective strategy to tune the geometry and hierarchical structure of polypeptoid-based nanomaterials.

99 GENERAL AND MISCELLANEOUS↗
DOE OSTI2025

Inverse Trans Influence and Uranium-Arene σ-Bonding Drive Molecular Geometry: Ligand Modification from Hard to Soft Flips the Oxide from Axial to Equatorial

A rare example of an equatorially bound terminal uranium(V) oxo complex in a chelating sulfur-based ligand environment, namely [(mes( Me,Ad ArS) 3 )U V (O eq )(THF)] (2), is presented. Octahedrally coordinated 2 is obtained by reaction of the mesitylene-anchored tris-thiophenolate-coordinated uranium(III) complex [U III ((SAr Ad,Me ) 3 mes)] (1) with the oxygen-atom transfer reagent N 2 O. The observed, equatorially bound oxo ligand in 2 is in stark contrast to its known tris-aryloxide analog, [(mes( Me,Ad ArO) 3 )U V (O ax )(THF)] (A), where the oxo ligand occupies the typically observed axial coordination site. Complexes 1 and 2 are characterized by single-crystal X-ray diffraction analyses and spectroscopic and magnetochemical methods, including 1 H NMR, UV/vis/NIR electronic absorption, as well as EPR spectroscopy and SQUID magnetometry, thus confirming the C S symmetry and the pentavalent oxidation state of 2. Encompassing quantum chemical calculations (DFT and CASPT2) on 2 and its tris-phenolate analog A, support and rationalize the structural and electronic differences. The molecular orbital pictures show that a stabilizing σ-bonding interaction arising from the U–O eq inverse trans influence (ITI) is present in 2 but missing in A. In 2, the sulfur 3p orbitals are closer in energy to the uranium 5f manifold than the arene π-system, leading to an ITI, while U–arene σ- or δ-bonding is not observed. Although the arene orbitals remain separated from the uranium 5f orbitals in A, the absence of an ITI allows the arene a 2u orbital to engage in a σ-type interaction with the metal. Thus, incorporating a tris-thiophenolate to an arene anchor introduces a new design concept in molecular f-element chemistry. This approach stabilizes an equatorially bound U(V) oxo center, contrasting with its tris-phenolate counterpart, where oxo coordination is axial. The observed geometric divergence, driven by competing ITI and U–arene interactions, not only tunes electronic structure but also leads to differentiated reactivity: only the phenolate analogs activate H 2 O, while the thiolates do not.

Hydrocarbons↗
DOE OSTI2021

Machine-learning accelerated geometry optimization in molecular simulation

Geometry optimization is an important part of both computational materials and surface science because it is the path to finding ground state atomic structures and reaction pathways. These properties are used in the estimation of thermodynamic and kinetic properties of molecular and crystal structures. This process is slow at the quantum level of theory because it involves an iterative calculation of forces using quantum chemical codes such as density functional theory (DFT), which are computationally expensive and which limit the speed of the optimization algorithms. It would be highly advantageous to accelerate this process because then one could do either the same amount of work in less time or more work in the same time. Here, we provide a neural network (NN) ensemble based active learning method to accelerate the local geometry optimization for multiple configurations simultaneously. We illustrate the acceleration on several case studies including bare metal surfaces, surfaces with adsorbates, and nudged elastic band for two reactions. In all cases, the accelerated method requires fewer DFT calculations than the standard method. In addition, we provide an Atomic Simulation Environment (ASE)-optimizer Python package to make the usage of the NN ensemble active learning for geometry optimization easier.

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

Non-equilibrium molecular dynamics simulation to evaluate the effect of confinement on fluid flow in silica nanopores

The flow of fluids in nano-confinement has applications in separations, water purification, medical systems and in the recovery of fluids in petroleum systems. It is believed that fluids do not obey continuum laws when flowing in nanopores (or in confinement). This paper attempts to shed light on various nano confinement effects such as fluid-wall interactions, pore size and molecular geometry, using molecular dynamic simulation. Here, water, hexane, and methanol flow behaviors were simulated for pore diameters ranging from 1 to 8 nm. In addition to the density analysis of the confined fluids, the fluid flows in saturated nanopores were simulated by the sectional flow method. Water and methanol molecules were completely stabilized in the 1 nm pore. Hexane molecules could not enter the 1 nm pore, due to geometric considerations. For the 2 – 8 nm pores, all fluids showed reduced flow rate compared to the Hagen–Poiseuille flow, due to an interfacial molecular layer stabilized at the pore surface. Flow reduction observed in these studies is contrary to significant flow enhancements observed for fluid flow in carbon nanotubes of similar dimensions. Water flow showed almost constant stick length for all the pore sizes. Methanol flow had the largest stick length. Hexane flow was reduced because of overcrowding of molecules at the pore surface. The effect of confinement diminished with an increase in pore diameter.

42 ENGINEERING↗
DOE OSTI2023

Nuclear–Electronic Orbital QM/MM Approach: Geometry Optimizations and Molecular Dynamics

Hybrid quantum mechanical/molecular mechanical (QM/MM) methods allow simulations of chemical reactions in atomistic solvent and heterogeneous environments such as proteins. Herein, the nuclear–electronic orbital (NEO) QM/MM approach is introduced to enable the quantization of specified nuclei, typically protons, in the QM region using a method such as NEO-density functional theory (NEO-DFT). This approach includes proton delocalization, polarization, anharmonicity, and zero-point energy in geometry optimizations and dynamics. Expressions for the energies and analytical gradients associated with the NEO-QM/MM method, as well as the previously developed polarizable continuum model (NEO-PCM), are provided. Geometry optimizations of small organic molecules hydrogen bonded to water in either dielectric continuum solvent or explicit atomistic solvent illustrate that aqueous solvation can strengthen hydrogen-bonding interactions for the systems studied, as indicated by shorter intermolecular distances at the hydrogen-bond interface. We then performed a real-time direct dynamics simulation of a phenol molecule in explicit water using the NEO-QM/MM method. Furthermore, these developments and initial examples provide the foundation for future studies of nuclear–electronic quantum dynamics in complex chemical and biological environments.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
NASA NTRS1978

H-H, C-H, and C-C NMR spin-spin coupling constants calculated by the FP-INDO method for aromatic hydrocarbons

The FP-INDO (finite perturbation-intermediate neglect of differential overlap) method is used to calculate the H-H, C-H, and C-C coupling constants in hertz for molecules of six different benzenoid hydrocarbons: benzene, naphthalene, biphenyl, anthracene, phenanthrene, and pyrene. The calculations are based on both the actual and the average molecular geometries. It is found that only the actual molecular geometries can always yield the correct relative order of values for the H-H coupling constants. For the calculated C-C coupling constants, as for the calculated C-H coupling constants, the signs are positive (negative) for an odd (even) number of bonds connecting the two nuclei. Agreements between the calculated and experimental values of the coupling constants for all six molecules are comparable to those reported previously for other molecules.

Long, S. A. T.↗
NASA NTRS1987

Structure, vibrational spectrum, and energetics of the CH5(+) ion - A theoretical investigation

Analytic derivative techniques are utilized to determine molecular geometries, vibration spectra, molecular force field, proton affinity, and energetics of clustering of protonated methane. It is observed that the geometries are gradient optimized, while the harmonic force fields and IR intensities are also determined analytically at the SCF level. It is determined that the frequency of the torsional model for rotation of the H2 about the pseudothreefold axis of the CH3(+) group is extremely low, and the proton affinity of CH4 is calculated as 129.0 kcal/mol.

Komornicki, Andrew↗
DOE OSTI2023

Controlling Molecular Structure and Spin with Multiconfigurational Quantum Chemistry (Final Technical Report)

For many first-row transition metal complexes, structure-property relationships can be obtained from high-level molecular geometry optimizations and subsequent electronic structure studies. The project funded under this award utilized newly implemented fully internally contracted (FIC) nuclear gradients for extended multi-state (XMS) complete active space second-order multireference perturbation theory (CASPT2) to explore geometry changes in first-row transition metal coordination complexes. At the start of the project, only one full geometry optimization using FIC-CASPT2 analytical gradients had been reported (J. Chem. Theory Comput., 2016, 12 (8), 3781), and many open-questions regarding the performance and achievable accuracy in applying such computations to larger complexes persisted. Key deliverables in the report include the implementation of the numerical Hessian and subsequent vibrational analysis in the BAGEL program package. This includes both full Hessian vibrational analysis (FHVA) and a partial Hessian vibrational analysis (PHVA). The latter of which allows us to reduce the significant cost of computing the vibrational frequencies in a large molecule by focusing on the modes of highest interest. The gradient code also proved useful in the development of an approach to improve the Hubbard-U correction by removing bias in predicting spin-splitting in Fe(II) complexes via plane-wave density functional theory (DFT). Finally, we showed for three families of complexes (spin-crossover (SCO) complexes, metallocorroles complexes, and those with metal-metal bonds) that CASPT2 can result in good molecular geometries and established best practices in undertaking this work. We are now using this approach in collaborative projects with experimental groups, which would not have been possible at the start of this project. A common theme also arose through this work that static and dynamic correlation must be recovered in a balanced way to yield quantitative results. By systematically varying how both effects were included, we were able to make key chemical insights. This work supported the training of two postdoctoral scholars, one graduate student, and one undergraduate student in applying high-level wavefunction based methods to challenging systems.

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

Augmenting Molecular Graphs with Geometries via Machine Learning Interatomic Potentials

Accurate molecular property predictions require 3D geometries, which are typically obtained using expensive methods such as density functional theory (DFT). Here, we attempt to obtain molecular geometries by relying solely on machine learning interatomic potential (MLIP) models. To this end, we first curate a large-scale molecular relaxation dataset comprising 3.5 million molecules and 300 million snapshots. Then MLIP pre-trained models are trained with supervised learning to predict energy and forces given 3D molecular structures. Once trained, we show that the pre-trained models can be used in different ways to obtain geometries either explicitly or implicitly. First, it can be used to obtain approximate low-energy 3D geometries via geometry optimization. While these geometries do not consistently reach DFT-level chemical accuracy or convergence, they can still improve downstream performance compared to non-relaxed structures. To mitigate potential biases and enhance downstream predictions, we introduce geometry fine-tuning based on the relaxed 3D geometries. Second, the pre-trained models can be directly fine-tuned for property prediction when ground truth 3D geometries are available. Our results demonstrate that MLIP pre-trained models trained on relaxation data can learn transferable molecular representations to improve downstream molecular property prediction and can provide practically valuable but approximate molecular geometries that benefit property predictions. Our code is publicly available at: https://github.com/divelab/AIRS/.

Fu, Cong [Texas A & M Univ., College Station, TX (↗
DOE OSTI2021

Quadrupole-mediated dielectric response and the charge-asymmetric solvation of ions in water

Treating water as a linearly responding dielectric continuum on molecular length scales allows very simple estimates of the solvation structure and thermodynamics for charged and polar solutes. While this approach can successfully account for basic length and energy scales of ion solvation, computer simulations indicate not only its quantitative inaccuracies but also its inability to capture some basic and important aspects of microscopic polarization response. Here, we consider one such shortcoming, a failure to distinguish the solvation thermodynamics of cations from that of otherwise-identical anions, and we pursue a simple, physically inspired modification of the dielectric continuum model to address it. The adaptation is motivated by analyzing the orientational response of an isolated water molecule whose dipole is rigidly constrained. Its free energy suggests a Hamiltonian for dipole fluctuations that accounts implicitly for the influence of higher-order multipole moments while respecting constraints of molecular geometry. Finally, we propose a field theory with the suggested form, whose nonlinear response breaks the charge symmetry of ion solvation. An approximate variational solution of this theory, with a single adjustable parameter, yields solvation free energies that agree closely with simulation results over a considerable range of solute size and charge.

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

Geometry optimization speedup through a geodesic approach to internal coordinates

We present a new geodesic-based method for geometry optimization in a basis set of redundant internal coordinates. Overall, our method updates the molecular geometry by following the geodesic generated by a displacement vector on the internal coordinate manifold, which dramatically reduces the number of steps required to converge to a minimum. Our method can be implemented in any existing optimization code, requiring only implementation of derivatives of the Wilson B-matrix and the ability to numerically solve an ordinary differential equation.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗