Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “complex system theory”

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 307 records · Page 17

Repartitioned Brillouin-Wigner perturbation theory with a size-consistent second-order correlation energy

Second-order Møller-Plesset perturbation theory (MP2) often breaks down catastrophically in small-gap systems, leaving much to be desired in its performance for myriad chemical applications such as noncovalent interactions, thermochemistry, and dative bonding in transition metal complexes. This divergence problem has reignited interest in Brillouin-Wigner perturbation theory (BWPT), which is regular at all orders but lacks size consistency and extensivity, severely limiting its application to chemistry. In this work, we propose an alternative partitioning of the Hamiltonian that leads to a regular BWPT perturbation series that, through the second order, is size-extensive, size-consistent (provided its Hartree–Fock reference is also), and orbital invariant. Here, our second-order size-consistent Brillouin-Wigner (BW-s2) approach can describe the exact dissociation limit of H 2 in a minimal basis set, regardless of the spin polarization of the reference orbitals. More broadly, we find that BW-s2 offers improvements relative to MP2 for covalent bond breaking, noncovalent interaction energies, and metal/organic reaction energies, although rivaling coupled-cluster with single and double substitutions for thermochemical properties.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Guest-Host Interactions in Clathrate Hydrates: Benchmark MP2 and CCSD(T)/CBS Binding Energies of CH4, CO2 and H2S in (H2O)20 Cages

We present benchmark binding energies of naturally occurring gas molecules CH4, CO2, and H2S in the small cage, namely the pentagonal dodecahedron (512) (H2O)20, which is one of the constituent cages of the 3 major lattices (structures I, II and H) of clathrate hydrates. These weak interactions require higher levels of electron correlation and converge slowly with increasing basis set to the Complete Basis Set (CBS) limit, necessitating the use of large basis sets up to the augcc- pV5Z and subsequent correction for Basis Set Superposition Error (BSSE). For the host hollow (H2O)20 cages we have identified a most stable isomer with binding energy of -200.8 ± 2.1 kcal/mol at the CCSD(T)/CBS limit (-199.2 ± 0.5 kcal/mol at the MP2/CBS limit). Additionally, we report converged second order Moller-Plesset (MP2) CBS binding energies for the encapsulation of guests in the (H2O)20 cage of -4.3 ± 0.1 for CH4@(H2O)20, -6.6 ± 0.1 for CO2@(H2O)20 and -8.5 ± 0.1 kcal/mol for H2S@(H2O)20, respectively. For CH4@(H2O)20, exhibiting the weakest encapsulation affinity among the three, we report CCSD(T)/aug-cc-pVTZ binding energies and, based on them, a CCSD(T)/CBS estimate of -4.75 ± 0.1 kcal/mol. To the best of our knowledge, the CCSD(T)/aug-cc-pVTZ calculation for CH4@(H2O)20 is the largest one reported to date (168 valence electrons, 1978 basis functions and the correlation of 84 doubly occupied and 1873 virtual orbitals) and required a scalable implementation of the (T) module on 6144 nodes (350208 cores) of the “Cori” supercomputer at the National Energy Research Supercomputing Center (NERSC) for a total execution time of 195 minutes (for the (T) part). These efficient scalable implementations of highly correlated methods offer the capability to obtain long-lasting benchmarks of intermolecular interactions in complex systems. They also provide a path towards parametrizing classical potentials needed to study the dynamical and transport properties in these complex systems as well as assess the accuracy of lower scaling electronic structure methods such as Density Functional Theory (DFT) and MP2 including its spin-biased variants.

Heindel, Joseph↗

Towards complete assignment of the infrared spectrum of the protonated water cluster H + (H 2 O) 21

The spectroscopic features of protonated water species in dilute acid solutions have been long sought after for understanding the microscopic behavior of the proton in water with gas-phase water clusters H + (H 2 O) n extensively studied as bottom-up model systems. We present a new protocol for the calculation of the infrared (IR) spectra of complex systems, which combines the fragment-based Coupled Cluster method and anharmonic vibrational quasi-degenerate perturbation theory, and demonstrate its accuracy towards the complete and accurate assignment of the IR spectrum of the H + (H 2 O) 21 cluster. The site-specific IR spectral signatures reveal two distinct structures for the internal and surface four-coordinated water molecules, which are ice-like and liquid-like, respectively. The effect of inter-molecular interaction between water molecules is addressed, and the vibrational resonance is found between the O-H stretching fundamental and the bending overtone of the nearest neighboring water molecule. The revelation of the spectral signature of the excess proton offers deeper insight into the nature of charge accommodation in the extended hydrogen-bonding network underpinning this aqueous cluster.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Self-consistent continuum solvation for optical absorption of complex molecular systems in solution

We introduce a new method to compute the optical absorption spectra of complex molecular systems in solution, based on the Liouville approach to time-dependent density-functional perturbation theory and the revised self-consistent continuum solvation model. The former allows one to obtain the absorption spectrum over a whole wide frequency range, using a recently proposed Lanczos-based technique, or selected excitation energies, using the Casida equation, without having to ever compute any unoccupied molecular orbitals. The latter is conceptually similar to the polarizable continuum model and offers the further advantages of allowing an easy computation of atomic forces via the Hellmann-Feynman theorem and a ready implementation in periodic-boundary conditions. The new method has been implemented using pseudopotentials and plane-wave basis sets, benchmarked against polarizable continuum model calculations on 4-aminophthalimide, alizarin, and cyanin and made available through the QUANTUM ESPRESSO distribution of open-source codes.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Systematic approach for simultaneously correcting the band-gap and p - d separation errors of common cation III-V or II-VI binaries in density functional theory calculations within a local density approximation

We propose a systematic approach that can empirically correct three major errors typically found in a density functional theory (DFT) calculation within the local density approximation (LDA) simultaneously for a set of common cation binary semiconductors, such as III-V compounds, (Ga or In)X with X = N,P,As,Sb, and II-VI compounds, (Zn or Cd)X, with X = O,S,Se,Te. By correcting (1) the binary band gaps at high-symmetry points , L, X, (2) the separation of p-and d-orbital-derived valence bands, and (3) conduction band effective masses to experimental values and doing so simultaneously for common cation binaries, the resulting DFT-LDA-based quasi-first-principles method can be used to predict the electronic structure of complex materials involving multiple binaries with comparable accuracy but much less computational cost than a GW level theory. This approach provides an efficient way to evaluate the electronic structures and other material properties of complex systems, much needed for material discovery and design.

Wang, Jianwei↗

Complex magnetic order in nickelate slabs

Here, magnetic ordering phenomena have a profound influence on the macroscopic properties of correlated-electron materials, but their realistic prediction remains a formidable challenge. An archetypical example is the ternary nickel oxide system RNiO 3 (R = rare earth), where the period-four magnetic order with proposals of collinear and non-collinear structures and the amplitude of magnetic moments on different Ni sublattices have been subjects of debate for decades. Here we introduce an elementary model system—NdNiO 3 slabs embedded in a non-magnetic NdGaO 3 matrix—and use polarized resonant X-ray scattering (RXS) to show that both collinear and non-collinear magnetic structures can be realized, depending on the slab thickness. The crossover between both spin structures is correctly predicted by density functional theory and can be qualitatively understood in a low-energy spin model. We further demonstrate that the amplitude ratio of magnetic moments in neighbouring NiO 6 octahedra can be accurately determined by RXS in combination with a correlated double cluster model. Targeted synthesis of model systems with controlled thickness and synergistic application of polarized RXS and ab initio theory thus provide new perspectives for research on complex magnetism, in analogy to two-dimensional materials created by exfoliation.

25 ENERGY STORAGE↗

Binding of carboxylate and water to monovalent cations

The interactions of carboxylate anions with water and cations are important for a wide variety of systems, both biological and synthetic. Here, in order to gain insight on properties of the local complexes, we apply density functional theory, to treat the complex electrostatic interactions, and investigate mixtures with varied numbers of carboxylate anions (acetate) and waters binding to monovalent cations, Li + , Na + and K + . The optimal structure with overall lowest free energy contains two acetates and two waters such that the cation is four-fold coordinated, similar to structures found earlier for pure water or pure carboxylate ligands. More generally, the complexes with two acetates have the lowest free energy. In transitioning from the overall optimal state, exchanging an acetate for water has a lower free energy barrier than exchanging water for an acetate. In most cases, the carboxylates are monodentate and in the first solvation shell. As water is added to the system, hydrogen bonding between waters and carboxylate O atoms further stabilizes monodentate structures. These structures, which have strong electrostatic interactions that involve hydrogen bonds of varying strength, are significantly polarized, with ChelpG partial charges that vary substantially as the bonding geometry varies. Overall, these results emphasize the increasing importance of water as a component of binding sites as the number of ligands increases, thus affecting the preferential solvation of specific metal ions and clarifying Hofmeister effects. Finally, structural analysis correlated with free energy analysis supports the idea that binding to more than the preferred number of carboxylates under architectural constraints are a key to ion transport.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Chemically controlled shape-morphing of elastic sheets

Two-dimensional responsive materials that change shape into complex three-dimensional structures are valuable for creating systems ranging from wearable electronics to soft robotics. Typically, the final 3D structure is unique and predetermined through the materials’ processing. Here, we use theory and simulation to devise a distinctive approach for driving shape changes of 2D elastic sheets in fluid-filled microchambers. The sheets are coated with catalyst to generate controllable fluid flows, which transform the sheets into complex 3D shapes. A given shape can be achieved by patterning the arrangement of the catalytic domains on the sheet and introducing the appropriate reactant to initiate a specific catalytic reaction. Moreover, a single sheet that encompasses multiple catalytic domains can be transformed into a variety of 3D shapes through the addition of one or more reactants. Materials systems that morph on-demand into a variety of distinct structures can simplify manufacturing processes and broaden the utility of soft materials.

36 MATERIALS SCIENCE↗

Emulating a System Dynamics Model with Agent-Based Models: A Methodological Case Study in Simulation of Diabetes Progression

An agent-based simulation model hierarchy emulating disease states and behaviors critical to progression of diabetes type 2 was designed and implemented in the DEVS framework. The models are translations of basic elements of an established system dynamics model of diabetes. In this model hierarchy, which mimics diabetes progression over an aggregated U.S. population, was dis-aggregated and reconstructed bottom-up at the individual (agent) level. Four levels of model complexity were defined in order to systematically evaluate which parameters are needed to mimic outputs of the system dynamics model. Moreover, the four estimated models attempted to replicate stock counts representing disease states in the system dynamics model, while estimating impacts of an elderliness factor, obesity factor and health-related behavioral parameters. Health-related behavior was modeled as a simple realization of the Theory of Planned Behavior, a joint function of individual attitude and diffusion of social norms that spread over each agent s social network. Although the most complex agent-based simulation model contained 31 adjustable parameters, all models were considerably less complex than the system dynamics model which required numerous time series inputs to make its predictions. In all three elaborations of the baseline model provided significantly improved fits to the output of the system dynamics model. The performances of the baseline agent-based model and its extensions illustrate a promising approach to translate complex system dynamics models into agent-based model alternatives that are both conceptually simpler and capable of capturing main effects of complex local agent-agent interactions.

59 BASIC BIOLOGICAL SCIENCES↗

Kinetics of the Strongly Correlated CH 3 O + O 2 Reaction: The Importance of Quadruple Excitations in Atmospheric and Combustion Chemistry

Kinetics measurements on radical–radical reactions are often unavailable experimentally, and obtaining quantitative rate constants for such reactions by theoretical methods is challenging because the transition states and the reactants are often strongly correlated. Treating strongly correlated systems by coupled cluster theory limited to single, double, and triple connected excitations is often inadequate. We therefore use a new method, called GMM(P), for extrapolation to the complete configuration interaction limit to go beyond triple excitations and in particular to approximate the CCSDTQ(P)/CBS limit. Here, we present this method and use it to investigate the CH 3 O + O 2 reaction. The contribution of connected quadruple excitations to the barrier height energy is found to be -3.13 kcal/mol, and adding a quasiperturbative calculation of the effect of connected pentuple excitations brings the post-connected-triples contributions to -3.44 kcal/mol, which corresponds to Boltzmann factors that increase calculated rate constants by factors of 1.0 × 10 3 , 3.3 × 10 2 , and 18 at 250, 298, and 600 K, respectively. We present rate constants for temperatures from 250 to 2000 K, and we find that the Arrhenius activation energy increases from 0.58 to 9.68 kcal/mol over this range. We also find reasonably good accuracy for the barrier height with the MN15-L exchange–correlation functional, and we calculate rate constants by a combination of GMM(P) and MN15-L electronic structure calculations and conventional and variational transition state theory, in particular canonical variational theory with small-curvature tunneling. The present findings have broad implications for obtaining quantitative rate constants for complex reaction systems in atmospheric and combustion chemistry.

Chemistry↗

Medial packing and elastic asymmetry stabilize the double-gyroid in block copolymers

Triply-periodic networks are among the most complex and functionally valuable self-assembled morphologies, yet they form in nearly every class of biological and synthetic soft matter building blocks. In contrast to simpler assembly motifs – spheres, cylinders, layers – networks require molecules to occupy variable local environments, confounding attempts to understand their formation. Here, we examine the double-gyroid network phase by using a geometric formulation of the strong stretching theory of block copolymer melts, a prototypical soft self-assembly system. The theory establishes the direct link between molecular packing, assembly thermodynamics and the medial map, a generic measure of the geometric center of complex shapes. We show that “medial packing” is essential for stability of double-gyroid in strongly-segregated melts, reconciling a long-standing contradiction between infinite- and finite-segregation theories. Additionally, we find a previously unrecognized non-monotonic dependence of network stability on the relative entropic elastic stiffness of matrix-forming to tubular-network forming blocks. The composition window of stable double-gyroid widens for both large and small elastic asymmetry, contradicting intuitive notions that packing frustration is localized to the tubular domains. This study demonstrates the utility of optimized medial tessellations for understanding soft-molecular assembly and packing frustration via an approach that is readily generalizable far beyond gyroids in neat block copolymers.

36 MATERIALS SCIENCE↗

Structural complexity of snapshots of two-dimensional Fermi-Hubbard systems

The development of quantum gas microscopy for two-dimensional optical lattices has provided an unparalleled tool to study the Fermi-Hubbard model (FHM) with ultracold atoms. Spin-resolved projective measurements, or snapshots, have played a significant role in quantifying correlation functions, theory verification, and thus the uncovering of underlying physical phenomena such as antiferromagnetism at commensurate filling on bipartite lattices and other charge and spin correlations, as well as dynamical properties at various densities. Here we employ a recent concept, the multiscale structural complexity, and show that when computed for the snapshots (of single spin species, local moments, or total density) it can provide a theory-free property, immediately accessible to experiments. Specifically, after benchmarking results for Ising and $XY$ models, we study the structural complexity of snapshots of the repulsive FHM in the two-dimensional square lattice as a function of doping and temperature. We generate projective measurements using determinant quantum Monte Carlo and compare their complexities against those from the experiment. We demonstrate that these complexities are linked to relevant physical observables such as the entropy and double occupancy. Their behaviors capture the development of correlations and relevant length scales in the system. Furthermore, we provide an open-source code in python which can be implemented into data analysis routines in experimental settings for the square lattice.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗

Dynamics of multiple viscoelastic carbon nanotube based nanocomposites with axial magnetic field

Nanocomposites and magnetic field effects on nanostructures have received great attention in recent years. A large amount of research work was focused on developing the proper theoretical framework for describing many physical effects appearing in structures on nanoscale level. Great step in this direction was successful application of nonlocal continuum field theory of Eringen. In the present paper, the free transverse vibration analysis is carried out for the system composed of multiple single walled carbon nanotubes (MSWCNT) embedded in a polymer matrix and under the influence of an axial magnetic field. Equivalent nonlocal model of MSWCNT is adopted as viscoelastically coupled multi-nanobeam system (MNBS) under the influence of longitudinal magnetic field. Governing equations of motion are derived using the Newton second low and nonlocal Rayleigh beam theory, which take into account small-scale effects, the effect of nanobeam angular acceleration, internal damping and Maxwell relation. Explicit expressions for complex natural frequency are derived based on the method of separation of variables and trigonometric method for the “Clamped-Chain” system. In addition, an analytical method is proposed in order to obtain asymptotic damped natural frequency and the critical damping ratio, which are independent of boundary conditions and a number of nanobeams in MNBS. The validity of obtained results is confirmed by comparing the results obtained for complex frequencies via trigonometric method with the results obtained by using numerical methods. The influence of the longitudinal magnetic field on the free vibration response of viscoelastically coupled MNBS is discussed in detail. In addition, numerical results are presented to point out the effects of the nonlocal parameter, internal damping, and parameters of viscoelastic medium on complex natural frequencies of the system. The results demonstrate the efficiency of the suggested methodology to find the closed form solutions for the free vibration response of multiple nanostructure systems under the influence of magnetic field.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

Measurement of pulsed-power-driven magnetic fields via proton deflectometry

Measuring magnetic field and current distribution in Z-pinch plasma systems is crucial to the validation of Z-pinch theory. In this letter, the demonstration of proton deflectometry to pulsed-power-driven loads at the mega-amp scale is presented, which is capable of making more detailed field maps in high-density regions of plasmas. In this method, a laser-driven, broad-spectrum, MeV-energy proton beam is directed through a pulsed-power-driven plasma system, and the resulting deflections are measured to examine configuration of magnetic fields and to infer the currents that support them. The technique was first demonstrated on simple short-circuit loads, and the results are in excellent agreement with numerical simulations providing reliable estimates of the field and current configurations. It was then applied to a more complex—radial foil—plasma load. The measurements show unexpected proton deflections that exhibit the complexity of the plasma load and that with further analysis will reveal details about the current and magnetic field topology in this complex configuration.

Mariscal, D.↗

A multi-species exchange model for fully fluctuating polymer field theory simulations

Field-theoretic models have been used extensively to study the phase behavior of inhomogeneous polymer melts and solutions, both in self-consistent mean-field calculations and in numerical simulations of the full theory capturing composition fluctuations. The models commonly used can be grouped into two categories, namely, species models and exchange models. Species models involve integrations of functionals that explicitly depend on fields originating both from species density operators and their conjugate chemical potential fields. In contrast, exchange models retain only linear combinations of the chemical potential fields. In the two-component case, development of exchange models has been instrumental in enabling stable complex Langevin (CL) simulations of the full complex-valued theory. No comparable stable CL approach has yet been established for field theories of the species type. Here, we introduce an extension of the exchange model to an arbitrary number of components, namely, the multi-species exchange (MSE) model, which greatly expands the classes of soft material systems that can be accessed by the complex Langevin simulation technique. We demonstrate the stability and accuracy of the MSE-CL sampling approach using numerical simulations of triblock and tetrablock terpolymer melts, and tetrablock quaterpolymer melts. This method should enable studies of a wide range of fluctuation phenomena in multiblock/multi-species polymer blends and composites.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

An introduction to orbit dynamics and its application to satellite-based earth monitoring systems

The long term behavior of satellites is studied at a level of complexity suitable for the initial planning phases of earth monitoring missions. First-order perturbation theory is used to describe in detail the basic orbit dynamics of satellite motion around the earth and relative to the sun. Surface coverage capabilities of satellite orbits are examined. Several examples of simulated observation and monitoring missions are given to illustrate representative applications of the theory. The examples stress the need for devising ways of maximizing total mission output in order to make the best possible use of the resultant data base as input to those large-scale, long-term earth monitoring activities which can best justify the use of satellite systems.

Brooks, D. R.↗

Orbital-Resolved Imaging of the Adsorbed State of Pyridine on GaP(110) Identifies Sites Susceptible to Nucleophilic Attack

Artificial photosynthesis by photoelectrocatalytic CO 2 reduction is dependent, as is natural photosynthesis, on interfacial electron transfer to couple light excitation energy to reaction centers. For heterogeneous systems, in the context of frontier orbital theory, artificial reaction centers are defined through the interactions of filled and empty orbitals within a few electron volts of the Fermi energy of the adsorbate complex. Here we report a scanning tunneling microscopy (STM) and density functional theory investigation of the orbital-resolved adsorption state defining the dative bonding interaction between a III-V semiconductor surface [GaP(110)] and a N containing heteroaromatic (pyridine). This system was selected for its relevance to photoelectrocatalysis utilizing heteroaromatic co-catalysts, which has been reported to yield highly selective CO 2 reduction to fuels. By examining the distribution of unoccupied molecular orbitals, we show that STM images can be used to positively identify the sites on pyridine susceptible to nucleophilic attack, consistent with frontier orbital theory. Furthermore, this indicates that STM can be used to explore the local reaction centers of adsorbed ambidentate electrophiles and nucleophiles relevant to artificial photosynthesis, and more broadly to generate critical mechanistic information for various heterogeneous acid-base reactions.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Theoretical study of chromophores for biological sensing: Understanding the mechanism of rhodol based multi-chromophoric systems

Development of two-photon fluorescent probes can aid in visualizing the cellular environment. Multi-chromophore systems display complex manifolds of electronic transitions, enabling their use for optical sensing applications. Time-Dependent Density Functional Theory (TDDFT) methods allow for accurate predictions of the optical properties. These properties are related to the electronic transitions in the molecules, which include two-photon absorption cross-sections. We use TDDFT to understand the mechanism of aza-crown based fluorescent probes for metals sensing applications. Here our findings suggest changes in local excitation in the rhodol chromophore between unbound form and when bound to the metal analyte. These changes are caused by a charge transfer from the aza-crown group and pyrazol units toward the rhodol unit. Understanding this mechanism leads to an optimized design with higher two-photon excited fluorescence to be used in medical applications.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗