Search NASA⌕ Search

SEARCH · Search NASA

Results for “molecular simulation”

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 343 records · Page 19
All sourcesNASA NTRSDOE OSTI
DOE OSTI2025

Fouling behavior of zwitterionic membranes compared to polyamide membranes

Membrane fouling remains a critical bottleneck for reverse osmosis (RO) desalination, driving energy consumption and reducing membrane lifetime. Here, we employ all-atom molecular dynamics simulations to investigate the antifouling behavior of random zwitterionic amphiphilic copolymer (r-ZAC) membranes composed of sulfobetaine methacrylate (SBMA) and allyl methacrylate (AMA), benchmarked against conventional polyamide (PA) RO membranes. Structural and dynamical analyses—including radial distribution functions, coordination numbers, tetrahedral order parameters, vector orientation, and residence-time correlation functions—reveal that r-ZAC surfaces sustain tightly bound, long-lived hydration layers with preserved tetrahedrality and anisotropic water orientation, in sharp contrast to the weak and disordered hydration of PA. Steered molecular dynamics simulations demonstrate that r-ZAC membranes impose substantial free-energy barriers to foulant approach (alginate ≈ 90 kcal/mol, sucrose ≈ 35 kcal/mol, humic acid ≈ 15 kcal/mol), whereas PA membranes exhibit negligible barriers (< 1 kcal/mol) and thermodynamically favorable adsorption. Detailed foulant–surface interaction analyses show that zwitterionic hydration and electrostatic heterogeneity in r-ZAC suppress adhesion, except in the case of amphiphilic humic acid, which exploits multiple binding modes. Together, these results establish molecular-level design principles for antifouling membranes: the combination of zwitterionic hydration, structured interfacial water, and controlled amphiphilic balance in r-ZAC membranes provides superior resistance to organic fouling relative to PA.

Cross-linked polyamide↗
DOE OSTI2025

Accurate and efficient parameterization of an atomic cluster expansion (ACE) potential for ammonia under extreme conditions

We present a machine learning interatomic potential for ammonia designed to capture its complex multiphase behavior, including both molecular and superionic phases. The potential is based on the atomic cluster expansion (ACE) formulation and has been parameterized to facilitate high-fidelity molecular dynamics simulations of ammonia under extreme conditions, for pressures up to 100 GPa and for temperatures above 500 K and up to 6000 K. A diverse range of configurations was generated through high-quality ab initio molecular dynamics simulations, covering insulating and superionic ice phases, liquid ammonia, molecular nitrogen (N 2 ) and hydrogen (H 2 ), and metastable compounds that form upon dissociation, including $NH^{+}_{4}$, $H^{+}_{3}$, N 2 H 4 , and N 3 H. We demonstrate that the ammonia ACE potential accurately reproduces experimental and density functional theory predicted isotherms and Hugoniots. Crucially, the potential is able to capture the intricate phase behavior of ammonia, including the transition from insulating molecular fluid to the superionic phase. This work provides a robust interatomic potential that can be used for large-scale, accurate simulations of ammonia under extreme thermodynamic conditions, offering a powerful tool for investigating its behavior in various phases and applications.

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

Ranges of applicability for the continuum beam model in the mechanics of carbon nanotubes and nanorods

Limitations in the validity of the continuum beam model for carbon nanotubes (NTs) and nanorods are examined. Applicability of all assumptions used in the model is restricted by the two criteria for geometric parameters that characterize the structure of NTs. The key non-dimensional parameters that control the NT buckling behavior are derived via dimensional analysis of the nanomechanical problem. A mechanical law of geometric similitude for NT buckling is extended from continuum mechanics for different molecular structures. A model applicability map, where two classes of beam-like NTs are identified, is constructed for distinct ranges of non-dimensional parameters. Expressions for the critical buckling loads and strains are tailored for two classes of NTs and compared with the data provided by the molecular dynamics simulations. copyright 2001 Elsevier Science Ltd. All rights reserved.

Continuum beam model↗
DOE OSTI2024

Melting point of iron at high pressure: An assessment of uncertainties and effect of electronic temperature

An accurate calculation of the melting point of iron at various pressures in the Earth's core is important for understanding the core structure, geodynamo, and the Earth's history. Previous studies have assessed the melt line of iron at these extreme conditions using various experimental measurement techniques as well as both ab initio and classic molecular dynamics simulations. However, experimental measurements have uncertainties up to several hundred Kelvin, and inconsistencies remain among simulation results. Here in this work, we propose an iterative framework that couples density functional theory (DFT) calculations and molecular dynamics simulations performed using an ensemble of interatomic potentials to assess the effect of electronic temperature on the melting point. We systematically validate the potentials by comparing lattice constants and phonon dispersion curves at 0 K and enthalpy differences between liquid and HCP, FCC, BCC phases of iron close to the melt line at 300 GPa with DFT. Our results show that HCP iron melts at 6144 K (at 300 GPa), BCC phase is thermodynamically unstable, and FCC is metastable at this temperature. The melting points of FCC and BCC phases at 300 GPa are 5858 and 5647 K, respectively.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
NASA NTRS1999

Bone Blood Flow During Simulated Microgravity: Physiological and Molecular Mechanisms

Blood flow to bone has been shown to affect bone mass and presumably bone strength. Preliminary data indicate that blood flow to the rat femur decreases after 14 days of simulated microgravity, using hindlimb suspension (HLS). If adult rats subjected to HLS are given dobutamine, a synthetic catecholamine which can cause peripheral vasodilation and increased blood flow, the loss of cortical bone area usually observed is prevented. Further, mechanisms exist at the molecular level to link changes in bone blood flow to changes in bone cell activity, particularly for vasoactive agents like nitric oxide (NO). The decreases in fluid shear stress created by fluid flow associated with the shifts of plasma volume during microgravity may result in alterations in expression of vasoactive agents such as NO, producing important functional effects on bone cells. The primary aim of this project is to characterize changes in 1) bone blood flow, 2) indices of bone mass, geometry, and strength, and 3) changes in gene expression for modulators of nitric oxide activity (e.g., nitric oxide synthase) and other candidate genes involved in signal transduction of mechanical loading after 3, 7, 14, 21, and 28 days of HLS in the adult rat. Using a rat of at least 5 months of age avoids inadvertently studying effects of simulated microgravity on growing, rather than adult, bone. Utilizing the results of these studies, we will then define how altered blood flow contributes to changes in bone with simulated microgravity by administering a vasodilatory agent (which increases blood flow to tissues) during hindlimb suspension. In all studies, responses in the unloaded hindlimb bones (tibial shaft, femoral neck) will be compared with those in the weightbearing humeral shaft and the non-weightbearing calvarium (skull) from the same animal. Bone volumetric mineral density and geometry will be quantified by peripheral quantitative CT; structural and material properties of the long bones will be determined by 3-point bending (tibia, humerus) or compression (femoral neck) testing to failure. A unique aspect of these studies will be defining the time course of changes in gene expression in bone cell populations with unloading, accomplished with Northern blots, in situ hybridization, and immunohistochemistry. These studies have high relevance for concurrent protocols being proposed by investigators on NSBRI Cardiovascular and Muscle teams, with blood flow data available on a number of tissues other than bone. Further, dobutamine and other Beta-agonists have been tested as countermeasures for altered muscle and cardiovascular function. Results of the intervention tested in our studies have potential relevance for a number of systemic changes seen with prolonged spaceflight.

Bloomfield, Susan A.↗
DOE OSTI2026

Custom-trained Machine-learning Interatomic Potentials: ZnCl2 Aqueous Solution

This dataset was generated using an iterative active-learning strategy implemented in the ArcaNN software package (https://github.com/arcann-chem/arcann_training) to train machine-learning interatomic potentials for aqueous ZnCl2 solutions. Each active-learning cycle consisted of three stages: training, exploration, and labeling. The initial training set combined configurations generated in this work from enhanced-sampling ab initio molecular dynamics simulations with configurations from a previously reported neural-network-potential study of aqueous ZnCl2. The enhanced-sampling ab initio molecular dynamics simulations involved Zn–Cl separation and the chloride coordination number around Zn²? as collective variables. These configurations served as the seed dataset. Subsequent active-learning cycles expanded the training set by identifying and labeling configurations that were poorly represented by the current models, thereby improving coverage of ion-association states and changes in local coordination and charge-state environments relevant to the solution free-energy landscape. For all selected configurations, single-point calculations of the total energies and atomic forces were performed within density functional theory using the CP2K Quickstep module. Reference calculations employed the revPBE-D3 and r2SCAN exchange-correlation functionals. Motivated by recent work on aqueous Zn²?, the main revPBE calculations omitted D3 dispersion contributions involving Zn²?, while retaining the D3 correction for water and chloride. For comparison, fully dispersion-corrected revPBE-D3 reference calculations were also performed, with D3 applied to all species, including Zn²?. Valence electrons were treated explicitly, while core electrons were represented using norm-conserving Goedecker–Teter–Hutter pseudopotentials. The wave functions were expanded using the mixed Gaussian-and-plane-wave scheme with TZV2P-MOLOPT basis sets for all elements and a 600 Ry auxiliary plane-wave cutoff for the electron density. Self-consistent-field convergence was accelerated using the orbital-transformation and Direct Inversion in the Iterative Subspace algorithms, with a convergence threshold of 10?6. All single-point calculations were performed in periodic orthorhombic cells. The CELL_REF keyword in CP2K was used to define a fixed reference cell with a box length of 25 Å. This treatment ensured a consistent reference for configurations extracted from NpT trajectories with fluctuating cell dimensions. The resulting DFT energies and atomic forces constitute the ground-truth labels used to train the MLIPs. The resulting MLIP was trained for aqueous ZnCl2 solutions spanning concentrations from 0 to 30 molal and a broad pH range, from strongly acidic to strongly basic conditions. Representative examples of configurations included in the MLIP training dataset are provided below. These include 1) Representative configurations from the dataset labeled at the revPBE-D3 level, with D3 dispersion interactions involving Zn2+ excluded (revPBE-wo-D3). 2) Representative configurations from the dataset labeled at the fully dispersion-corrected revPBE-D3 level, with D3 interactions applied to all species, including Zn2+ (revPBE-D3). 3) Representative configurations from the dataset labeled at the r2SCAN level of theory (r2SCAN).

Dinpajooh, Mohammadhasan [Pacific Northwest Nation↗
DOE OSTI2024

Beyond nothingness in the formation and functional relevance of voids in polymer films

Abstract Voids—the nothingness—broadly exist within nanomaterials and impact properties ranging from catalysis to mechanical response. However, understanding nanovoids is challenging due to lack of imaging methods with the needed penetration depth and spatial resolution. Here, we integrate electron tomography, morphometry, graph theory and coarse-grained molecular dynamics simulation to study the formation of interconnected nanovoids in polymer films and their impacts on permeance and nanomechanical behaviour. Using polyamide membranes for molecular separation as a representative system, three-dimensional electron tomography at nanometre resolution reveals nanovoid formation from coalescence of oligomers, supported by coarse-grained molecular dynamics simulations. Void analysis provides otherwise inaccessible inputs for accurate fittings of methanol permeance for polyamide membranes. Three-dimensional structural graphs accounting for the tortuous nanovoids within, measure higher apparent moduli with polyamide membranes of higher graph rigidity. Our study elucidates the significance of nanovoids beyond the nothingness, impacting the synthesis‒morphology‒function relationships of complex nanomaterials.

36 MATERIALS SCIENCE↗
DOE OSTI2026

Ionic Liquid-Enhanced Interfaces to Boost Reactive C O2 Capture

The addition of ionic liquids (ILs) to a mixture containing a molecular solvent and other ionic species can induce the heterogeneous redistribution of cations and anions at the gas–liquid interface. This nonuniform redistribution of cations and anions driven by the differences in the solvophilicity of ions can improve the thermophysical and interfacial properties of such mixtures, creating a local chemical environment that is conducive to some reactions. In this work, ILs are added to a mixture of potassium hydroxide (KOH) and ethylene glycol (EG), used as a reactive absorbent and electrolyte in the migration-assisted moisture-gradient (MAMG) process for CO 2 capture. Molecular dynamics (MD) simulations are employed to probe into the effects of complex ion–ion and ion–solvent interactions and to examine the chemical composition at the gas–liquid interface. A total of 12 systems are investigated using molecular simulations to identify trends in the performance of IL additives based on the choice of cation, anion, and IL concentration. The cation effects are studied using IL additives based on 1-ethyl-3-methylimidazolium ([EMIM] + ) and 1-butyl-3-methylimidazolium ([BMIM] + ), while the impact of anions is examined using additives based on dicyanamide [DCA] − , triflate [TfO] − , bistriflimide [NTf 2 ] − , and hexafluorophosphate [PF 6 ] − anions, respectively. The influence of the IL concentration is also evaluated at molar concentrations between 1% and 4%. The simulation results indicate that the use of IL additives can affect the physical CO 2 solubility, surface tension, and the localization of CO 2 around the [OH] − ions at the gas–liquid interface. It is also evident that the choice of cations, anions, and IL concentration determines the extent to which the IL additives impact the local physicochemical properties. Physical dissolution, diffusive transport, and interaction with [OH] − are critical intermediate steps toward reactive CO 2 capture using a liquid absorbent. Hence, the improvement in one or more of these properties, aided by IL additives, is expected to improve the overall CO 2 capture performance. Experiments reaffirmed the impact of IL additives on CO 2 capture performance and the sensitivity to the choice of the cation, anion, and concentration of the IL additive.

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

Release dynamics of nanodiamonds created by laser-driven shock-compression of polyethylene terephthalate

Laser-driven dynamic compression experiments of plastic materials have found surprisingly fast formation of nanodiamonds (ND) via X-ray probing. This mechanism is relevant for planetary models, but could also open efficient synthesis routes for tailored NDs. We investigate the release mechanics of compressed NDs by molecular dynamics simulation of the isotropic expansion of finite size diamond from different P-T states. Analysing the structural integrity along different release paths via molecular dynamic simulations, we found substantial disintegration rates upon shock release, increasing with the on-Hugnoiot shock temperature. We also find that recrystallization can occur after the expansion and hence during the release, depending on subsequent cooling mechanisms. Our study suggests higher ND recovery rates from off-Hugoniot states, e.g., via double-shocks, due to faster cooling. Laser-driven shock compression experiments of polyethylene terephthalate (PET) samples with in situ X-ray probing at the simulated conditions found diamond signal that persists up to 11 ns after breakout. In the diffraction pattern, we observed peak shifts, which we attribute to thermal expansion of the NDs and thus a total release of pressure, which indicates the stability of the released NDs.

36 MATERIALS SCIENCE↗
DOE OSTI2025

Structure and Dynamics of CO 2 at the Air–Water Interface from Classical and Neural Network Potentials

The accurate description of the structure and dynamics of CO 2 at the instantaneous air–water interface, along with the effects of surface fluctuations on the CO 2 -transport processes, is essential for the development of negative emission technologies aimed at minimizing climate change. In this study, we performed molecular dynamics simulations of CO 2 at the air–water interface using neural network potentials (NNPs) trained on ab initio data generated through density-functional-theory-based molecular dynamics simulations. We compared these results with classical force fields to assess their performance in modeling interfacial CO 2 behavior. Our findings revealed that the asymmetric interactions, coupled with thermal surface fluctuations at the air–water interface, significantly influence CO 2 transport into the aqueous phase. The simulations demonstrate that classical force fields underestimate both the free energy of CO 2 transport and the strength of its interactions at the interface compared with the neural network potentials. In conclusion, the free energy and the interfacial dynamics of CO 2 are primarily influenced by the distribution of water within the instantaneous interfacial water layer, responsible for creating an asymmetric intermolecular interaction environment within the interfacial region.

Ab initio molecular dynamics↗
DOE OSTI2025

Long timescale solvation dynamics and confinement: The case of non-ionic deep eutectic solvents of lauric acid and N-methylacetamide

Microscopic segregation and molecular heterogeneities in complex liquids are the result of the interplay between different intermolecular forces, all of which contribute to the energy landscape of the system. A consequence of the intricate energy landscape is the nontrivial effect on the solvation dynamics. Here, the impact of molecular heterogeneities on the solvation dynamics is studied using infrared spectroscopies and molecular dynamics simulations. In particular, this study focuses on the dynamical effect of nanoscopic heterogeneities present in deep eutectic solvents (DESs) composed of lauric acid (LA) and N-methylacetamide (NMA). To this end, a molecular probe containing a carbon triple bond is used as an infrared reporter. The results show that the vibrational probe is likely to be located in the NMA polar domains. Furthermore, the probe solvation dynamics derived from the 2DIR spectra presents a slowdown of its timescale with increasing LA concentration in the DES. Kubo modeling of the probe solvation dynamics shows a correlation between the amplitude of its long time component and the presence of molecular heterogeneities in the sample. Semiclassical modeling of the vibrational line shape of the triple bond stretch demonstrates that the heterogeneities affect the whole solvation dynamics of the system through the amplitudes of the frequency fluctuations. Molecular dynamics simulations confirm the experimental results and their interpretation by showing a slowdown of the solvation dynamics when the LA heterogeneities are present. Overall, the study presents a molecular framework to explain the effect of confinement created by nanoscopic LA heterogeneities on the solvation dynamics of the system.

Chemical compounds↗
DOE OSTI2024

Anisotropic interactions for continuum modeling of protein–membrane systems

In this work, a model for anisotropic interactions between proteins and cellular membranes is proposed for large-scale continuum simulations. The framework of the model is based on dynamic density functional theory, which provides a formalism to describe the lipid densities within the membrane as continuum fields while still maintaining the fidelity of the underlying molecular interactions. Within this framework, we extend recent results to include the anisotropic effects of protein–lipid interactions. As applications, we consider two membrane proteins of biological interest: a RAS–RAF complex tethered to the membrane and a membrane embedded G protein-coupled receptor. A strong qualitative and quantitative agreement is found between the numerical results and the corresponding molecular dynamics simulations. Combining the scope of continuum level simulations with the details from molecular level particle simulations enables research into protein–membrane behaviors at a more biologically relevant scale, which crucially can also be accessed via experiment.

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

Calculation of Non-Bonded Forces Due to Sliding of Bundled Carbon Nanotubes

An important consideration for load transfer in bundles of single-walled carbon nanotubes is the nonbonded (van der Waals) forces between the nanotubes and their effect on axial sliding of the nanotubes relative to each other. In this research, the non-bonded forces in a bundle of seven hexagonally packed (10,10) single-walled carbon nanotubes are represented as an axial force applied to the central nanotube. A simple model, based on momentum balance, is developed to describe the velocity response of the central nanotube to the applied force. The model is verified by comparing its velocity predictions with molecular dynamics simulations that were performed on the bundle with different force histories applied to the central nanotube. The model was found to quantitatively predict the nanotube velocities obtained from the molecular dynamics simulations. Both the model and the simulations predict a threshold force at which the nanotube releases from the bundle. This force converts to a shear yield strength of 10.5-11.0 MPa for (10,10) nanotubes in a bundle.

Frankland, S. J. V.↗
DOE OSTI2024

Atomistic simulations to reveal HIP-bonding mechanisms of Al6061/Al6061

Molecular dynamics simulations were employed to understand the diffusion bonding process during hot isostatic pressing (HIP) of Al6061/Al6061 alloy. Simulations of the HIP process reveal atomistic phenomena that are difficult or unlikely to be observed experimentally and provide useful insights into the mechanism of diffusion and bonding. Here, the results reveal that at the start of the HIP process, a massive incursion of oxygen atoms occurs from the pre-existing γ-Al 2 O 3 to the 6061 region across the interphase interface. These oxygen atoms interact with the enriched Mg atom layer present at the existing γ-Al 2 O 3 and 6061 matrix to form a secondary complex Mg 2 Al 2 O 5 phase. Diffusion calculations also show that transport of atoms due to the applied pressure is 4–5 orders of magnitude higher than would occur in the absence of HIP conditions. The Mg 2 Al 2 O 5 phase also provides efficient pathways for the rapid transport of Mg atoms. Because of the higher diffusion coefficients observed for Mg within the phase, Mg atoms can move more swiftly compared to their diffusion within other phases such as γ-Al 2 O 3 . This accelerated mobility facilitates the rapid movement of Mg atoms across the interface, leading to changes in the local composition and the potential growth of the Mg 2 Al 2 O 5 phase.

36 MATERIALS SCIENCE↗
DOE OSTI2025

Investigating Morphology and Diffusion in Simulations of Precise Anion-Conducting Polymers

Using atomistic molecular dynamics simulations, we investigate the morphology and transport properties of a new class of polymers which are functionalized with quaternary ammonium groups for use as anion exchange membranes. The polymers are precision polyolefins with either a trimethylammonium (p5CNMe3) or a dimethyl-hexyl ammonium (p5CNMe2Hx) pendant group at every fifth carbon along a polyethylene backbone. Simulations are performed at hydration levels of 5, 10, 15, and 20 water molecules per ammonium group. The hydrated polymers form nanoscale, percolated hydrophilic domains (water channels) in the hydrophobic polymer matrix that become wider with increasing water content. Water and hydroxide anion diffusion coefficients also increase with increasing water content. The morphology of the water domains is similar in both polymers, while the diffusion coefficients are somewhat lower in p5CNMe2Hx at fixed water content. Furthermore, the diffusion coefficients in both polymers fall on the same curve as a function of the fractal dimension of the percolated water channels, which appears to be a useful scalar measure of the effects of the nanoscale morphology on water and hydroxide anion transport.

Anion exchange membrane↗
NASA NTRS2003

[Results of Simulation Studies]

Lattice Monte Carlo and off-lattice molecular dynamics simulations of h(sub 1)t(sub 4) and h(sub 4)t(sub l) (head/tail) amphiphile solutions have been performed as a function of surfactant concentration and temperature. The lattice and off-lattice systems exhibit quite different self-assembly behavior at equivalent thermodynamic conditions. We found that in the weakly aggregating regime (no preferred-size micelles), all models yield similar micelle size distributions at the same average aggregation number, albeit at different thermodynamic conditions (temperatures). In the strongly aggregating regime, this mapping between models (through temperature adjustment) fails, and the models exhibit qualitatively different micellization behavior. Incipient micellization in a model self-associating telechelic polymer solution results in a network with a transient elastic response that decays by a two-step relaxation: the first is due to a heterogeneous jump-diffusion process involving entrapment of end-groups within well-defined clusters and this is followed by rapid diffusion to neighboring clusters and a decay (terminal relaxation) due to cluster disintegration. The viscoelastic response of the solution manifests characteristics of a glass transition and entangled polymer network.

Source record↗
NASA NTRSDate not supplied

Efficient Implementation for Unitary Coupled Cluster State Preparation for Near-Term Quantum Computers

Unitary coupled cluster theory (UCC) is a common wave function ansatz for quantum simulation of molecular electronic structure using the variational quantum eigenvalue solver (VQE). Even for small molecules using a double-ζ basis, the number of variational parameters required to minimize the electronic energy (i.e., optimize the circuit) is large and beyond the reach of current quantum computers. For example, a circuit simulating C2 using the UCCSD ansatz and the cc-pVDZ basis set with frozen-core will require over 10,000 variational parameters and a Hilbert space of over 10^8 determinants. To make progress on simulating such molecular systems on near-term quantum computers, we explore how much of the optimization can be approximately prepared with classical simulation while reducing the number of optimization steps performed on a quantum device. Recently, Chen, Cheng, and Freericks [J. Chem. Theory Comput. 2021, 17, 841-847] presented an algorithm for the factorized form of the UCC ansatz that allows for efficient UCC optimizations on classical hardware. We flip the algorithm around and use it to prepare approximate quantum circuits for systems that require a large number of qubits to represent. We will present results from our implementation and discuss strategies for incorporating this implementation for algorithms involving near-term quantum computers.

J Wayne Mullinax↗
NASA NTRSDate not supplied

Efficient Implementation for Unitary Coupled Cluster State Preparation for Near-Term Quantum Computers

Unitary coupled cluster theory (UCC) is a common wave function ansatz for quantum simulation of molecular electronic structure using the variational quantum eigenvalue solver (VQE). Even for small molecules using a double-ζ basis, the number of variational parameters required to minimize the electronic energy (i.e., optimize the circuit) is large and beyond the reach of current quantum computers. For example, a circuit simulating C2 using the UCCSD ansatz and the cc-pVDZ basis set with frozen-core will require over 10,000 variational parameters and a Hilbert space of over 10^(8) determinants. To make progress on simulating such molecular systems on near-term quantum computers, we explore how much of the optimization can be approximately prepared with classical simulation while reducing the number of optimization steps performed on a quantum device. Recently, Chen, Cheng, and Freericks [J. Chem. Theory Comput. 2021, 17, 841-847] presented an algorithm for the factorized form of the UCC ansatz that allows for efficient UCC optimizations on classical hardware. We flip the algorithm around and use it to prepare approximate quantum circuits for systems that require a large number of qubits to represent. We will present results from our implementation and discuss strategies for incorporating this implementation for algorithms involving near-term quantum computers.

Quantum Computing↗