Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Structured liquids”

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 379 records · Page 21
All sourcesNASA NTRSDOE OSTI
DOE OSTI2024

Retrofitting Holcim Ste. Genevieve Cement Plant with CO2 Capture Plant Using Air Liquide Cryocap™ FG Technology

The global cement manufacturing industry is a major contributor to carbon dioxide emissions. The International Energy Agency's "Net Zero Emissions by 2050 Scenario" identifies CCS as a major strategy for meeting that goal. This project is among the first attempts to transfer capture technology developed at coal-fired power plants to the cement industry. The main objective of the project is to execute and complete a front-end engineering and design (FEED) studies for commercial-scale, carbon capture projects that separates 95% of the total CO2 emissions at the Holcim (US) Ste. Genevieve cement manufacturing facility using Air Liquide’s Pressure Swing Adsorption system (PSA) assisted Cryocap™ technology. The Holcim Ste. Genevieve cement plant in Missouri, US, boasts one of the largest single cement production lines in the world, with a capacity of approximately 12,000 t/day. The plant currently uses traditional fuels, namely coal and petcoke. The captured CO2 will be pipeline and geological storage grade. The industrial host site emits approximately 3.0 million tonne CO2/yr. Air Liquide’s Cryocap™ technology has been developed over the last 18+ years for CO2 capture applications. It has been shown to be applicable to a variety of industrial applications (e.g., steel, cement, SMR, Fluidized Catalytic Crackers (FCCs)). Cryocap™ FG consists of a Pressure Swing Adsorption (PSA) unit coupled with a Cryogenic System. The PSA pre-concentrates the CO2 from the flue gas, while the cryogenic unit enables the CO2 purity to be increased to the desired level. The project team is led by the Prairie Research Institute at the University of Illinois at Urbana-Champaign. The tasks include: complete FEED study for retrofitting the industrial facility with a carbon capture system to support developing a detailed cost estimate; business case analysis outlining the anticipated revenue and credits if projects was built and operated; technoeconomic analysis (TEA) outlining how capture system achieves DOE capture goals; and life cycle (LCA) analysis demonstrating zero net carbon emissions. The FEED study was successfully completed. This includes completing the process basis of design; preliminary engineering; outside battery limits (OSBL) detailed engineering including a Zero Liquid Discharge (ZLD) wastewater treatment system; inside battery limits (ISBL) detailed engineering [1]. An overall project capital cost estimate within a -20%/+30% accuracy was developed. The major contributors to the Total Plant Cost (TPC), by system, are the costs associated with the Outside Battery Limit (OSBL) section of the plant which includes a new river water intake structure and a Zero Liquid Discharge (ZLD) system. By cost category, the major contributors to the TPC are equipment and subcontractor costs, followed closely by engineering, construction management, home office and contractor fees. The TEA has been created to reflect the findings of the project. It analyzes the economic performance of the Cryocap™ technology by reviewing the estimated capital costs, operating cost, and revenue. The Cost of Capture (COC) associated with the Cryocap™ technology for 95% CO2 capture, when considering NETL 2018 economic assumptions (42/58 debt/equity ratio, 5.15% interest on debt and 1.42% return on equity in real dollars) and 2022 economic assumptions (42/58 debt/equity ratio, 8.82% interest on debt and 4.90% return on equity in real dollars) was found to be much lower than that for the DOE-NETL’s base-line cases. The highest contributors to the COC are annualized capital expenditures (CAPEX) and electricity consumption which can be offset by using lower cost renewable sources. The LCA was conducted using OpenLCA which is an open-source software that is recommended by NETL. The database utilized for this study was a modified version of TRACI 2.1 (developed by the US. Environmental Protection Agency’s National Risk Management Research Laboratory and modified by NETL). The Cryocap™ FG technology does not consume fuels in significant quantities and does not utilize specialized chemical solvents subject to decomposition, such as those utilized in amine-based carbon capture systems. The Cryocap™ FG technology mainly utilizes electricity as its energy input; hence, its calculated emissions are mainly associated with the generation of electricity offsite and are dependent on the energy matrix of the grid at the time of project implementation. The water consumption impact of the Cryocap™ FG is mostly for makeup of the water lost by evaporation in the cooling tower; however, the carbon capture plant will be equipped with a ZLD system to avoid effluent streams and minimize water consumption. The successful construction and operation of this plant based on this study results will provide a means to demonstrate an economically attractive and transformational capture technology that can be used to retrofit existing plants and be deployed at new plants.

01 COAL, LIGNITE, AND PEAT↗
DOE OSTI2024

Attosecond-pump attosecond-probe x-ray spectroscopy of liquid water

Attosecond-pump/attosecond-probe experiments have long been sought as the most straightforward method for observing electron dynamics in real time. Although there has been much success with overlapped near-infrared femtosecond and extreme ultraviolet attosecond pulses combined with theory, true attosecond-pump/attosecond-probe experiments have been limited. Here we used a synchronized attosecond x-ray pulse pair from an x-ray free-electron laser to study the electronic response to valence ionization in liquid water through all x-ray attosecond transient absorption spectroscopy (AX-ATAS). Our analysis showed that the AX-ATAS response is confined to the subfemtosecond timescale, eliminating any hydrogen atom motion and demonstrating experimentally that the 1b 1 splitting in the x-ray emission spectrum is related to dynamics and is not evidence of two structural motifs in ambient liquid water.

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

Unraveling the Molecular Origin of Prey-Wrapping Spider Silk's Unique Mechanical Properties and Assembly Process Using NMR

Prey wrapping spider silk's unique mechanical properties are investigated confirming the silk's high degree of extensibility and superior toughness compared to other types of spider silk. For the first time, the pre-spinning dope phase is studied in isotope-enriched intact aciniform (AC) silk glands using solution NMR that reveals a combination of α-helical domains linked by disordered random coil chains consistent with previously proposed “beads-on-a-string” models. The model is further refined through the AlphaFold2 protein structure prediction tool. Finally, extensive magic angle spinning (MAS) solid-state (SS) NMR data for isotopically-enriched fibers is used to refine the structural model for AC silk from two species, A. aurantia and A. argentata. The SSNMR data shows that the AC silk fibers are highly α-helical, coiled-coil in structure but, also exhibit significant β-sheet components that can be traced back to the Gly-rich disordered linker regions in the pre-spinning dope phase that are converted to β-sheet structures during fiber formation. This combination of mechanical and structural characterization enhances the understanding of AC silk's liquid-to-solid transition and structure-mechanics relationship. In conclusion, these prey wrap silk results and models will provide the basis for the design of biomimetic materials inspired by the AC spider silk system.

36 MATERIALS SCIENCE↗
DOE OSTI2022

Ultra-selective membrane composed of charge-stabilized fixed carrier and amino acid-based ionic liquid mobile carrier for highly efficient carbon capture

Membrane technology has been extensively studied for CO 2 capture applications, especially for flue gas sources. In the past decades, facilitated transport membranes (FTMs) have made breakthroughs overcoming the permeance-selectivity trade-off upper bound restricting traditional CO 2 separation membranes, but are still facing challenges towards practical applications, including limited performance, such as insufficient CO 2 /N 2 selectivity to achieve 95 % CO 2 dry-base purity by one step separation, and long-term stability. Herein, we designed and fabricated a novel FTM structure containing an ionic liquid (1-ethyl-3-methylimidazolium aminoacetate, [Emim][Gly]) as mobile CO 2 -carrier and a polymeric amine (polyethyleneimine, PEI) as fixed CO 2 -carrier. In this study, the fixed carrier is confined within a carbon nanotube (CNT) framework of 230 nm thickness via electrostatic forces adjusted by a polyelectrolyte (polystyrene sulfonate, PSS), while the mobile carrier diffuses freely within the CNT framework. After optimization of the membrane recipe and spray-coating fabrication procedure, following our previous work, the resulting CNT-PSS-PEI ~ IL membranes demonstrated an ultra-high CO 2 /N 2 selectivity up to 1,000 with CO 2 permeance up to 2,400 GPU (Gas Permeation Unit, 1 GPU = 3.348 × 10 -10 mol·s -1 ·m -2 ·Pa -1 ) under vacuum operation condition. Furthermore, one 100-cm 2 flat sheet membrane sample was prepared and exhibited one-stage CO 2 enrichment from 15 % to 95 % purity (dry-base) for the first time amongst all reported CO 2 separation membranes. The membrane retained a stable performance over 50-h operation period under vacuum condition. The extraordinary CO 2 separation performance illustrates the great potential of the CNT-PSS-PEI ~ IL membranes for flue gas carbon capture application.

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

Structural characterization of sugarcane lignins extracted from different protic ionic liquid pretreatments

Previous studies based on protic ionic liquids have shown them to be effective for extracting lignins from biomass. Moreover, ILs synthesized with amine cations can promote the insertion of nitrogen into lignin, favoring its valorization for industrial applications. Here, lignins extracted from sugarcane bagasse by eight different PILs were analyzed. Several techniques of characterization were employed (2D HSQC NMR, 31P NMR, GPC, Py-GC/MS, TGA, FTIR, and elemental analysis), seeking a deeper understanding of lignin structure post-PILS treatment. In particular, lignin obtained with acetate and lactate anions with monoethanolamonium cations showed similar structural compositions. Lignins obtained from ethylenediamonium cation combined anions acetate and lactate exhibited relatively higher molecular weights compared to the other lignins studied. The elemental analysis of these lignins further revealed the presence of nitrogen in their structures, suggesting that amination occurred. In this work we used PILs not yet reported in the literature for this purpose, and consequently obtained lignins with specific structures. Overall, our results demonstrated that depending on the PILs used, lignins with distinct properties could be obtained which might be used for value-added applications. This work covered new advances in elucidation of lignin chemical structure and opened a new path for lignin valorization.

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

Identification of efficient and stable solid-solid electrochemical interfaces for Li-S batteries. Final Report

The objective of this work is to develop a better understanding of the interactions between battery components that drive the observed battery performance. Using computational tools at the atomistic level, we elucidate molecular, crystal, and amorphous structures, interfacial structures and reactions, and electrochemical behavior. We focus on Li metal anodes, liquid electrolytes, solid electrolytes, S cathodes, and metal-oxide cathodes, and their corresponding interfaces. Within liquid electrolytes, we study their bulk structure and dynamics, as a function of salt concentration. The analysis allows to understand the formation of localized high concentrated regions as the solvent is substituted by a weakly interacting diluent species. Analysis of the electron distribution yields the reactive behavior at interfaces. Moreover, we employ tools that allow the characterization of ion and electron transport at interfaces, during the processes of solvation and desolvation occurring while cycling. For Sulfur cathodes we elucidate synthesis processes and the role of carbon-sulfur and nitrogen-sulfur interactions not only to provide electronic conductivity but also to define alternative reaction pathways during discharge of Li-S batteries. Full-cell simulations offer new insights regarding the integration of all the battery components. Our studies also include solid electrolytes and their interfaces with Li metal, S-C, and metal-oxide electrodes.

25 ENERGY STORAGE↗
DOE OSTI2021

Evidence of a liquid–liquid transition in a glass-forming ionic liquid

A liquid–liquid transition (LLT) is a transformation from one liquid to another through a first-order transition. The LLT is fundamental to the understanding of the liquid state and has been reported in a few materials such as silicon, phosphorus, triphenyl phosphite, and water. Furthermore, it has been suggested that the unique properties of materials such as water, which is critical for life on the planet, are linked to the existence of the LLT. However, the experimental evidence for the existence of an LLT in many molecular liquids remains controversial, due to the prevalence and high propensity of the materials to crystallize. In this work, we show evidence of an LLT in a glass-forming trihexyltetradecylphosphonium borohydride ionic liquid that shows no tendency to crystallize under normal laboratory conditions. We observe a step-like increase in the static dielectric permittivity at the transition. Furthermore, the sizes of nonpolar local domains and ion-coordination numbers deduced from wide-angle X-ray scattering also change abruptly at the LLT. We independently corroborate these changes in local organization using Raman spectroscopy. The experimental access to the evolution of local order and structural dynamics across a liquid–liquid transition opens up unprecedented possibilities to understand the nature of the liquid state.

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

Material Model Parameters Optimization in Liquid Mercury Target Dynamics Simulation With Machine Learning Surrogates

A pulsed spallation target is subjected to very short (∼0.7μs) but intense loads (23.3 kJ) from repeated proton pulses, which knock away neutrons from the mercury atoms’ nuclei for a wide range application in physics, engineering, medicine, petroleum exploration, biology, chemistry, etc. The effect of this pulsed loading on the stainless-steel target module which contains the flowing mercury target material is difficult to predict not only due to its short but intense explosive-like physical reaction, but also the nonlinear material behavior of the liquid mercury in the structure. Injecting small helium bubbles in the mercury has been an efficient method of mitigating the pressure wave at high power level stage. However, prediction of the resultant loading on the target is more difficult when helium gas is intentionally injected into the mercury. A 2-phase material model that incorporates the Rayleigh-Plesset (R-P) model is expected to address this complex multi-physics dynamics problem by including the bubble dynamics in the liquid mercury. A parameter sensitivity study was firstly employed to understand their impact on the simulation strains. The investigated parameters included E, μ, γ, σ, n, VFgas, and gas cumulative volume curve control parameters a and b. Verification and validation results from sparse polynomial expansions (SPE) method and directional Gaussian smoothing (DGS) optimization show that the surrogate model had training error of ∼7% and validation error of ∼15%, indicating that machine learning methods and surrogate models can help optimize the uncertain parameters in the complex 2-phase material model. This approach is expected to fill the knowledge gap between unknown liquid-gas mixture material model and measured vessel strain responses.

Lin, Lianshan↗
DOE OSTI2022

MISPR : an open-source package for high-throughput multiscale molecular simulations

Computational tools provide a unique opportunity to study and design optimal materials by enhancing our ability to comprehend the connections between their atomistic structure and functional properties. However, designing materials with tailored functionalities is complicated due to the necessity to integrate various computational-chemistry software (not necessarily compatible with one another), the heterogeneous nature of the generated data, and the need to explore vast chemical and parameter spaces. The latter is especially important to avoid bias in scattered data points-based models and derive statistical trends only accessible by systematic datasets. Here, we introduce a robust high-throughput multi-scale computational infrastructure coined MISPR (Materials Informatics for Structure–Property Relationships) that seamlessly integrates classical molecular dynamics (MD) simulations with density functional theory (DFT). By enabling high-performance data analytics and coupling between different methods and scales, MISPR addresses critical challenges arising from the needs of automated workflow management and data provenance recording. The major features of MISPR include automated DFT and MD simulations, error handling, derivation of molecular and ensemble properties, and creation of output databases that organize results from individual calculations to enable reproducibility and transparency. In this work, we describe fully automated DFT workflows implemented in MISPR to compute various properties such as nuclear magnetic resonance chemical shift, binding energy, bond dissociation energy, and redox potential with support for multiple methods such as electron transfer and proton-coupled electron transfer reactions. The infrastructure also enables the characterization of large-scale ensemble properties by providing MD workflows that calculate a wide range of structural and dynamical properties in liquid solutions. MISPR employs the methodologies of materials informatics to facilitate understanding and prediction of phenomenological structure–property relationships, which are crucial to designing novel optimal materials for numerous scientific applications and engineering technologies.

36 MATERIALS SCIENCE↗
DOE OSTI2024

Universality of Collective Density Fluctuations in Liquids at and away from Equilibrium – An Integrated Neutron Scattering, Theoretical, and Computational Study (Final Technical Report)

Liquids, ubiquitous on earth, are prototypical disordered condensed matter. As the mediate phase between non-structured gases and ordered solids, liquids p ossess intriguing complexities that are absent on either extreme counter phase, especially when they are driven away from equilibrium. Despite decades of studies, the physics of liquids still keeps surprising us and challenges our understanding of condensed matter, emergent phenomena, and complex systems. Among the many interesting properties, the nature of collective dynamics in liquids remains elusive because of both the presence of strong interactions and the absence of translational invariance. This project studied the universal features of collective density fluctuations (or collective modes), which include both phonon-like excitations and collective relaxations, in the Q-dependent generalized hydrodynamic regime in liquids and liquid-like matter both at and away from equilibrium. We performed systematic studies of the fast collective dynamics (including both phonon-like excitations and collective relaxations) of four classes of liquids characterized by different fragilities and interactions using synergistically integrated coherent neutron scattering experiments (Inelastic Neutron Scattering (INS), Quasi-Elastic Neutron Scattering (QENS), and Neutron Spin Echo (NSE)), a ViscoElastic Hydrodynamics theory we have been developing, andMolecular Dynamics (MD) simulations and our analysis package LiquidLib. The transformative knowledge of the collective dynamics in liquids is important to understand and ultimately control the transport of mass, energy, and charge in liquids, both at and away from thermal equilibrium, which may will likely impact a wide range of disciplines.

36 MATERIALS SCIENCE↗
DOE OSTI2021

MB-Fit: Software infrastructure for data-driven many-body potential energy functions

Many-body potential energy functions (MB-PEFs), which integrate data-driven representations of many-body short-range quantum mechanical interactions with physics-based representations of many-body polarization and long-range interactions, have recently been shown to provide high accuracy in the description of molecular interactions from the gas to the condensed phase. Here, we present MB-Fit, a software infrastructure for the automated development of MB-PEFs for generic molecules within the TTM-nrg (Thole-type model energy) and MB-nrg (many-body energy) theoretical frameworks. Besides providing all the necessary computational tools for generating TTM-nrg and MB-nrg PEFs, MB-Fit provides a seamless interface with the MBX software, a many-body energy and force calculator for computer simulations. Given the demonstrated accuracy of the MB-PEFs, particularly within the MB-nrg framework, we believe that MB-Fit will enable routine predictive computer simulations of generic (small) molecules in the gas, liquid, and solid phases, including, but not limited to, the modeling of quantum isomeric equilibria in molecular clusters, solvation processes, molecular crystals, and phase diagrams.

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

Secondary electron emission measurements from imidazolium-based ionic liquids

The electron-induced secondary electron emission (SEE) yields of imidazolium-based ionic liquids are presented for primary electron beam energies between 30 and 1000 eV. These results are important for understanding plasma synthesis of nanoparticles in plasma discharges with an ionic liquid electrode. Due to their low vapor pressure and high conductivity, ionic liquids can produce metal nanoparticles in low-pressure plasmas through reduction of dissolved metal salts. In this work, the low vapor pressure of ionic liquids is exploited to directly measure SEE yields by bombarding the liquid with electrons and measuring the resulting currents. The ionic liquids studied are [BMIM][Ac], [EMIM][Ac], and [BMIM][BF 4 ]. The SEE yields vary significantly over the energy range, with maximum yields of around 2 at 200 eV for [BMIM][Ac] and [EMIM][Ac], and 1.8 at 250 eV for [BMIM][BF 4 ]. Molecular orbital calculations indicate that the acetate anion is the likely electron donor for [BMIM][Ac] and [EMIM][Ac], while in [BMIM][BF 4 ], the electrons likely originate from the [BMIM] + cation. The differences in SEE yields are attributed to varying ionization potentials and molecular structures of the ionic liquids. These findings are essential for accurate modeling of plasma discharges and understanding SEE mechanisms in ionic liquids.

ionic liquids↗
DOE OSTI2025

Development, Verification, and Validation of an OpenFOAM-Based Solver for Modeling Inertial Fusion Energy Chambers

Our work seeks to introduce a computational tool tailored to the physics of inertial fusion energy chambers, in particular, those concepts based on thick liquid walls. In this approach, the structural materials are protected by several neutron mean-free-paths of renewable liquid and thus will be able to survive much longer than un-shielded walls, with virtually all structures lasting for the life of the plant and enabling the use of commercially available and qualified materials. The OpenFOAM-based solver named rhoCentralFoam has been used as a starting point. rhoCentralFoam belongs to the standard OpenFOAM solver toolset. It is a high-speed, explicit compressible flow solver with shock-capturing capability. While the main features have been retained, the solver had to be restructured to make use of tabular data for equations of states, a necessary addition to model the complex thermo-physical properties of ionized gasses. This entailed the need to change the independent state variables used by the solver, resulting in a new thermodynamic library and slightly different solution algorithm. Moreover, a radiation heat transfer model based on the P-1 approximation was added to the solver. The solver is verified against an analytical solution from the Sedov-Taylor-Neumann test problem to showcase the ability of the hydrodynamic solvers to handle strong shocks, whereas the P-1 model was verified using a simple one-dimensional problem with an analytical solution. Additionally, a validation case involving shock-wave propagation through jet array is presented, and the results are compared with experimental data from the open literature. Lastly, in order to showcase the utility of the solver for practical cases, we applied the refined solver to two representative scenarios: gas venting within the HYLIFE-II chamber and the compression of the gas following the partial ablation of the liquid wall.

Chamber dynamics↗
DOE OSTI2023

World beyond the nearest neighbors

The structure beyond the nearest neighbor atoms in liquid and glass is characterized by the medium-range order (MRO). In the conventional approach, the MRO is considered to result directly from the short-range order (SRO) in the nearest neighbors. To this bottom–up approach starting with the SRO, we propose to add a top–down approach in which global collective forces drive liquid to form density waves. The two approaches are in conflict with each other, and the compromise produces the structure with the MRO. The driving force to produce density waves provides the stability and stiffness to the MRO, and controls various mechanical properties. This dual framework provides a novel perspective for description of the structure and dynamics of liquid and glass.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2020

Laminar drag reduction in microchannels with liquid infused textured surfaces

Liquid infused surfaces (LIS) achieve non-wetting properties through asperity structures containing pockets of a lubricating liquid rather than air. Although studies have demonstrated several potential applications of LIS, the effect of liquid infused texture on the flow of liquid in microchannels has not been systematically addressed. The present work focuses on the effect of liquid-infused texture geometry and the infused liquid properties on the flow characteristics of a bulk fluid in rectangular and round microchannel geometries. Considering four different textured geometries of transverse ridges, longitudinal ridges, in-lined square posts and staggered squares, a detailed parametric study is presented to assess the resulting drag reduction and friction factor, expressed in terms of the Poiseuille number, on the flow in microchannels. The effect of geometrical parameters is investigated for different asperity solid fraction (Φ s ), asperity height (h s ) and channel thickness (H) or tube diameter (D), while the effect of the infused fluid is studied in terms of a ratio of its viscosity relative to the viscosity of the flowing bulk fluid. Based on the results of the parametric study an extension of the conventional Moody diagram for non-wetting surfaces is developed. Finally, a design of experiments study is conducted to relate the drag reduction and Poiseuille number to the various governing parameters, using which conditions that maximize drag reduction are identified.

42 ENGINEERING↗
DOE OSTI2024

Liquid Crystal Orientation and Shape Optimization for the Active Response of Liquid Crystal Elastomers

Liquid crystal elastomers (LCEs) are responsive materials that can undergo large reversible deformations upon exposure to external stimuli, such as electrical and thermal fields. Controlling the alignment of their liquid crystals mesogens to achieve desired shape changes unlocks a new design paradigm that is unavailable when using traditional materials. While experimental measurements can provide valuable insights into their behavior, computational analysis is essential to exploit their full potential. Accurate simulation is not, however, the end goal; rather, it is the means to achieve their optimal design. Such design optimization problems are best solved with algorithms that require gradients, i.e., sensitivities, of the cost and constraint functions with respect to the design parameters, to efficiently traverse the design space. In this work, a nonlinear LCE model and adjoint sensitivity analysis are implemented in a scalable and flexible finite element-based open source framework and integrated into a gradient-based design optimization tool. To display the versatility of the computational framework, LCE design problems that optimize both the material, i.e., liquid crystal orientation, and structural shape to reach a target actuated shapes or maximize energy absorption are solved. Multiple parameterizations, customized to address fabrication limitations, are investigated in both 2D and 3D. The case studies are followed by a discussion on the simulation and design optimization hurdles, as well as potential avenues for improving the robustness of similar computational frameworks for applications of interest.

42 ENGINEERING↗
DOE OSTI2020

Temperature dependence of nuclear quantum effects on liquid water via artificial neural network model based on SCAN meta-GGA functional

We investigate the temperature dependence of nuclear quantum effects (NQEs) on structural and dynamic properties of liquid water by training a neural network force field using first-principles molecular dynamics (FPMD) based on the strongly constrained and appropriately normed meta-generalized gradient approximation exchange-correlation approximation. The FPMD simulation based on density functional theory has become a powerful computational approach for studying a wide range of condensed phase systems. However, its large computational cost makes it difficult to incorporate NQEs in the simulation and investigate temperature dependence of various properties. To circumvent this difficulty, we use an artificial neural network model and employ the thermostatted ring polymer MD approach for studying the temperature dependence of NQEs on various properties. The NQEs generally bring the radial distribution functions closer to the experimental measurements. Translational diffusivity and rotational dynamics of water molecules are both slowed down by the NQEs. The competing inter-molecular and intra-molecular quantum effects on hydrogen bonds, as discussed by Habershon, Markland, and Manolopoulos [J. Chem. Phys. 131(2), 024501 (2019)], can explain the observed temperature dependence of the NQEs on the dynamical properties in our simulation.

Yao, Yi (ORCID:0000000299576979)↗
DOE OSTI2020

Fermi-liquid theory and divergences of the two-particle irreducible vertex in the periodic Anderson lattice

Here we analyze the divergences of the irreducible vertex function in dynamical mean field theory, which may indicate either a nonphysical breakdown of the perturbation theory or a response to some physical phenomenon. To investigate this question, we construct a quasiparticle vertex from the diverging irreducible vertex functions. This vertex describes the scattering between quasiparticles and quasiholes in a Fermi liquid. We show that the quasparticle vertex does not diverge in the charge channel, wherein the irreducible vertex does diverge; and we show that the quasiparticle vertex does diverge in the spin channel, wherein the irreducible vertex does not diverge. This divergence occurs at the Mott transition, wherein the Fermi-liquid theory breaks down. Both the half filled Hubbard and Anderson lattices are investigated. In general, our results support that the divergences of the irreducible vertex function do not indicate a nonphysical failure of the perturbation theory. Instead, the divergences are the mathematical consequence of inverting a matrix (the local charge susceptibility) which accumulates increasingly negative diagonal elements as the Hubbard interaction suppresses charge fluctuations. Indeed, we find that the first divergences of the irreducible vertex in both Hubbard and Anderson lattices occurs near the maximum magnitude of the (negative) vertex-connected part of the charge susceptibility.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗