Browse
Submit Data
Databases
API
Help

Dataset Information

7 Views

0 Connections

0 Citations

0 Reanalyses

0 Downloads

Omics score: 0

Equation-of-Motion MLCCSD and CCSD-in-HF Oscillator Strengths and Their Application to Core Excitations.

ABSTRACT: We present an implementation of equation-of-motion oscillator strengths for the multilevel CCSD (MLCCSD) model where CCS is used as the lower level method (CCS/CCSD). In this model, the double excitations of the cluster operator are restricted to an active orbital space, whereas the single excitations are unrestricted. Calculated nitrogen K-edge spectra of adenosine, adenosine triphosphate (ATP), and an ATP-water system are used to demonstrate the performance of the model. Projected atomic orbitals (PAOs) are used to partition the virtual space into active and inactive orbital sets. Cholesky decomposition of the Hartree-Fock density is used to partition the occupied orbitals. This Cholesky-PAO partitioning is cheap, scaling as O(N3), and is suitable for the calculation of core excitations, which are localized in character. By restricting the single excitations of the cluster operator to the active space, as well as the double excitations, the CCSD-in-HF model is obtained. A comparison of the two models-MLCCSD and CCSD-in-HF-is presented for the core excitation spectra of the adenosine and ATP systems.

SUBMITTER: Folkestad SD

PROVIDER: S-EPMC8011930 | biostudies-literature |

REPOSITORIES: biostudies-literature

ACCESS DATA

Json Xml

Similar Datasets

Multilevel CC2 and CCSD in Reduced Orbital Spaces: Electronic Excitations in Large Molecular Systems.

Project description:We present efficient implementations of the multilevel CC2 (MLCC2) and multilevel CCSD (MLCCSD) models. As the system size increases, MLCC2 and MLCCSD exhibit the scaling of the lower-level coupled cluster model. To treat large systems, we combine MLCC2 and MLCCSD with a reduced-space approach in which the multilevel coupled cluster calculation is performed in a significantly truncated molecular orbital basis. The truncation scheme is based on the selection of an active region of the molecular system and the subsequent construction of localized Hartree-Fock orbitals. These orbitals are used in the multilevel coupled cluster calculation. The electron repulsion integrals are Cholesky decomposed using a screening protocol that guarantees accuracy in the truncated molecular orbital basis and reduces computational cost. The Cholesky factors are constructed directly in the truncated basis, ensuring low storage requirements. Systems for which Hartree-Fock is too expensive can be treated by using a multilevel Hartree-Fock reference. With the reduced-space approach, we can handle systems with more than a thousand atoms. This is demonstrated for paranitroaniline in aqueous solution.

| S-EPMC8016205 | biostudies-literature

The realization of the dipole (γ, γ) method and its application to determine the absolute optical oscillator strengths of helium.

Project description:The dipole (γ, γ) method, which is the inelastic x-ray scattering operated at a negligibly small momentum transfer, is proposed and realized to determine the absolute optical oscillator strengths of the vanlence-shell excitations of atoms and molecules. Compared with the conventionally used photoabsorption method, this new method is free from the line saturation effect, which can seriously limit the accuracies of the measured photoabsorption cross sections for discrete transitions with narrow natural linewidths. Furthermore, the Bethe-Born conversion factor of the dipole (γ, γ) method varies much more slowly with the excitation energy than does that of the dipole (e, e) method. Absolute optical oscillator strengths for the excitations of 1s(2) → 1 snp(n = 3-7) of atomic helium have been determined using the high-resolution dipole (γ, γ) method, and the excellent agreement of the present measurements with both those measured by the dipole (e, e) method and the previous theoretical calculations indicates that the dipole (γ, γ) method is a powerful tool to measure the absolute optical oscillator strengths of the valence-shell excitations of atoms and molecules.

| S-EPMC4683586 | biostudies-other

Spin-Conserved and Spin-Flip Optical Excitations from the Bethe-Salpeter Equation Formalism.

Project description:Like adiabatic time-dependent density-functional theory (TD-DFT), the Bethe-Salpeter equation (BSE) formalism of many-body perturbation theory, in its static approximation, is "blind" to double (and higher) excitations, which are ubiquitous, for example, in conjugated molecules like polyenes. Here, we apply the spin-flip ansatz (which considers the lowest triplet state as the reference configuration instead of the singlet ground state) to the BSE formalism in order to access, in particular, double excitations. The present scheme is based on a spin-unrestricted version of the GW approximation employed to compute the charged excitations and screened Coulomb potential required for the BSE calculations. Dynamical corrections to the static BSE optical excitations are taken into account via an unrestricted generalization of our recently developed (renormalized) perturbative treatment. The performance of the present spin-flip BSE formalism is illustrated by computing excited-state energies of the beryllium atom, the hydrogen molecule at various bond lengths, and cyclobutadiene in its rectangular and square-planar geometries.

| S-EPMC8154368 | biostudies-literature

Benzene, Toluene, and Monosubstituted Derivatives: Diabatic Nature of the Oscillator Strengths of S1 ← S0 Transitions.

Project description:For benzene, toluene, aniline, fluorobenzene, and phenol, even sophisticated treatments of electron correlation, such as MRCI and XMS-CASPT2 calculations, show oscillator strengths typically lower than experiment. Inclusion of a simple pseudo-diabatization approach to perturb the S1 state with approximate vibronic coupling to the S2 state for each molecule results in more accurate oscillator strengths. Their absolute values agree better with experiment for all molecules except aniline. When the coupling between the S1 and S2 states is strong at the S0 geometry, the simple diabatization scheme performs less well with respect to the oscillator strengths relative to the adiabatic values. However, we expect the scheme to be useful in many cases where the coupling is weak to moderate (where the maximum component of the coupling has a magnitude less than 1.5 au). Such calculations give an insight into the effects of vibronic coupling of excited states on UV/vis spectra.

| S-EPMC8279645 | biostudies-literature

Prediction and classification in equation-free collective motion dynamics.

Project description:Modeling the complex collective behavior is a challenging issue in several material and life sciences. The collective motion has been usually modeled by simple interaction rules and explained by global statistics. However, it remains difficult to bridge the gap between the dynamic properties of the complex interaction and the emerging group-level functions. Here we introduce decomposition methods to directly extract and classify the latent global dynamics of nonlinear dynamical systems in an equation-free manner, even including complex interaction in few data dimensions. We first verified that the basic decomposition method can extract and discriminate the dynamics of a well-known rule-based fish-schooling (or bird-flocking) model. The method extracted different temporal frequency modes with spatial interaction coherence among three distinct emergent motions, whereas these wave properties in multiple spatiotemporal scales showed similar dispersion relations. Second, we extended the basic method to map high-dimensional feature space for application to actual small-dimensional systems complexly changing the interaction rules. Using group sports human data, we classified the dynamics and predicted the group objective achievement. Our methods have a potential for classifying collective motions in various domains which obey in non-trivial dominance law known as active matters.

| S-EPMC6237418 | biostudies-literature

The equation of motion for supershear frictional rupture fronts.

Project description:The rupture fronts that mediate the onset of frictional sliding may propagate at speeds below the Rayleigh wave speed or may surpass the shear wave speed and approach the longitudinal wave speed. While the conditions for the transition from sub-Rayleigh to supershear propagation have been studied extensively, little is known about what dictates supershear rupture speeds and how the interplay between the stresses that drive propagation and interface properties that resist motion affects them. By combining laboratory experiments and numerical simulations that reflect natural earthquakes, we find that supershear rupture propagation speeds can be predicted and described by a fracture mechanics-based equation of motion. This equation of motion quantitatively predicts rupture speeds, with the velocity selection dictated by the interface properties and stress. Our results reveal a critical rupture length, analogous to Griffith's length for sub-Rayleigh cracks, below which supershear propagation is impossible. Above this critical length, supershear ruptures can exist, once excited, even for extremely low preexisting stress levels. These results significantly improve our fundamental understanding of what governs the speed of supershear earthquakes, with direct and important implications for interpreting their unique supershear seismic radiation patterns.

| S-EPMC6051736 | biostudies-literature

Constrained Fourth Order Latent Differential Equation Reduces Parameter Estimation Bias for Damped Linear Oscillator Models.

Project description:Second order linear differential equations can be used as models for regulation since under a range of parameter values they can account for return to equilibrium as well as potential oscillations in regulated variables. One method that can estimate parameters of these equations from intensive time series data is the method of Latent Differential Equations (LDE). However, the LDE method can exhibit bias in its parameters if the dimension of the time delay embedding and thus the width of the convolution kernel is not chosen wisely. This article presents a simulation study showing that a constrained fourth order Latent Differential Equation (FOLDE) model for the second order system almost completely eliminates bias as long as the width of the convolution kernel is less than two thirds the period of oscillations in the data. The FOLDE model adds two degrees of freedom over the standard LDE model but significantly improves model fit.

| S-EPMC7518521 | biostudies-literature

Multidimensional x-ray spectroscopy of valence and core excitations in cysteine.

Project description:Several nonlinear spectroscopy experiments which employ broadband x-ray pulses to probe the coupling between localized core and delocalized valence excitation are simulated for the amino acid cysteine at the K-edges of oxygen and nitrogen and the K- and L-edges of sulfur. We focus on two-dimensional (2D) and 3D signals generated by two- and three-pulse stimulated x-ray Raman spectroscopy (SXRS) with frequency-dispersed probe. We show how the four-pulse x-ray signals [Formula: see text] and [Formula: see text] can give new 3D insight into the SXRS signals. The coupling between valence- and core-excited states can be visualized in three-dimensional plots, revealing the origin of the polarizability that controls the simpler pump-probe SXRS signals.

| S-EPMC3637361 | biostudies-other

Assessment of the CCSD and CCSD(T) coupled-cluster methods in calculating heats of formation for Zn complexes.

Project description:Heats of formation were calculated using coupled-cluster methods for a series of zinc complexes. The calculated values were evaluated against previously conducted computational studies using density functional methods as well as experimental values. Heats of formation for nine neutral ZnX(n) complexes [X = -Zn, -H, -O, -F2, -S, -Cl, -Cl2, -CH3, (-CH3)2] were determined at the CCSD and CCSD(T) levels using the 6-31G** and TZVP basis sets as well as the LANL2DZ-6-31G** (LACVP**) and LANL2DZ-TZVP hybrid basis sets. The CCSD(T)/6-31G** level of theory was found to predict the heat of formation for the nonalkyl Zn complexes most accurately. The alkyl Zn species were problematic in that none of the methods that were tested accurately predicted the heat of formation for these complexes. In instances where experimental geometric parameters were available, these were most accurately predicted by the CCSD/6-31G** level of theory; going to CCSD(T) did not improve agreement with the experimental values. Coupled-cluster methods did not offer a systemic improvement over DFT calculations for a given functional/basis set combination. With the exceptions of ZnH and ZnF2, there are multiple density functionals that outperform coupled-cluster calculations with the 6-31G** basis set.

| S-EPMC2903681 | biostudies-literature

Assessment of the CCSD and CCSD(T) Coupled-Cluster Methods in Calculating Heats of Formation for Cu Complexes.

Project description:Heats of formation for nine complexes of the form CuX(n) (X = Cu, H, O, OH, S, F, F(2), Cl, Cl(2)) were calculated using the CCSD and CCSD(T) coupled cluster methods with the 6-31G** and TZVP basis sets as well as the LANL2DZ basis set/pseudopotential on Cu with both the 6-31G** and TZVP basis sets applied to the nonmetal atoms. These values were compared with literature heat of formation values. A second order Douglas-Kroll-Hess relativistic correction was applied at the CCSD/TZVP and CCSD(T)/TZVP levels of theory. Overall, the CCSD(T)/TZVP level of theory with the relativistic correction was most suited for the heat of formation calculations possessing low absolute average error and RMSD and the ability to analyze each copper complex, except for the problematic case of copper(II) fluoride. Finally, experimental geometric parameters were compared with the calculated structures in such cases where these data were available. None of the investigated levels of theory predicted bond lengths consistently better than other methods, and it was determined that the most accurate bond length does not necessarily result in the most accurate calculated heat of formation value for a given complex.

| S-EPMC2885807 | biostudies-literature

OmicsDI is part of the ELIXIR infrastructure

OmicsDI is an Elixir interoperability service. Learn more ›

OmicsDI Databases

PRIDE
PeptideAtlas
MassIVE
JPOST Repository
Physiome Model Repository

EGA
EVA
ENA
LINCS
PAXDB
Cell Collective

MetaboLights
Metabolomics Workbench
MetabolomeExpress
GNPS
BioModels
FAIRDOMHub

ArrayExpress
dbGaP
ExpressionAtlas
GEO
NODE

Information

Databases
Help
API
Contact us
Code on GitHub
Terms of use
Submit Data