Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “consistency”

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 541 records · Page 30
All sourcesNASA NTRSDOE OSTI
DOE OSTI2013

A self-consistent field study of diblock copolymer/charged particle system morphologies for nanofiltration membranes

A combination of self-consistent field theory and density functional theory was used to examine the stable, 3-dimensional equilibrium morphologies formed by diblock copolymers with a tethered nanoparticle attached either between the two blocks or at the end of one of the blocks. Both neutral and interacting particles were examined, with and without favorable/unfavorable energetic potentials between the particles and the block segments. The phase diagrams of the various systems were constructed, allowing the identification of three types of ordered mesophases composed of lamellae, hexagonally packed cylinders, and spheroids. In particular, we examined the conditions under which the mesophases could be generated wherein the tethered particles were primarily located within the interface between the two blocks of the copolymer. Key factors influencing these properties were determined to be the particle position along the diblock chain, the interaction potentials of the blocks and particles, the block copolymer composition, and molecular weight of the copolymer.

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

Cubic-scaling algorithm and self-consistent field for the random-phase approximation with second-order screened exchange

The random-phase approximation with second-order screened exchange (RPA+SOSEX) is a model of electron correlation energy with two caveats: its accuracy depends on an arbitrary choice of mean field, and it scales as O(n{sup 5}) operations and O(n{sup 3}) memory for n electrons. We derive a new algorithm that reduces its scaling to O(n{sup 3}) operations and O(n{sup 2}) memory using controlled approximations and a new self-consistent field that approximates Brueckner coupled-cluster doubles theory with RPA+SOSEX, referred to as Brueckner RPA theory. The algorithm comparably reduces the scaling of second-order Møller-Plesset perturbation theory with smaller cost prefactors than RPA+SOSEX. Within a semiempirical model, we study H{sub 2} dissociation to test accuracy and H{sub n} rings to verify scaling.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2014

First principles molecular dynamics without self-consistent field optimization

We present a first principles molecular dynamics approach that is based on time-reversible extended Lagrangian Born-Oppenheimer molecular dynamics [A. M. N. Niklasson, Phys. Rev. Lett. 100, 123004 (2008)] in the limit of vanishing self-consistent field optimization. The optimization-free dynamics keeps the computational cost to a minimum and typically provides molecular trajectories that closely follow the exact Born-Oppenheimer potential energy surface. Only one single diagonalization and Hamiltonian (or Fockian) construction are required in each integration time step. The proposed dynamics is derived for a general free-energy potential surface valid at finite electronic temperatures within hybrid density functional theory. Even in the event of irregular functional behavior that may cause a dynamical instability, the optimization-free limit represents a natural starting guess for force calculations that may require a more elaborate iterative electronic ground state optimization. Our optimization-free dynamics thus represents a flexible theoretical framework for a broad and general class of ab initio molecular dynamics simulations.

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

Linking lipid architecture to bilayer structure and mechanics using self-consistent field modelling

To understand how lipid architecture determines the lipid bilayer structure and its mechanics, we implement a molecularly detailed model that uses the self-consistent field theory. This numerical model accurately predicts parameters such as Helfrichs mean and Gaussian bending modulus k{sub c} and k{sup ¯} and the preferred monolayer curvature J{sub 0}{sup m}, and also delivers structural membrane properties like the core thickness, and head group position and orientation. We studied how these mechanical parameters vary with system variations, such as lipid tail length, membrane composition, and those parameters that control the lipid tail and head group solvent quality. For the membrane composition, negatively charged phosphatidylglycerol (PG) or zwitterionic, phosphatidylcholine (PC), and -ethanolamine (PE) lipids were used. In line with experimental findings, we find that the values of k{sub c} and the area compression modulus k{sub A} are always positive. They respond similarly to parameters that affect the core thickness, but differently to parameters that affect the head group properties. We found that the trends for k{sup ¯} and J{sub 0}{sup m} can be rationalised by the concept of Israelachivili's surfactant packing parameter, and that both k{sup ¯} and J{sub 0}{sup m} change sign with relevant parameter changes. Although typically k{sup ¯}<0, membranes can form stable cubic phases when the Gaussian bending modulus becomes positive, which occurs with membranes composed of PC lipids with long tails. Similarly, negative monolayer curvatures appear when a small head group such as PE is combined with long lipid tails, which hints towards the stability of inverse hexagonal phases at the cost of the bilayer topology. To prevent the destabilisation of bilayers, PG lipids can be mixed into these PC or PE lipid membranes. Progressive loading of bilayers with PG lipids lead to highly charged membranes, resulting in J{sub 0}{sup m}≫0, especially at low ionic strengths. We anticipate that these changes lead to unstable membranes as these become vulnerable to pore formation or disintegration into lipid disks.

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

Time-dependent restricted-active-space self-consistent-field singles method for many-electron dynamics

The time-dependent restricted-active-space self-consistent-field singles (TD-RASSCF-S) method is presented for investigating TD many-electron dynamics in atoms and molecules. Adopting the SCF notion from the muticonfigurational TD Hartree-Fock (MCTDHF) method and the RAS scheme (single-orbital excitation concept) from the TD configuration-interaction singles (TDCIS) method, the TD-RASSCF-S method can be regarded as a hybrid of them. We prove that, for closed-shell N{sub e}-electron systems, the TD-RASSCF-S wave function can be fully converged using only N{sub e}/2 + 1 ⩽ M ⩽ N{sub e} spatial orbitals. Importantly, based on the TD variational principle, the converged TD-RASSCF-S wave function with M = N{sub e} is more accurate than the TDCIS wave function. The accuracy of the TD-RASSCF-S approach over the TDCIS is illustrated by the calculation of high-order harmonic generation spectra for one-dimensional models of atomic helium, beryllium, and carbon in an intense laser pulse. The electronic dynamics during the process is investigated by analyzing the behavior of electron density and orbitals. The TD-RASSCF-S method is accurate, numerically tractable, and applicable for large systems beyond the capability of the MCTDHF method.

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

Analytical gradients of complete active space self-consistent field energies using Cholesky decomposition: Geometry optimization and spin-state energetics of a ruthenium nitrosyl complex

We present a formulation of analytical energy gradients at the complete active space self-consistent field (CASSCF) level of theory employing density fitting (DF) techniques to enable efficient geometry optimizations of large systems. As an example, the ground and lowest triplet state geometries of a ruthenium nitrosyl complex are computed at the DF-CASSCF level of theory and compared with structures obtained from density functional theory (DFT) using the B3LYP, BP86, and M06L functionals. The average deviation of all bond lengths compared to the crystal structure is 0.042 Å at the DF-CASSCF level of theory, which is slightly larger but still comparable with the deviations obtained by the tested DFT functionals, e.g., 0.032 Å with M06L. Specifically, the root-mean-square deviation between the DF-CASSCF and best DFT coordinates, delivered by BP86, is only 0.08 Å for S{sub 0} and 0.11 Å for T{sub 1}, indicating that the geometries are very similar. While keeping the mean energy gradient errors below 0.25%, the DF technique results in a 13-fold speedup compared to the conventional CASSCF geometry optimization algorithm. Additionally, we assess the singlet-triplet energy vertical and adiabatic differences with multiconfigurational second-order perturbation theory (CASPT2) using the DF-CASSCF and DFT optimized geometries. It is found that the vertical CASPT2 energies are relatively similar regardless of the geometry employed whereas the adiabatic singlet-triplet gaps are more sensitive to the chosen triplet geometry.

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

Structural dynamics of phenylisothiocyanate in the light-absorbing excited states: Resonance Raman and complete active space self-consistent field calculation study

The excited state structural dynamics of phenyl isothiocyanate (PITC) after excitation to the light absorbing S{sub 2}(A′), S{sub 6}(A′), and S{sub 7}(A′) excited states were studied by using the resonance Raman spectroscopy and complete active space self-consistent field method calculations. The UV absorption bands of PITC were assigned. The vibrational assignments were done on the basis of the Fourier transform (FT)-Raman and FT-infrared measurements, the density-functional theory computations, and the normal mode analysis. The A-, B-, and C-bands resonance Raman spectra in cyclohexane, acetonitrile, and methanol solvents were, respectively, obtained at 299.1, 282.4, 266.0, 252.7, 228.7, 217.8, and 208.8 nm excitation wavelengths to probe the corresponding structural dynamics of PITC. The results indicated that the structural dynamics in the S{sub 2}(A′), S{sub 6}(A′), and S{sub 7}(A′) excited states were very different. The conical intersection point CI(S{sub 2}/S{sub 1}) were predicted to play important role in the low-lying excited state decay dynamics. Two major decay channels were predicted for PITC upon excitation to the S{sub 2}(A′) state: the radiative S{sub 2,min} → S{sub 0} transition and the nonradiative S{sub 2} → S{sub 1} internal conversion via CI(S{sub 2}/S{sub 1}). The differences in the decay dynamics between methyl isothiocyanate and PITC in the first light absorbing excited state were discussed. The role of the intersystem crossing point ISC(S{sub 1}/T{sub 1}) in the excited state decay dynamics of PITC is evaluated.

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

On the full exploitation of symmetry in periodic (as well as molecular) self-consistent-field ab initio calculations

Use of symmetry can dramatically reduce the computational cost (running time and memory allocation) of self-consistent-field ab initio calculations for molecular and crystalline systems. Crucial for running time is symmetry exploitation in the evaluation of one- and two-electron integrals, diagonalization of the Fock matrix at selected points in reciprocal space, reconstruction of the density matrix. As regards memory allocation, full square matrices (overlap, Fock, and density) in the Atomic Orbital (AO) basis are avoided and a direct transformation from the packed AO to the symmetry adapted crystalline orbital basis is performed, so that the largest matrix to be handled has the size of the largest sub-block in the latter basis. Quantitative examples, referring to the implementation in the CRYSTAL code, are given for high symmetry families of compounds such as carbon fullerenes and nanotubes.

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

Self-consistent field theory of tethered polymers: One dimensional, three dimensional, strong stretching theories and the effects of excluded-volume-only interactions

We examine end-tethered polymers in good solvents, using one- and three-dimensional self-consistent field theory, and strong stretching theories. We also discuss different tethering scenarios, namely, mobile tethers, fixed but random ones, and fixed but ordered ones, and the effects and important limitations of including only binary interactions (excluded volume terms). We find that there is a “mushroom” regime in which the layer thickness is independent of the tethering density, σ, for systems with ordered tethers, but we argue that there is no such plateau for mobile or disordered anchors, nor is there one in the 1D theory. In the other limit of brushes, all approaches predict that the layer thickness scales linearly with N. However, the σ{sup 1/3} scaling is a result of keeping only excluded volume interactions: when the full potential is included, the dependence is faster and more complicated than σ{sup 1/3}. In fact, there does not appear to be any regime in which the layer thickness scales in the combination Nσ{sup 1/3}. We also compare the results for two different solvents with each other, and with earlier Θ solvent results.

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

Investigation of electron-atom/molecule scattering resonances: Two complex multiconfigurational self-consistent field approaches

Resonances are temporarily bound states which lie in the continuum part of the Hamiltonian. If the electronic coordinates of the Hamiltonian are scaled (“dilated”) by a complex parameter, η = αe{sup iθ} (α, θ real), then its complex eigenvalues represent the scattering states (resonant and non-resonant) while the eigenvalues corresponding to the bound states and the ionization and the excitation thresholds remain real and unmodified. These make the study of these transient species amenable to the bound state methods. We developed a quadratically convergent multiconfigurational self-consistent field method (MCSCF), a well-established bound-state technique, combined with a dilated Hamiltonian to investigate resonances. This is made possible by the adoption of a second quantization algebra suitable for a set of “complex conjugate biorthonormal” spin orbitals and a modified step-length constraining algorithm to control the walk on the complex energy hypersurface while searching for the stationary point using a multidimensional Newton-Raphson scheme. We present our computational results for the {sup 2}PBe{sup −} shape resonances using two different computationally efficient methods that utilize complex scaled MCSCF (i.e., CMCSCF). These two methods are to straightforwardly use CMCSCF energy differences and to obtain energy differences using an approximation to the complex multiconfigurational electron propagator. It is found that, differing from previous computational studies by others, there are actually two {sup 2}PBe{sup −} shape resonances very close in energy. In addition, N{sub 2} resonances are examined using one of these methods.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2015

A reduced-scaling density matrix-based method for the computation of the vibrational Hessian matrix at the self-consistent field level

An analytical method to calculate the molecular vibrational Hessian matrix at the self-consistent field level is presented. By analysis of the multipole expansions of the relevant derivatives of Coulomb-type two-electron integral contractions, we show that the effect of the perturbation on the electronic structure due to the displacement of nuclei decays at least as r{sup −2} instead of r{sup −1}. The perturbation is asymptotically local, and the computation of the Hessian matrix can, in principle, be performed with O(N) complexity. Our implementation exhibits linear scaling in all time-determining steps, with some rapid but quadratic-complexity steps remaining. Sample calculations illustrate linear or near-linear scaling in the construction of the complete nuclear Hessian matrix for sparse systems. For more demanding systems, scaling is still considerably sub-quadratic to quadratic, depending on the density of the underlying electronic structure.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2015

A hybrid framework of first principles molecular orbital calculations and a three-dimensional integral equation theory for molecular liquids: Multi-center molecular Ornstein–Zernike self-consistent field approach

In this study, we reported the development of a new quantum mechanics/molecular mechanics (QM/MM)-type framework to describe chemical processes in solution by combining standard molecular-orbital calculations with a three-dimensional formalism of integral equation theory for molecular liquids (multi-center molecular Ornstein–Zernike (MC-MOZ) method). The theoretical procedure is very similar to the 3D-reference interaction site model self-consistent field (RISM-SCF) approach. Since the MC-MOZ method is highly parallelized for computation, the present approach has the potential to be one of the most efficient procedures to treat chemical processes in solution. Benchmark tests to check the validity of this approach were performed for two solute (solute water and formaldehyde) systems and a simple S{sub N}2 reaction (Cl{sup −} + CH{sub 3}Cl → ClCH{sub 3} + Cl{sup −}) in aqueous solution. The results for solute molecular properties and solvation structures obtained by the present approach were in reasonable agreement with those obtained by other hybrid frameworks and experiments. In particular, the results of the proposed approach are in excellent agreements with those of 3D-RISM-SCF.

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

Analytical gradients of the state-average complete active space self-consistent field method with density fitting

An efficient implementation of the state-averaged complete active space self-consistent field (SA-CASSCF) gradients employing density fitting (DF) is presented. The DF allows a reduction both in scaling and prefactors of the different steps involved. The performance of the algorithm is demonstrated on a set of molecules ranging up to an iron-Heme b complex which with its 79 atoms and 811 basis functions is to our knowledge the largest SA-CASSCF gradient computed. For smaller systems where the conventional code could still be used as a reference, both the linear response calculation and the gradient formation showed a clear timing reduction and the overall cost of a geometry optimization is typically reduced by more than one order of magnitude while the accuracy loss is negligible.

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

Numerical self-consistent field theory study of the response of strong polyelectrolyte brushes to external electric fields

The response of strong polyelectrolyte (PE) brushes grafted on an electrode to electric fields generated by opposite surface charges on the PE-grafted electrode and a second parallel electrode has been numerically investigated by self-consistent field theory. The influences of grafting density, average charge fraction, salt concentration, and mobile ion size on the variation of the brush height against an applied voltage bias were investigated. In agreement with molecular dynamics simulation results, a higher grafting density requires a larger magnitude of voltage bias to achieve the same amount of relative change in the brush height. In the experimentally relevant parameter regime of the applied voltage, the brush height becomes insensitive to the voltage bias when the grafting density is high. Including the contribution of surface charges on the grafting electrode, overall charge neutrality inside the PE brushes is generally maintained, especially for PE brushes with high grafting density and high average charge fraction. Our numerical study further reveals that the electric field across the two electrodes is highly non-uniform because of the complex interplay between the surface charges on the electrodes, the charges on the grafted PE chains, and counterions.

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

Onsets of hierarchy truncation and self–consistent Born approximation with quantum mechanics prescriptions invariance

The issue of efficient hierarchy truncation is related to many approximate theories. In this paper, we revisit this issue from both the numerical efficiency and quantum mechanics prescription invariance aspects. The latter requires that the truncation approximation made in Schrödinger picture, such as the quantum master equations and their self–consistent–Born–approximation improvements, should be transferable to their Heisenberg–picture correspondences, without further approximations. We address this issue with the dissipaton equation of motion (DEOM), which is a unique theory for the dynamics of not only reduced systems but also hybrid bath environments. We also highlight the DEOM theory is not only about how its dynamical variables evolve in time, but also the underlying dissipaton algebra. We demonstrate this unique feature of DEOM with model systems and report some intriguing nonlinear Fano interferences characteristics that are experimentally measurable.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2015

Communication: An efficient approach to compute state-specific nuclear gradients for a generic state-averaged multi-configuration self consistent field wavefunction

We present a new, very efficient semi-numerical approach for the computation of state-specific nuclear gradients of a generic state-averaged multi-configuration self consistent field wavefunction. Our approach eliminates the costly coupled-perturbed multi-configuration Hartree-Fock step as well as the associated integral transformation stage. The details of the implementation within the Firefly quantum chemistry package are discussed and several sample applications are given. The new approach is routinely applicable to geometry optimization of molecular systems with 1000+ basis functions using a standalone multi-core workstation.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2016

Time-dependent density functional theory (TD-DFT) coupled with reference interaction site model self-consistent field explicitly including spatial electron density distribution (RISM-SCF-SEDD)

Theoretical approach to design bright bio-imaging molecules is one of the most progressing ones. However, because of the system size and computational accuracy, the number of theoretical studies is limited to our knowledge. To overcome the difficulties, we developed a new method based on reference interaction site model self-consistent field explicitly including spatial electron density distribution and time-dependent density functional theory. We applied it to the calculation of indole and 5-cyanoindole at ground and excited states in gas and solution phases. The changes in the optimized geometries were clearly explained with resonance structures and the Stokes shift was correctly reproduced.

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

Self-interaction-corrected Kohn–Sham effective potentials using the density-consistent effective potential method

Density functional theory (DFT) and beyond-DFT methods are often used in combination with photoelectron spectroscopy to obtain physical insights into the electronic structure of molecules and solids. The Kohn–Sham eigenvalues are not electron removal energies except for the highest occupied orbital. The eigenvalues of the highest occupied molecular orbitals often underestimate the electron removal or ionization energies due to the self-interaction (SI) errors in approximate density functionals. Here, we adapt and implement the density-consistent effective potential method of Kohut, Ryabinkin, and Staroverov (2014) to obtain SI-corrected local effective potentials from the SI-corrected Fermi–Löwdin orbitals and density in the Fermi–Löwdin orbital self-interaction correction scheme. The implementation is used to obtain the density of states (photoelectron spectra) and HOMO–LUMO gaps for a set of molecules and polyacenes. Good agreement with experimental values is obtained compared to a range of SI uncorrected density functional approximations.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗