Structure and function of LCI1: a plasma membrane CO[subscript 2] channel in the Chlamydomonas CO[subscript 2] concentrating mechanism
Abstract not provided
Kono, Alfredo↗
SEARCH · Engineering Papers
Results for “molecular functionalization”
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.
Velocity structure functions in multiphase turbulence: interpreting kinematics of Hα filaments in cool-core clusters
ABSTRACT The central regions of cool-core galaxy clusters harbour multiphase gas, with gas temperatures ranging from $10$ to $10^7\, \mathrm{K}$. Feedback from active galactic nuclei jets prevents the gas from undergoing a catastrophic cooling flow. However, the exact mechanism of this feedback energy input is unknown, mainly due to the lack of velocity measurements of the hot-phase gas. However, recent observations have measured the velocity structure functions (VSFs) of the cooler molecular (${\sim} 10\, \mathrm{K}$) and Hα filaments (${\sim} 10^4\, \mathrm{K}$) and used them to indirectly estimate the motions of the hot phase. In the first part of this study, we conduct high-resolution (3843–15363 resolution elements) simulations of homogeneous isotropic subsonic turbulence, without radiative cooling. We analyse the second-order velocity structure functions (VSF2) in these simulations and study the effects of varying spatial resolution, the introduction of magnetic fields, and the effect of projection along the line of sight (LOS) on it. In the second part of the study, we analyse high-resolution (7683 resolution elements) idealized simulations of multiphase turbulence in the intracluster medium from the companion study Mohapatra et al. We compare the VSF2 for both the hot ($T\sim 10^7\, \mathrm{K}$) and cold ($T\sim 10^4\, \mathrm{K}$) phases and find that their amplitude depends on the density contrast between the phases. They have similar scaling with separation, but introducing magnetic fields steepens the VSF2 of only the cold phase. We also find that projection along the LOS steepens the VSF2 for the hot phase and mostly flattens it for the cold phase.
Mohapatra, Rajsekhar (ORCID:0000000216007552)↗
Editorial: Structure and mechanism of microbial membrane active transporters
Membrane active transporters play essential roles in microbial physiology. They couple energy transduction to conformational changes that drive translocation of nutrients, substrates and ions, as well as molecular communication. The structure and function of microbial membrane active transporters are highly diverse. Typical examples include the primary active transporters in the ATP-binding cassette (ABC) superfamily (Thomas and Tampé, 2020; Davidson et al., 2008; Locher et al., 2002), the secondary active transporters in the Major Facilitator Superfamily (MFS) (Drew et al., 2021; Kaback and Guan, 2019), and the ligand-gated porins in the TonB-dependent transporter (TBDT) family (Klebba et al., 2021). As structural, proteogenomic, and computational methods advance, active transporters are increasingly recognized as dynamic molecular machines whose mechanisms can now be visualized and modeled with remarkable precision, building on decades of biochemical and biophysical discovery that established the foundations of this field. The transporter studies recruited in this Research Topic provide us with new insights into the field including structure-function of sugar transporters in yeast, structural prediction and classification of ABC complexes in Bacillus subtilis, Type VI Secretion System (T6SS) in Bacteroides fragilis, amino acids uptake in Escherichia coli and bacterial spore germination.
mechanism↗
DeePKS + ABACUS as a Bridge between Expensive Quantum Mechanical Models and Machine Learning Potentials
Recently, the development of machine learning (ML) potentials has made it possible to perform large-scale and long-time molecular simulations with the accuracy of quantum mechanical (QM) models. However, for different levels of QM methods, such as density functional theory (DFT) at the meta-GGA level and/or with exact exchange, quantum Monte Carlo, etc., generating a sufficient amount of data for training an ML potential has remained computationally challenging due to their high cost. In this work, we demonstrate that this issue can be largely alleviated with Deep Kohn–Sham (DeePKS), an ML-based DFT model. DeePKS employs a computationally efficient neural network-based functional model to construct a correction term added upon a cheap DFT model. Upon training, DeePKS offers closely matched energies and forces compared with high-level QM method, but the number of training data required is orders of magnitude less than that required for training a reliable ML potential. As such, DeePKS can serve as a bridge between expensive QM models and ML potentials: one can generate a decent amount of high-accuracy QM data to train a DeePKS model and then use the DeePKS model to label a much larger amount of configurations to train an ML potential. Further, this scheme for periodic systems is implemented in a DFT package ABACUS, which is open source and ready for use in various applications.
71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
Statistical analysis of HAADF-STEM images to determine the surface coverage and distribution of immobilized molecular complexes
The surface immobilization of molecular catalysts is attractive because it combines the benefits of homogeneous and heterogeneous catalysis. However, determining the surface coverage and distribution of a molecular catalyst on a solid support is often challenging, inhibiting our ability to design improved catalytic systems. Here, in this work, we demonstrate that the combination of scanning transmission electron microscopy (STEM) and image analysis of the individual positions of heavy atoms in transition metal complexes via a convolutional neural network (CNN) allows statistically robust determination of the surface coverage and distribution of immobilized molecular catalysts. These observations provide information about how changes in the functionalization conditions, attachment group, and structure of the molecular catalyst affect the surface coverage and distribution, providing insight into the chemical mechanism of surface immobilization. The method could be generally valuable for correlating the surface coverage and distribution to the activity, selectivity, and stability of a catalytic system.
HAADF-STEM↗
Enhanced liquid metal wetting on oxide surfaces via patterned particles
Liquid metals typically do not intrinsically wet oxide surfaces due to chemical dissimilarity, preventing the formation of well-bonded interfaces in many applications. This paper investigates how the wetting and spreading of liquid silver (Ag) droplets on oxide surfaces can be enhanced via patterned nickel (Ni) particles using multiscale modeling. First, a force field (FF) for an Ag/yttria-stabilized-zirconia system was developed by fitting the interfacial binding energy to density-functional theory (DFT) calculations. Then, molecular dynamics (MD) simulations using these FF parameters revealed that the good wetting of Ag on Ni caused Ag droplets to be attracted to the Ni particles, promoting spreading. Thus, the optimal hexagonal Ni pattern was shown to enhance the wetting area by 224% compared to a bare oxide surface. Second, an analytical model predicting the maximum possible enhancement in wetting area as a function of the intrinsic wetting angle, size of the Ag droplet, and size of the Ni particles was derived and validated with MD simulations, to generalize this method to other substrates. Here, both the analytical model and MD results suggested that the ability of Ni patterns to enhance the Ag wetting area becomes more significant as the wettability of substrate becomes worse.
36 MATERIALS SCIENCE↗
First-principles investigation of structure-property relationships in stable and metastable MXenes
Understanding the structure–property relationships in layered transition-metal carbides or nitrides, known as MXenes, is of critical importance for their rational design, synthesis, and application. However, the vast chemical and structural diversity of MXenes, stemming from their wide range of M and X elements, surface terminations, and different atomic coordination environments, makes it challenging to clearly understand these structure–property relationships. Here, in this work, we perform first-principles density functional theory (DFT) calculations and molecular dynamics (MD) simulations to comprehensively investigate the stability and a variety of physical properties of MXenes with different coordination environments. Using Ti- and Mo-based carbide MXenes as model systems, energetic calculations reveal that Ti-based MXenes are most stable in octahedral coordination, whereas Mo-based MXenes preferentially adopt prismatic coordination. This fundamental difference in preferred atomic coordination gives rise to markedly distinct properties between these two systems as a function of the fraction of octahedral and prismatic sites. For instance, the in-plane stiffness of Ti-based MXenes increases as octahedral coordination becomes dominant, but it decreases in the Mo-based MXenes under the same conditions. Additional stability analyses based on mechanical, lattice-dynamical, and temperature-dependent thermodynamic properties demonstrate that many metastable MXenes not only satisfy the strict stability criteria but can also undergo phase transitions among different structures and even become stabilized at elevated temperatures. Although surface terminations, such as F and O atoms, do not alter the energetic ordering or the overall stiffness trends among stable and metastable MXenes, they influence other material properties. For instance, O termination can induce semiconducting behavior in both stable and metastable Ti 2 CO 2 MXenes. This study significantly advances the fundamental understanding of structure–property relationships in MXenes and provides valuable guidance for developing coordination-based design principles to precisely engineer MXenes with improved properties.
Oyeniran, Noah [University of Alabama, Tuscaloosa,↗
Modeling Disorder in Pyrochlores and Other Anion-Deficient Fluorite Structural Derivative Oxides
Their very flexible chemistry gives oxide materials a richness in functionality and wide technological application. A specific group of oxides that have a structure related to fluorite but with less oxygen, termed anion-deficient fluorite structural derivatives and with pyrochlores being the most notable example, has been shown to exhibit a diversity of useful properties. For example, the possibility to undergo a transition from an ordered to disordered state allows these oxides to have high radiation tolerance. Atomistic-scale calculations in the form of molecular dynamics (MD) and density functional theory (DFT) have been extensively used to understand what drives this order/disorder transition. Here we give a brief overview of how atomistic-scale calculations are utilized in modeling disorder in pyrochlores and other anion-deficient fluorite structural derivatives. We discuss the modeling process from simple point defects to completely disordered structures, the dynamics during the disordering process, and the use of mathematical models to generate ordered solid-solution configurations. We also attempt to identify the challenges in modeling short range order and discuss future directions to more comprehensive models of the disordered structures.
37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
Insights into Probing the Effect of Molecular Weight of Poly(carboxybetaine methacrylate) on the Performance of Forward Osmosis Desalination
Zwitterionic polymers have proven to be a promising nonfouling material that can be applied in the design of selective layers of thin film composite (TFC) membranes. Extending the permeability and usage of TFC membranes have attracted increasing interest in membrane-based desalination processes since water-flux reduction associated with biofouling nowadays persists as a common challenge. By virtue of its strong hydration, this polymer category is very useful to counteract biofouling in marine and biomedical systems, but the benefits from their application in membrane technology are still emerging. The efficacy of the nonfouling property as a function of the polymer’s molecular weight remains unknown. In pursuit of that vision, this study fosters new scientific insights via probing different molecular weights of poly(carboxybetain methacrylate) (PCBMA) coated on the surface as a selective layer for the prepared TFC membranes. The coated zwitterionic membranes (zM) exhibited excellent performance in preventing water flux decay in a bench-scale forward osmosis system. The prepared zM membranes revealed enhanced hydrophilic properties and retained their operational water flux when compared to the control. Our results suggest that using an intermediate-size molecular weight (PCBMA M n 50,000) will result in the best operational performance. The intermediate size resulted in the lowest flux decline rate (R t ) of 0.01 ± 0.001 (zM-50) when compared to the unmodified control membrane 0.56 ± 0.071 (M0) after using a model BSA foulant solution. Furthermore, all coated membranes exhibited similar trends in the observed reverse salt flux profiles, as well. In conclusion, the constructed zM membranes will serve as a model to develop further selective layers in the construction of TFC membranes.
36 MATERIALS SCIENCE↗
Ab initio generalized Langevin equation
We introduce a machine learning–based approach called ab initio generalized Langevin equation (AIGLE) to model the dynamics of slow collective variables (CVs) in materials and molecules. In this scheme, the parameters are learned from atomistic simulations based on ab initio quantum mechanical models. Force field, memory kernel, and noise generator are constructed in the context of the Mori–Zwanzig formalism, under the constraint of the fluctuation–dissipation theorem. Combined with deep potential molecular dynamics and electronic density functional theory, this approach opens the way to multiscale modeling in a variety of situations. Here, we demonstrate this capability with a study of two mesoscale processes in crystalline lead titanate, namely the field-driven dynamics of a planar ferroelectric domain wall, and the dynamics of an extensive lattice of coarse-grained electric dipoles. In the first case, AIGLE extends the reach of ab initio simulations to a regime of noise-driven motions not accessible to molecular dynamics. In the second case, AIGLE deals with an extensive set of CVs by adopting a local approximation for the memory kernel and retaining only short-range noise correlations. The scheme is computationally more efficient than molecular dynamics by several orders of magnitude and mimics the microscopic dynamics at low frequencies where it reproduces accurately the dominant far-infrared absorption frequency.
97 MATHEMATICS AND COMPUTING↗
Renormalized site density functional theory for models of ion hydration
Development of accurate statistical mechanics models of molecular liquid systems is a problem of great practical and fundamental importance. Site-density functional theory (SDFT) is one of the promising directions in this area, but its success hinges upon the ability to efficiently reconcile co-existence of two distinct intra- and inter-molecular interaction regimes in a molecular liquid. Renormalized formulation of site-density functional theory (RSDFT), that we have recently developed, resolves this problem by introducing an additional potential field variable that decouples two interaction scales and maps the molecular liquid problem onto the effective simple liquid mixture. This work provides a critical assessment of RSDFT for the hydrated ion system - a problem, which historically has always been one of the most difficult cases for SDFT applications. Using two-site model of water, we perform a comprehensive analysis of hydrated alkali metal and halogen ions, including both structural and free energy based characteristics. The results indicate that RSDFT provides a significant improvement over conventional 3D-RISM implementations, and may prove useful in coarse grained simulations based on two-site solvent models.
Chuev, Gennady N.↗
Improving scalability of electronic structure code for molecular simulations in the presence of environment
A scalable density functional electronic code with Gaussian basis set, called UTEP-NRLMOL, is developed to perform simulations of molecular systems in the presence of the environment with particular attention to the memory requirements. In the electronic structure calculations, the memory and computation time are proportional to the number of atoms. Memory requirements for density functional calculations scale as N*N, where N is the number of atoms. While the recent advances in HPC offer platforms with large numbers of cores, the limited amount of memory available on a given node and poor scalability of the electronic structure codes hinder their efficient usage of these platforms. We have introduced new scaling and parallelization paradigms using MPI-3 shared-memory functionality combined with usage of sparse algebra and storage of matrices in sparse format. This extends the range of applicability of the UTEP-NRLMOL code to large systems over 10,000 atoms, or using up to 67,000 basis functions, and making use of HPC architectures using over 6,000 processors utilizing all available cores. We have also interfaced code with effective fragment potential and polarizable continuum model libraries. The code was used in simulations of several applications which are published in reputed scientific journals.
74 ATOMIC AND MOLECULAR PHYSICS↗
First-principles molten salt phase diagrams through thermodynamic integration
Precise prediction of phase diagrams in molecular dynamics simulations is challenging due to the simultaneous need for long time and large length scales and accurate interatomic potentials. Here, we show that thermodynamic integration from low-cost force fields to neural network potentials trained using density-functional theory (DFT) enables rapid first-principles prediction of the solid–liquid phase boundary in the model salt NaCl. We use this technique to compare the accuracy of several DFT exchange–correlation functionals for predicting the NaCl phase boundary and find that the inclusion of dispersion interactions is critical to obtain good agreement with experiment. Importantly, our approach introduces a method to predict solid–liquid phase boundaries for any material at an ab initio level of accuracy, with the majority of the computational cost at the level of classical potentials.
37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
Atomic structure and grain shape evolution of nanodiamond during annealing in oxidizing atmosphere from neutron diffraction and MD simulations
In this work, real space pair-distribution function (PDF) analysis in conjunction with molecular dynamics (MD) simulations was applied to characterize crystal structure, shape and surface structure of diamond nanoparticles annealed in air up to 880 °C. MD simulated models of individual nanodiamond particles of various shapes and sizes were used to calculate theoretical radial distribution functions G(r) for powders with different grain size distributions. The calculated curves, showing distinct differences for the explored parameter space, were compared to the experimental G(r) functions obtained from neutron scattering for annealed detonated nanodiamond samples in a search for a model best describing the actual powders. It was found that as-synthesized nanodiamond grains are terminated by equally abundant (100), (110) and (111) crystallographic faces and their structure does not change upon annealing in air up to 280 °C. The shape and size of grains annealed above 480 °C changes due to surface etching. The etching rate differs between crystals faces and depends on temperature, resulting in the observed differences in grain shape. Below 680 °C the share of the (110) facets gradually increases and at 680 °C the grains become rhombic dodecahedra terminated solely by (110) surfaces. Further increase of the annealing temperature promotes formation of (111) surfaces and at 880 °C the preferred grain shape becomes an octahedron terminated by a variant of (111) surfaces possessing three dangling sp3 bonds per atom. At lower temperatures all grains become smaller and the size distribution moves towards the lower values. Furthermore, as the annealing temperature increases, the smallest grains disappear entirely and the average grain size grows back.
36 MATERIALS SCIENCE↗
Boron adatom adsorption on graphene: A case study in computational chemistry methods for surface interactions
Though weak surface interactions and adsorption can play an important role in plasma processing and materials science, they are not necessarily simple to model. A boron adatom adsorbed on a graphene sheet serves as a case study for how carefully one must select the correct technique from a toolbox of computational chemistry methods. Using a variety of molecular dynamics potentials and density functional theory functionals, we evaluate the adsorption energy, investigate barriers to adsorption and migration, calculate corresponding reaction rates, and show that a surprisingly high level of theory may be necessary to verify that the system is described correctly.
37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
The symmetric quasi-classical model using on-the-fly time-dependent density functional theory within the Tamm–Dancoff approximation
The primary computational challenge when simulating nonadiabatic ab initio molecular dynamics is the unfavourable compute costs of electronic structure calculations with molecular size. Simple electronic structure theories, like time-dependent density functional theory within the Tamm–Dancoff approximation (TDDFT/TDA), alleviate this cost for moderately sized molecular systems simulated on realistic time scales. Although TDDFT/TDA does have some limitations in accuracy, an appealing feature is that, in addition to including electron correlation through the use of a density functional, the cost of calculating analytic nuclear gradients and nonadiabatic coupling vectors is often computationally feasible even for moderately sized basis sets. Here in this work, some of the benefits and limitations of TDDFT/TDA are discussed and analysed with regard to its applicability as a ‘back-end’ electronic structure method for the symmetric quasi-classical Meyer–Miller model (SQC/MM). In order to investigate the benefits and limitations of TDDFT/TDA, SQC/MM is employed to predict and analyse a prototypical example of excited-state hydrogen transfer in gas-phase malonaldehyde. Then, the ring-opening dynamics of selenophene are simulated, which highlight some of the deficiencies of TDDFT/TDA. Additionally, some new algorithms are proposed that speed up the calculation of analytic nuclear gradients and nonadiabatic coupling vectors for a set of excited electronic states.
molecular dynamics↗
Liquid–Vapor Phase Equilibrium in Molten Aluminum Chloride (AlCl 3 ) Enabled by Machine Learning Interatomic Potentials
Molten salts are promising candidates in numerous clean energy applications, where knowledge of thermophysical properties and vapor pressure across their operating temperature ranges is critical for safe operations. Due to challenges in evaluating these properties using experimental methods, fast and scalable molecular simulations are essential to complement the experimental data. In this study, we developed machine learning interatomic potentials (MLIP) to study the AlCl 3 molten salt across varied thermodynamic conditions (T = 473–613 K and P = 2.7–23.4 bar), which allowed us to predict temperature-surface tension correlations and liquid–vapor phase diagram from direct simulations of two-phase coexistence in this molten salt. Two MLIP architectures, a Kernel-based potential and neural network interatomic potential (NNIP), were considered to benchmark their performance for AlCl 3 molten salt using experimental structure and density values. The NNIP potential employed in two-phase equilibrium simulations yields the critical temperature and critical density of AlCl 3 that are within 10 K (∼3%) and 0.03 g/cm 3 (∼7%) of the reported experimental values. An accurate correlation between temperature and viscosities is obtained as well. In doing so, we report that the inclusion of low-density configurations in their training is critical to more accurately represent the AlCl 3 system across a wide phase-space. The MLIP trained using PBE-D3 functional in the ab initio molecular dynamics (AIMD) simulations (120 atoms) also showed close agreement with experimentally determined molten salt structure comprising Al 2 Cl 6 dimers, as validated using Raman spectra and neutron structure factor. Furthermore, the PBE-D3 as well as its trained MLIP showed better liquid density and temperature correlation for AlCl 3 system when compared to several other density functionals explored in this work. Overall, the demonstrated approach to predict temperature correlations for liquid and vapor densities in this study can be employed to screen nuclear reactors-relevant compositions, helping to mitigate safety concerns.
Ab initio molecular dynamics↗
Solution-processable and functionalizable ultra-high molecular weight polymers via topochemical synthesis
Abstract Topochemical polymerization reactions hold the promise of producing ultra-high molecular weight crystalline polymers. However, the totality of topochemical polymerization reactions has failed to produce ultra-high molecular weight polymers that are both soluble and display variable functionality, which are restrained by the crystal-packing and reactivity requirements on their respective monomers in the solid state. Herein, we demonstrate the topochemical polymerization reaction of a family of para -azaquinodimethane compounds that undergo facile visible light and thermally initiated polymerization in the solid state, allowing for the first determination of a topochemical polymer crystal structure resolved via the cryoelectron microscopy technique of microcrystal electron diffraction. The topochemical polymerization reaction also displays excellent functional group tolerance, accommodating both solubilizing side chains and reactive groups that allow for post-polymerization functionalization. The thus-produced soluble ultra-high molecular weight polymers display superior capacitive energy storage properties. This study overcomes several synthetic and characterization challenges amongst topochemical polymerization reactions, representing a critical step toward their broader application.
37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗