Two-body density matrix for closed s-d shell nuclei

International Nuclear Information System (INIS)

Dimitrova, S.S.; Kadrev, D.N.; Antonov, A.N.; Stoitsov, M.V.

2000-01-01

The two-body density matrix for 4 He, 16 O and 40 Ca within the Low-order approximation of the Jastrow correlation method is considered. Closed analytical expressions for the two-body density matrix, the center of mass and relative local densities and momentum distributions are presented. The effects of the short-range correlations on the two-body nuclear characteristics are investigated. (orig.)

Variational minimization of atomic and molecular ground-state energies via the two-particle reduced density matrix

International Nuclear Information System (INIS)

Mazziotti, David A.

2002-01-01

Atomic and molecular ground-state energies are variationally determined by constraining the two-particle reduced density matrix (2-RDM) to satisfy positivity conditions. Because each positivity condition corresponds to correcting the ground-state energies for a class of Hamiltonians with two-particle interactions, these conditions collectively provide a new approach to many-body theory that, unlike perturbation theory, can capture significantly correlated phenomena including the multireference effects of potential-energy surfaces. The D, Q, and G conditions for the 2-RDM are extended through generalized lifting operators inspired from the formal solution of N-representability. These lifted conditions agree with the hierarchy of positivity conditions presented by Mazziotti and Erdahl [Phys. Rev. A 63, 042113 (2001)]. The connection between positivity and the formal solution explains how constraining higher RDMs to be positive semidefinite improves the N representability of the 2-RDM and suggests using pieces of higher positivity conditions that computationally scale like the D condition. With the D, Q, and G conditions as well as pieces of higher positivity the electronic energies for Be, LiH, H 2 O, and BH are computed through a primal-dual interior-point algorithm for positive semidefinite programming. The variational method produces potential-energy surfaces that are highly accurate even far from the equilibrium geometry where single-reference perturbation-based methods often fail to produce realistic energies

Density matrix renormalization group with efficient dynamical electron correlation through range separation

DEFF Research Database (Denmark)

Hedegård, Erik D.; Knecht, Stefan; Kielberg, Jesper Skau

2015-01-01

We present a new hybrid multiconfigurational method based on the concept of range-separation that combines the density matrix renormalization group approach with density functional theory. This new method is designed for the simultaneous description of dynamical and static electroncorrelation...... effects in multiconfigurational electronic structure problems....

A J matrix engine for density functional theory calculations

International Nuclear Information System (INIS)

White, C.A.; Head-Gordon, M.

1996-01-01

We introduce a new method for the formation of the J matrix (Coulomb interaction matrix) within a basis of Cartesian Gaussian functions, as needed in density functional theory and Hartree endash Fock calculations. By summing the density matrix into the underlying Gaussian integral formulas, we have developed a J matrix open-quote open-quote engine close-quote close-quote which forms the exact J matrix without explicitly forming the full set of two electron integral intermediates. Several precomputable quantities have been identified, substantially reducing the number of floating point operations and memory accesses needed in a J matrix calculation. Initial timings indicate a speedup of greater than four times for the (pp parallel pp) class of integrals with speedups increasing to over ten times for (ff parallel ff) integrals. copyright 1996 American Institute of Physics

Quasi-particle entanglement: redefinition of the vacuum and reduced density matrix approach

International Nuclear Information System (INIS)

Samuelsson, P; Sukhorukov, E V; Buettiker, M

2005-01-01

A scattering approach to entanglement in mesoscopic conductors with independent fermionic quasi-particles is discussed. We focus on conductors in the tunnelling limit, where a redefinition of the quasi-particle vacuum transforms the wavefunction from a many-body product state of non-interacting particles to a state describing entangled two-particle excitations out of the new vacuum (Samuelsson, Sukhorukov and Buettiker 2003 Phys. Rev. Lett. 91 157002). The approach is illustrated with two examples: (i) a normal-superconducting system, where the transformation is made between Bogoliubov-de Gennes quasi-particles and Cooper pairs, and (ii) a normal system, where the transformation is made between electron quasi-particles and electron-hole pairs. This is compared to a scheme where an effective two-particle state is derived from the manybody scattering state by a reduced density matrix approach

Spectral function from Reduced Density Matrix Functional Theory

Science.gov (United States)

Romaniello, Pina; di Sabatino, Stefano; Berger, Jan A.; Reining, Lucia

2015-03-01

In this work we focus on the calculation of the spectral function, which determines, for example, photoemission spectra, from reduced density matrix functional theory. Starting from its definition in terms of the one-body Green's function we derive an expression for the spectral function that depends on the natural occupation numbers and on an effective energy which accounts for all the charged excitations. This effective energy depends on the two-body as well as higher-order density matrices. Various approximations to this expression are explored by using the exactly solvable Hubbard chains.

Excitation energies with linear response density matrix functional theory along the dissociation coordinate of an electron-pair bond in N-electron systems

International Nuclear Information System (INIS)

Meer, R. van; Gritsenko, O. V.; Baerends, E. J.

2014-01-01

Time dependent density matrix functional theory in its adiabatic linear response formulation delivers exact excitation energies ω α and oscillator strengths f α for two-electron systems if extended to the so-called phase including natural orbital (PINO) theory. The Löwdin-Shull expression for the energy of two-electron systems in terms of the natural orbitals and their phases affords in this case an exact phase-including natural orbital functional (PILS), which is non-primitive (contains other than just J and K integrals). In this paper, the extension of the PILS functional to N-electron systems is investigated. With the example of an elementary primitive NO functional (BBC1) it is shown that current density matrix functional theory ground state functionals, which were designed to produce decent approximations to the total energy, fail to deliver a qualitatively correct structure of the (inverse) response function, due to essential deficiencies in the reconstruction of the two-body reduced density matrix (2RDM). We now deduce essential features of an N-electron functional from a wavefunction Ansatz: The extension of the two-electron Löwdin-Shull wavefunction to the N-electron case informs about the phase information. In this paper, applications of this extended Löwdin-Shull (ELS) functional are considered for the simplest case, ELS(1): one (dissociating) two-electron bond in the field of occupied (including core) orbitals. ELS(1) produces high quality ω α (R) curves along the bond dissociation coordinate R for the molecules LiH, Li 2 , and BH with the two outer valence electrons correlated. All of these results indicate that response properties are much more sensitive to deficiencies in the reconstruction of the 2RDM than the ground state energy, since derivatives of the functional with respect to both the NOs and the occupation numbers need to be accurate

Joint refinement model for the spin resolved one-electron reduced density matrix of YTiO3 using magnetic structure factors and magnetic Compton profiles data.

Science.gov (United States)

Gueddida, Saber; Yan, Zeyin; Kibalin, Iurii; Voufack, Ariste Bolivard; Claiser, Nicolas; Souhassou, Mohamed; Lecomte, Claude; Gillon, Béatrice; Gillet, Jean-Michel

2018-04-28

In this paper, we propose a simple cluster model with limited basis sets to reproduce the unpaired electron distributions in a YTiO 3 ferromagnetic crystal. The spin-resolved one-electron-reduced density matrix is reconstructed simultaneously from theoretical magnetic structure factors and directional magnetic Compton profiles using our joint refinement algorithm. This algorithm is guided by the rescaling of basis functions and the adjustment of the spin population matrix. The resulting spin electron density in both position and momentum spaces from the joint refinement model is in agreement with theoretical and experimental results. Benefits brought from magnetic Compton profiles to the entire spin density matrix are illustrated. We studied the magnetic properties of the YTiO 3 crystal along the Ti-O 1 -Ti bonding. We found that the basis functions are mostly rescaled by means of magnetic Compton profiles, while the molecular occupation numbers are mainly modified by the magnetic structure factors.

Decay of autoionizing states in time-dependent density functional and reduced density matrix functional theory

Energy Technology Data Exchange (ETDEWEB)

Kapoor, Varun; Brics, Martins; Bauer, Dieter [Institut fuer Physik, Universitaet Rostock, 18051 Rostock (Germany)

2013-07-01

Autoionizing states are inaccessible to time-dependent density functional theory (TDDFT) using known, adiabatic Kohn-Sham (KS) potentials. We determine the exact KS potential for a numerically exactly solvable model Helium atom interacting with a laser field that is populating an autoionizing state. The exact single-particle density of the population in the autoionizing state corresponds to that of the energetically lowest quasi-stationary state in the exact KS potential. We describe how this exact potential controls the decay by a barrier whose height and width allows for the density to tunnel out and decay with the same rate as in the ab initio time-dependent Schroedinger calculation. However, devising a useful exchange-correlation potential that is capable of governing such a scenario in general and in more complex systems is hopeless. As an improvement over TDDFT, time-dependent reduced density matrix functional theory has been proposed. We are able to obtain for the above described autoionization process the exact time-dependent natural orbitals (i.e., the eigenfunctions of the exact, time-dependent one-body reduced density matrix) and study the potentials that appear in the equations of motion for the natural orbitals and the structure of the two-body density matrix expanded in them.

Comparison of Conjugate Gradient Density Matrix Search and Chebyshev Expansion Methods for Avoiding Diagonalization in Large-Scale Electronic Structure Calculations

Science.gov (United States)

Bates, Kevin R.; Daniels, Andrew D.; Scuseria, Gustavo E.

1998-01-01

We report a comparison of two linear-scaling methods which avoid the diagonalization bottleneck of traditional electronic structure algorithms. The Chebyshev expansion method (CEM) is implemented for carbon tight-binding calculations of large systems and its memory and timing requirements compared to those of our previously implemented conjugate gradient density matrix search (CG-DMS). Benchmark calculations are carried out on icosahedral fullerenes from C60 to C8640 and the linear scaling memory and CPU requirements of the CEM demonstrated. We show that the CPU requisites of the CEM and CG-DMS are similar for calculations with comparable accuracy.

Kohn–Sham exchange-correlation potentials from second-order reduced density matrices

Energy Technology Data Exchange (ETDEWEB)

Cuevas-Saavedra, Rogelio; Staroverov, Viktor N., E-mail: vstarove@uwo.ca [Department of Chemistry, The University of Western Ontario, London, Ontario N6A 5B7 (Canada); Ayers, Paul W. [Department of Chemistry and Chemical Biology, McMaster University, Hamilton, Ontario L8S 4M1 (Canada)

2015-12-28

We describe a practical algorithm for constructing the Kohn–Sham exchange-correlation potential corresponding to a given second-order reduced density matrix. Unlike conventional Kohn–Sham inversion methods in which such potentials are extracted from ground-state electron densities, the proposed technique delivers unambiguous results in finite basis sets. The approach can also be used to separate approximately the exchange and correlation potentials for a many-electron system for which the reduced density matrix is known. The algorithm is implemented for configuration-interaction wave functions and its performance is illustrated with numerical examples.

The tensor hypercontracted parametric reduced density matrix algorithm: coupled-cluster accuracy with O(r(4)) scaling.

Science.gov (United States)

Shenvi, Neil; van Aggelen, Helen; Yang, Yang; Yang, Weitao; Schwerdtfeger, Christine; Mazziotti, David

2013-08-07

Tensor hypercontraction is a method that allows the representation of a high-rank tensor as a product of lower-rank tensors. In this paper, we show how tensor hypercontraction can be applied to both the electron repulsion integral tensor and the two-particle excitation amplitudes used in the parametric 2-electron reduced density matrix (p2RDM) algorithm. Because only O(r) auxiliary functions are needed in both of these approximations, our overall algorithm can be shown to scale as O(r(4)), where r is the number of single-particle basis functions. We apply our algorithm to several small molecules, hydrogen chains, and alkanes to demonstrate its low formal scaling and practical utility. Provided we use enough auxiliary functions, we obtain accuracy similar to that of the standard p2RDM algorithm, somewhere between that of CCSD and CCSD(T).

Reduced density matrix functional theory at finite temperature

Energy Technology Data Exchange (ETDEWEB)

Baldsiefen, Tim

2012-10-15

Density functional theory (DFT) is highly successful in many fields of research. There are, however, areas in which its performance is rather limited. An important example is the description of thermodynamical variables of a quantum system in thermodynamical equilibrium. Although the finite-temperature version of DFT (FT-DFT) rests on a firm theoretical basis and is only one year younger than its brother, groundstate DFT, it has been successfully applied to only a few problems. Because FT-DFT, like DFT, is in principle exact, these shortcomings can be attributed to the difficulties of deriving valuable functionals for FT-DFT. In this thesis, we are going to present an alternative theoretical description of quantum systems in thermal equilibrium. It is based on the 1-reduced density matrix (1RDM) of the system, rather than on its density and will rather cumbersomly be called finite-temperature reduced density matrix functional theory (FT-RDMFT). Its zero-temperature counterpart (RDMFT) proved to be successful in several fields, formerly difficult to address via DFT. These fields include, for example, the calculation of dissociation energies or the calculation of the fundamental gap, also for Mott insulators. This success is mainly due to the fact that the 1RDM carries more directly accessible ''manybody'' information than the density alone, leading for example to an exact description of the kinetic energy functional. This sparks the hope that a description of thermodynamical systems employing the 1RDM via FT-RDMFT can yield an improvement over FT-DFT. Giving a short review of RDMFT and pointing out difficulties when describing spin-polarized systems initiates our work. We then lay the theoretical framework for FT-RDMFT by proving the required Hohenberg-Kohn-like theorems, investigating and determining the domain of FT-RDMFT functionals and by deriving several properties of the exact functional. Subsequently, we present a perturbative method to

Reduced density matrix functional theory at finite temperature

International Nuclear Information System (INIS)

Baldsiefen, Tim

2012-10-01

Density functional theory (DFT) is highly successful in many fields of research. There are, however, areas in which its performance is rather limited. An important example is the description of thermodynamical variables of a quantum system in thermodynamical equilibrium. Although the finite-temperature version of DFT (FT-DFT) rests on a firm theoretical basis and is only one year younger than its brother, groundstate DFT, it has been successfully applied to only a few problems. Because FT-DFT, like DFT, is in principle exact, these shortcomings can be attributed to the difficulties of deriving valuable functionals for FT-DFT. In this thesis, we are going to present an alternative theoretical description of quantum systems in thermal equilibrium. It is based on the 1-reduced density matrix (1RDM) of the system, rather than on its density and will rather cumbersomly be called finite-temperature reduced density matrix functional theory (FT-RDMFT). Its zero-temperature counterpart (RDMFT) proved to be successful in several fields, formerly difficult to address via DFT. These fields include, for example, the calculation of dissociation energies or the calculation of the fundamental gap, also for Mott insulators. This success is mainly due to the fact that the 1RDM carries more directly accessible ''manybody'' information than the density alone, leading for example to an exact description of the kinetic energy functional. This sparks the hope that a description of thermodynamical systems employing the 1RDM via FT-RDMFT can yield an improvement over FT-DFT. Giving a short review of RDMFT and pointing out difficulties when describing spin-polarized systems initiates our work. We then lay the theoretical framework for FT-RDMFT by proving the required Hohenberg-Kohn-like theorems, investigating and determining the domain of FT-RDMFT functionals and by deriving several properties of the exact functional. Subsequently, we present a perturbative method to iteratively construct

Quasi-particle energy spectra in local reduced density matrix functional theory.

Science.gov (United States)

Lathiotakis, Nektarios N; Helbig, Nicole; Rubio, Angel; Gidopoulos, Nikitas I

2014-10-28

Recently, we introduced [N. N. Lathiotakis, N. Helbig, A. Rubio, and N. I. Gidopoulos, Phys. Rev. A 90, 032511 (2014)] local reduced density matrix functional theory (local RDMFT), a theoretical scheme capable of incorporating static correlation effects in Kohn-Sham equations. Here, we apply local RDMFT to molecular systems of relatively large size, as a demonstration of its computational efficiency and its accuracy in predicting single-electron properties from the eigenvalue spectrum of the single-particle Hamiltonian with a local effective potential. We present encouraging results on the photoelectron spectrum of molecular systems and the relative stability of C20 isotopes. In addition, we propose a modelling of the fractional occupancies as functions of the orbital energies that further improves the efficiency of the method useful in applications to large systems and solids.

Comparing two iteration algorithms of Broyden electron density mixing through an atomic electronic structure computation

International Nuclear Information System (INIS)

Zhang Man-Hong

2016-01-01

By performing the electronic structure computation of a Si atom, we compare two iteration algorithms of Broyden electron density mixing in the literature. One was proposed by Johnson and implemented in the well-known VASP code. The other was given by Eyert. We solve the Kohn-Sham equation by using a conventional outward/inward integration of the differential equation and then connect two parts of solutions at the classical turning points, which is different from the method of the matrix eigenvalue solution as used in the VASP code. Compared to Johnson’s algorithm, the one proposed by Eyert needs fewer total iteration numbers. (paper)

Reduced density-matrix functional theory: Correlation and spectroscopy.

Science.gov (United States)

Di Sabatino, S; Berger, J A; Reining, L; Romaniello, P

2015-07-14

In this work, we explore the performance of approximations to electron correlation in reduced density-matrix functional theory (RDMFT) and of approximations to the observables calculated within this theory. Our analysis focuses on the calculation of total energies, occupation numbers, removal/addition energies, and spectral functions. We use the exactly solvable Hubbard dimer at 1/4 and 1/2 fillings as test systems. This allows us to analyze the underlying physics and to elucidate the origin of the observed trends. For comparison, we also report the results of the GW approximation, where the self-energy functional is approximated, but no further hypothesis is made concerning the approximations of the observables. In particular, we focus on the atomic limit, where the two sites of the dimer are pulled apart and electrons localize on either site with equal probability, unless a small perturbation is present: this is the regime of strong electron correlation. In this limit, using the Hubbard dimer at 1/2 filling with or without a spin-symmetry-broken ground state allows us to explore how degeneracies and spin-symmetry breaking are treated in RDMFT. We find that, within the used approximations, neither in RDMFT nor in GW, the signature of strong correlation is present, when looking at the removal/addition energies and spectral function from the spin-singlet ground state, whereas both give the exact result for the spin-symmetry broken case. Moreover, we show how the spectroscopic properties change from one spin structure to the other.

Optical excitation and electron relaxation dynamics at semiconductor surfaces: a combined approach of density functional and density matrix theory applied to the silicon (001) surface

Energy Technology Data Exchange (ETDEWEB)

Buecking, N

2007-11-05

In this work a new theoretical formalism is introduced in order to simulate the phononinduced relaxation of a non-equilibrium distribution to equilibrium at a semiconductor surface numerically. The non-equilibrium distribution is effected by an optical excitation. The approach in this thesis is to link two conventional, but approved methods to a new, more global description: while semiconductor surfaces can be investigated accurately by density-functional theory, the dynamical processes in semiconductor heterostructures are successfully described by density matrix theory. In this work, the parameters for density-matrix theory are determined from the results of density-functional calculations. This work is organized in two parts. In Part I, the general fundamentals of the theory are elaborated, covering the fundamentals of canonical quantizations as well as the theory of density-functional and density-matrix theory in 2{sup nd} order Born approximation. While the formalism of density functional theory for structure investigation has been established for a long time and many different codes exist, the requirements for density matrix formalism concerning the geometry and the number of implemented bands exceed the usual possibilities of the existing code in this field. A special attention is therefore attributed to the development of extensions to existing formulations of this theory, where geometrical and fundamental symmetries of the structure and the equations are used. In Part II, the newly developed formalism is applied to a silicon (001)surface in a 2 x 1 reconstruction. As first step, density-functional calculations using the LDA functional are completed, from which the Kohn-Sham-wave functions and eigenvalues are used to calculate interaction matrix elements for the electron-phonon-coupling an the optical excitation. These matrix elements are determined for the optical transitions from valence to conduction bands and for electron-phonon processes inside the

Density-matrix formalism for the photoion-electron entanglement in atomic photoionization

International Nuclear Information System (INIS)

Radtke, T.; Fritzsche, S.; Surzhykov, A.

2006-01-01

The density-matrix theory, based on Dirac's relativistic equation, is applied for studying the entanglement between the photoelectron and residual ion in the course of the photoionization of atoms and ions. In particular, emphasis is placed on deriving the final-state density matrix of the overall system 'photoion+electron', including interelectronic effects and the higher multipoles of the radiation field. This final-state density matrix enables one immediately to analyze the change of entanglement as a function of the energy, angle and the polarization of the incoming light. Detailed computations have been carried out for the 5s photoionization of neutral strontium, leading to a photoion in a 5s 2 S J f =1/2 level. It is found that the photoion-electron entanglement decreases significantly near the ionization threshold and that, in general, it depends on both the photon energy and angle. The possibility to extract photoion-electron pairs with a well-defined degree of entanglement may have far-reaching consequences for quantum information and elsewhere

Reduction of multipartite qubit density matrixes to bipartite qubit density matrixes and criteria of partial separability of multipartite qubit density matrixes

OpenAIRE

Zhong, Zai-Zhe

2004-01-01

The partial separability of multipartite qubit density matrixes is strictly defined. We give a reduction way from N-partite qubit density matrixes to bipartite qubit density matrixes, and prove a necessary condition that a N-partite qubit density matrix to be partially separable is its reduced density matrix to satisfy PPT condition.

The density matrix renormalization group method. Application to the EPP model of a cyclic polyene chain

International Nuclear Information System (INIS)

Fano, G.; Ortolani, F.; Ziosi, L.

1997-10-01

The density matrix renormalization group (DMRG) method introduced by White for the study of strongly interacting electron systems is reviewed; the method is variational and considers a system of localized electrons as the union of two adjacent fragments A,B. A density matrix ρ is introduced, whose eigenvectors corresponding to the largest eigenvalues are the most significant, the most probable states of A in the presence of B; these states are retained, while states corresponding to small eigenvalues of ρ are neglected. It is conjectured that the decreasing behaviour of the eigenvalues is gaussian. The DMRG method is tested on the Pariser-Parr-Pople Hamiltonian of a cyclic polyene (CH) N up to N = 34. A Hilbert space of dimension 5. x 10 18 is explored. The ground state energy is 10 -3 eV within the full Cl value in the case N = 18. The DMRG method compares favourably also with coupled cluster approximations. The unrestricted Hartree-Fock solution (which presents spin density waves) is briefly reviewed, and a comparison is made with the DMRG energy values. Finally, the spin-spin and density-density correlation functions are computed; the results suggest that the antiferromagnetic order of the exact solution does not extend up to large distances but exists locally. No charge density waves are present. (author)

Time-dependent occupation numbers in reduced-density-matrix-functional theory: Application to an interacting Landau-Zener model

International Nuclear Information System (INIS)

Requist, Ryan; Pankratov, Oleg

2011-01-01

We prove that if the two-body terms in the equation of motion for the one-body reduced density matrix are approximated by ground-state functionals, the eigenvalues of the one-body reduced density matrix (occupation numbers) remain constant in time. This deficiency is related to the inability of such an approximation to account for relative phases in the two-body reduced density matrix. We derive an exact differential equation giving the functional dependence of these phases in an interacting Landau-Zener model and study their behavior in short- and long-time regimes. The phases undergo resonances whenever the occupation numbers approach the boundaries of the interval [0,1]. In the long-time regime, the occupation numbers display correlation-induced oscillations and the memory dependence of the functionals assumes a simple form.

Hartree--Fock density matrix equation

International Nuclear Information System (INIS)

Cohen, L.; Frishberg, C.

1976-01-01

An equation for the Hartree--Fock density matrix is discussed and the possibility of solving this equation directly for the density matrix instead of solving the Hartree--Fock equation for orbitals is considered. Toward that end the density matrix is expanded in a finite basis to obtain the matrix representative equation. The closed shell case is considered. Two numerical schemes are developed and applied to a number of examples. One example is given where the standard orbital method does not converge while the method presented here does

Reduced one-body density matrix of Tonks–Girardeau gas at finite temperature

International Nuclear Information System (INIS)

Fu Xiao-Chen; Hao Ya-Jiang

2015-01-01

With thermal Bose–Fermi mapping method, we investigate the Tonks–Girardeau gas at finite temperature. It is shown that at low temperature, the Tonks gas displays the Fermi-like density profiles, and with the increase in temperature, the Tonks gas distributes in wider region. The reduced one-body density matrix is diagonal dominant in the whole temperature region, and the off-diagonal elements shall vanish rapidly with the deviation from the diagonal part at high temperature. (paper)

Communication: satisfying fermionic statistics in the modeling of open time-dependent quantum systems with one-electron reduced density matrices.

Science.gov (United States)

Head-Marsden, Kade; Mazziotti, David A

2015-02-07

For an open, time-dependent quantum system, Lindblad derived the most general modification of the quantum Liouville equation in the Markovian approximation that models environmental effects while preserving the non-negativity of the system's density matrix. While Lindblad's modification is correct for N-electron density matrices, solution of the Liouville equation with a Lindblad operator causes the one-electron reduced density matrix (1-RDM) to violate the Pauli exclusion principle. Consequently, after a short time, the 1-RDM is not representable by an ensemble N-electron density matrix (not ensemble N-representable). In this communication, we derive the necessary and sufficient constraints on the Lindbladian matrix within the Lindblad operator to ensure that the 1-RDM remains N-representable for all time. The theory is illustrated by considering the relaxation of an excitation in several molecules F2, N2, CO, and BeH2 subject to environmental noise.

Rank restriction for the variational calculation of two-electron reduced density matrices of many-electron atoms and molecules

International Nuclear Information System (INIS)

Naftchi-Ardebili, Kasra; Hau, Nathania W.; Mazziotti, David A.

2011-01-01

Variational minimization of the ground-state energy as a function of the two-electron reduced density matrix (2-RDM), constrained by necessary N-representability conditions, provides a polynomial-scaling approach to studying strongly correlated molecules without computing the many-electron wave function. Here we introduce a route to enhancing necessary conditions for N representability through rank restriction of the 2-RDM. Rather than adding computationally more expensive N-representability conditions, we directly enhance the accuracy of two-particle (2-positivity) conditions through rank restriction, which removes degrees of freedom in the 2-RDM that are not sufficiently constrained. We select the rank of the particle-hole 2-RDM by deriving the ranks associated with model wave functions, including both mean-field and antisymmetrized geminal power (AGP) wave functions. Because the 2-positivity conditions are exact for quantum systems with AGP ground states, the rank of the particle-hole 2-RDM from the AGP ansatz provides a minimum for its value in variational 2-RDM calculations of general quantum systems. To implement the rank-restricted conditions, we extend a first-order algorithm for large-scale semidefinite programming. The rank-restricted conditions significantly improve the accuracy of the energies; for example, the percentages of correlation energies recovered for HF, CO, and N 2 improve from 115.2%, 121.7%, and 121.5% without rank restriction to 97.8%, 101.1%, and 100.0% with rank restriction. Similar results are found at both equilibrium and nonequilibrium geometries. While more accurate, the rank-restricted N-representability conditions are less expensive computationally than the full-rank conditions.

Kohn-Sham potentials from electron densities using a matrix representation within finite atomic orbital basis sets

Science.gov (United States)

Zhang, Xing; Carter, Emily A.

2018-01-01

We revisit the static response function-based Kohn-Sham (KS) inversion procedure for determining the KS effective potential that corresponds to a given target electron density within finite atomic orbital basis sets. Instead of expanding the potential in an auxiliary basis set, we directly update the potential in its matrix representation. Through numerical examples, we show that the reconstructed density rapidly converges to the target density. Preliminary results are presented to illustrate the possibility of obtaining a local potential in real space from the optimized potential in its matrix representation. We have further applied this matrix-based KS inversion approach to density functional embedding theory. A proof-of-concept study of a solvated proton transfer reaction demonstrates the method's promise.

Quantum Crystallography: Density Matrix-Density Functional Theory and the X-Ray Diffraction Experiment

Science.gov (United States)

Soirat, Arnaud J. A.

Density Matrix Theory is a Quantum Mechanical formalism in which the wavefunction is eliminated and its role taken over by reduced density matrices. The interest of this is that, it allows one, in principle, to calculate any electronic property of a physical system, without having to solve the Schrodinger equation, using only two entities much simpler than an N-body wavefunction: first and second -order reduced density matrices. In practice, though, this very promising possibility faces the tremendous theoretical problem of N-representability, which has been solved for the former, but, until now, voids any hope of theoretically determining the latter. However, it has been shown that single determinant reduced density matrices of any order may be recovered from coherent X-ray diffraction data, if one provides a proper Quantum Mechanical description of the Crystallography experiment. A deeper investigation of this method is the purpose of this work, where we, first, further study the calculation of X-ray reduced density matrices N-representable by a single Slater determinant. In this context, we independently derive necessary and sufficient conditions for the uniqueness of the method. We then show how to account for electron correlation in this model. For the first time, indeed, we derive highly accurate, yet practical, density matrices approximately N-representable by correlated-determinant wavefunctions. The interest of such a result lies in the Quantum Mechanical validity of these density matrices, their property of being entirely obtainable from X-ray coherent diffraction data, their very high accuracy conferred by this known property of the N-representing wavefunction, as well as their definition as explicit functionals of the density. All of these properties are finally used in both a theoretical and a numerical application: in the former, we show that these density matrices may be used in the context of Density Functional Theory to highly accurately determine

Symmetry properties of the electron density and following from it limits on the KS-DFT applications

Science.gov (United States)

Kaplan, Ilya G.

2018-03-01

At present, the Density Functional Theory (DFT) approach elaborated by Kohn with co-authors more than 50 years ago became the most widely used method for study molecules and solids. Using modern computation facilities, it can be applied to systems with million atoms. In the atmosphere of such great popularity, it is particularly important to know the limits of the applicability of DFT methods. In this report, I will discuss two cases when the conventional DFT approaches, using only electron density ρ and its gradients, cannot be applied (I will not consider the Ψ-versions of DFT). The first case is quite evident. In the degenerated states, the electron density may not be defined, since electronic and nuclear motions cannot be separated, the vibronic interaction mixed them. The second case is related to the spin of the state. As it was rigorously proved by group theoretical methods at the theorem level, the electron density does not depend on the total spin S of the arbitrary N-electron state. It means that the Kohn-Sham equations have the same form for states with different S. The critical survey of elaborated DFT procedures, taking into account spin, shows that they modified only exchange functionals, the correlation functionals do not correspond to the spin of the state. The point is that the conception of spin cannot be defined in the framework of the electron density formalism, which corresponds to the one-particle reduced density matrix. This is the main reason of the problems arising in the study by DFT of magnetic properties of the transition metals. The possible way of resolving these problems can be found in the two-particle reduced density matrix formulation of DFT.

Partial differential equation for the idempotent Dirac density matrix characterized solely by the exact non-relativistic ground-state electron density for spherical atomic ions

International Nuclear Information System (INIS)

March, N.H.

2009-08-01

In this Journal, March and Suhai have earlier set up a first-order Dirac idempotent density matrix theory for one- and two-level occupancy in which the only input required is the nonrelativistic ground-state electron density. Here, an analytic generalization is provided for the case of spherical electron densities for arbitrary level occupancy. Be-like atomic ions are referred to as an example, but 'almost spherical' molecules like SiH 4 and GeH 4 also become accessible. (author)

Electronic method for autofluorography of macromolecules on two-d matrices

International Nuclear Information System (INIS)

Case, A.L.; Davidson, J.B.

1983-01-01

A method for detecting, localizing, and quantifying macromolecules contained in a two-dimensional matrix is provided which employs a television-based position sensitive detection system. A molecule-containing matrix may be produced by conventional means to produce spots of light at the molecule locations which are detected by the television system. The matrix , such as a gel matrix, is exposed to an electronic camera system including an image-intensifier and secondary electron conduction camera capable of light integrating times of many minutes. A light image stored in the form of a charge image on the camera tube target is scanned by conventional television techniques, digitized, and stored in a digital memory. Intensity of any point on the image may be determined from the number at the memory address of the point. The entire image may be displayed on a television monitor for inspection and photographing or individual spots may be analyzed through selected readout of the memory locations. Compared to conventional film exposure methods, the exposure time may be reduced 100-1000 times

The problem of the universal density functional and the density matrix functional theory

International Nuclear Information System (INIS)

Bobrov, V. B.; Trigger, S. A.

2013-01-01

The analysis in this paper shows that the Hohenberg-Kohn theorem is the constellation of two statements: (i) the mathematically rigorous Hohenberg-Kohn lemma, which demonstrates that the same ground-state density cannot correspond to two different potentials of an external field, and (ii) the hypothesis of the existence of the universal density functional. Based on the obtained explicit expression for the nonrel-ativistic particle energy in a local external field, we prove that the energy of the system of more than two non-interacting electrons cannot be a functional of the inhomogeneous density. This result is generalized to the system of interacting electrons. It means that the Hohenberg-Kohn lemma cannot provide justification of the universal density functional for fermions. At the same time, statements of the density functional theory remain valid when considering any number of noninteracting ground-state bosons due to the Bose condensation effect. In the framework of the density matrix functional theory, the hypothesis of the existence of the universal density matrix functional corresponds to the cases of noninteracting particles and to interaction in the Hartree-Fock approximation.

Charge-constrained auxiliary-density-matrix methods for the Hartree–Fock exchange contribution

DEFF Research Database (Denmark)

Merlot, Patrick; Izsak, Robert; Borgoo, Alex

2014-01-01

Three new variants of the auxiliary-density-matrix method (ADMM) of Guidon, Hutter, and VandeVondele [J. Chem. Theory Comput. 6, 2348 (2010)] are presented with the common feature thatthey have a simplified constraint compared with the full orthonormality requirement of the earlier ADMM1 method. ....... All ADMM variants are tested for accuracy and performance in all-electron B3LYP calculations with several commonly used basis sets. The effect of the choice of the exchange functional for the ADMM exchange–correction term is also investigated....

The density-matrix renormalization group: a short introduction.

Science.gov (United States)

Schollwöck, Ulrich

2011-07-13

The density-matrix renormalization group (DMRG) method has established itself over the last decade as the leading method for the simulation of the statics and dynamics of one-dimensional strongly correlated quantum lattice systems. The DMRG is a method that shares features of a renormalization group procedure (which here generates a flow in the space of reduced density operators) and of a variational method that operates on a highly interesting class of quantum states, so-called matrix product states (MPSs). The DMRG method is presented here entirely in the MPS language. While the DMRG generally fails in larger two-dimensional systems, the MPS picture suggests a straightforward generalization to higher dimensions in the framework of tensor network states. The resulting algorithms, however, suffer from difficulties absent in one dimension, apart from a much more unfavourable efficiency, such that their ultimate success remains far from clear at the moment.

The time-dependent density matrix renormalisation group method

Science.gov (United States)

Ma, Haibo; Luo, Zhen; Yao, Yao

2018-04-01

Substantial progress of the time-dependent density matrix renormalisation group (t-DMRG) method in the recent 15 years is reviewed in this paper. By integrating the time evolution with the sweep procedures in density matrix renormalisation group (DMRG), t-DMRG provides an efficient tool for real-time simulations of the quantum dynamics for one-dimensional (1D) or quasi-1D strongly correlated systems with a large number of degrees of freedom. In the illustrative applications, the t-DMRG approach is applied to investigate the nonadiabatic processes in realistic chemical systems, including exciton dissociation and triplet fission in polymers and molecular aggregates as well as internal conversion in pyrazine molecule.

Effect of morphology and defect density on electron transfer of electrochemically reduced graphene oxide

Energy Technology Data Exchange (ETDEWEB)

Zhang, Yan, E-mail: yanzhang@sues.edu.cn [School of Material Engineering, Shanghai University of Engineering Science, Shanghai 201620 (China); Hao, Huilian, E-mail: huilian.hao@sues.edu.cn [School of Material Engineering, Shanghai University of Engineering Science, Shanghai 201620 (China); Wang, Linlin, E-mail: wlinlin@mail.ustc.edu.cn [College of Chemistry and Chemical Engineering, Shanghai University of Engineering Science, Shanghai 201620 (China)

2016-12-30

Highlights: • Different morphologies of ERGO on the surface of GCE were prepared via different methods. • The defect densities of ERGO were controlled by tuning the mass or concentration of GO. • A higher defect density of ERGO accelerates electron transfer rate. • ERGO with more exposed edge planes shows significantly higher electron transfer kinetics. • Both edge planes and defect density contribute to electron transfer of ERGO. - Abstract: Electrochemically reduced graphene oxide (ERGO) is widely used to construct electrochemical sensors. Understanding the electron transfer behavior of ERGO is essential for its electrode material applications. In this paper, different morphologies of ERGO were prepared via two different methods. Compared to ERGO/GCEs prepared by electrochemical reduction of pre-deposited GO, more exposed edge planes of ERGO are observed on the surface of ERGO-GCE that was constructed by electrophoretic deposition of GO. The defect densities of ERGO were controlled by tuning the mass or concentration of GO. The electron transfer kinetics (k{sup 0}) of GCE with different ERGOs was comparatively investigated. Owing to increased surface areas and decreased defect density, the k{sup 0} values of ERGO/GCE initially increase and then decrease with incrementing of GO mass. When the morphology and surface real areas of ERGO-GCE are the same, an increased defect density induces an accelerated electron transfer rate. k{sup 0} valuesof ERGO-GCEs are about 1 order of magnitude higher than those of ERGO/GCEs due to the difference in the amount of edge planes. This work demonstrates that both defect densities and edge planes of ERGO play crucial roles in electron transfer kinetics.

Effect of morphology and defect density on electron transfer of electrochemically reduced graphene oxide

International Nuclear Information System (INIS)

Zhang, Yan; Hao, Huilian; Wang, Linlin

2016-01-01

Highlights: • Different morphologies of ERGO on the surface of GCE were prepared via different methods. • The defect densities of ERGO were controlled by tuning the mass or concentration of GO. • A higher defect density of ERGO accelerates electron transfer rate. • ERGO with more exposed edge planes shows significantly higher electron transfer kinetics. • Both edge planes and defect density contribute to electron transfer of ERGO. - Abstract: Electrochemically reduced graphene oxide (ERGO) is widely used to construct electrochemical sensors. Understanding the electron transfer behavior of ERGO is essential for its electrode material applications. In this paper, different morphologies of ERGO were prepared via two different methods. Compared to ERGO/GCEs prepared by electrochemical reduction of pre-deposited GO, more exposed edge planes of ERGO are observed on the surface of ERGO-GCE that was constructed by electrophoretic deposition of GO. The defect densities of ERGO were controlled by tuning the mass or concentration of GO. The electron transfer kinetics (k"0) of GCE with different ERGOs was comparatively investigated. Owing to increased surface areas and decreased defect density, the k"0 values of ERGO/GCE initially increase and then decrease with incrementing of GO mass. When the morphology and surface real areas of ERGO-GCE are the same, an increased defect density induces an accelerated electron transfer rate. k"0 valuesof ERGO-GCEs are about 1 order of magnitude higher than those of ERGO/GCEs due to the difference in the amount of edge planes. This work demonstrates that both defect densities and edge planes of ERGO play crucial roles in electron transfer kinetics.

Dynamical simulation of electron transfer processes in self-assembled monolayers at metal surfaces using a density matrix approach.

Science.gov (United States)

Prucker, V; Bockstedte, M; Thoss, M; Coto, P B

2018-03-28

A single-particle density matrix approach is introduced to simulate the dynamics of heterogeneous electron transfer (ET) processes at interfaces. The characterization of the systems is based on a model Hamiltonian parametrized by electronic structure calculations and a partitioning method. The method is applied to investigate ET in a series of nitrile-substituted (poly)(p-phenylene)thiolate self-assembled monolayers adsorbed at the Au(111) surface. The results show a significant dependence of the ET on the orbital symmetry of the donor state and on the molecular and electronic structure of the spacer.

Dynamical simulation of electron transfer processes in self-assembled monolayers at metal surfaces using a density matrix approach

Science.gov (United States)

Prucker, V.; Bockstedte, M.; Thoss, M.; Coto, P. B.

2018-03-01

A single-particle density matrix approach is introduced to simulate the dynamics of heterogeneous electron transfer (ET) processes at interfaces. The characterization of the systems is based on a model Hamiltonian parametrized by electronic structure calculations and a partitioning method. The method is applied to investigate ET in a series of nitrile-substituted (poly)(p-phenylene)thiolate self-assembled monolayers adsorbed at the Au(111) surface. The results show a significant dependence of the ET on the orbital symmetry of the donor state and on the molecular and electronic structure of the spacer.

Feasibility Studies of the Two Filters Method in TJ-II for Electron Temperature Measurements in High Density Plasmas

International Nuclear Information System (INIS)

Baiao, D.; Medina, F.; Ochando, M.; Varandas, C.

2009-01-01

The TJ-II plasma soft X-ray emission was studied in order to establish an adequate setup for an electron temperature diagnostic suitable for high density, with spatial and temporal resolutions, based on the two-filters method. The preliminary experimental results reported were obtained with two diagnostics (an X-ray PHA based on a Ge detector and a tomography system) already installed in TJ-II stellarator. These results lead to the conclusion that the two-filters method was a suitable option for an electron temperature diagnostic for high-density plasmas in TJ-II. We present the design and fi rst results obtained with a prototype for the measurement of electron temperature in TJ-II plasmas heated with energetic neutral beams. This system consists in two AXUV20A detectors which measure the soft X-ray plasma emissivity trough beryllium filters of different thickness. From the two-filters technique it is possible to estimate the electron temperature. The analyses carried out allowed concluding which filter thicknesses are most suited for TJ-II plasmas, and enhanced the need of a computer code to simulate signals and plasma compositions. (Author) 7 refs.

Feasibility Studies of the Two Filters Method in TJ-II for Electron Temperature Measurements in High Density Plasmas

Energy Technology Data Exchange (ETDEWEB)

Baiao, D.; Medina, F.; Ochando, M.; Varandas, C.

2009-07-01

The TJ-II plasma soft X-ray emission was studied in order to establish an adequate setup for an electron temperature diagnostic suitable for high density, with spatial and temporal resolutions, based on the two-filters method. The preliminary experimental results reported were obtained with two diagnostics (an X-ray PHA based on a Ge detector and a tomography system) already installed in TJ-II stellarator. These results lead to the conclusion that the two-filters method was a suitable option for an electron temperature diagnostic for high-density plasmas in TJ-II. We present the design and fi rst results obtained with a prototype for the measurement of electron temperature in TJ-II plasmas heated with energetic neutral beams. This system consists in two AXUV20A detectors which measure the soft X-ray plasma emissivity trough beryllium filters of different thickness. From the two-filters technique it is possible to estimate the electron temperature. The analyses carried out allowed concluding which filter thicknesses are most suited for TJ-II plasmas, and enhanced the need of a computer code to simulate signals and plasma compositions. (Author) 7 refs.

Imaginary time density-density correlations for two-dimensional electron gases at high density

Energy Technology Data Exchange (ETDEWEB)

Motta, M.; Galli, D. E. [Dipartimento di Fisica, Università degli Studi di Milano, Via Celoria 16, 20133 Milano (Italy); Moroni, S. [IOM-CNR DEMOCRITOS National Simulation Center and SISSA, Via Bonomea 265, 34136 Trieste (Italy); Vitali, E. [Department of Physics, College of William and Mary, Williamsburg, Virginia 23187-8795 (United States)

2015-10-28

We evaluate imaginary time density-density correlation functions for two-dimensional homogeneous electron gases of up to 42 particles in the continuum using the phaseless auxiliary field quantum Monte Carlo method. We use periodic boundary conditions and up to 300 plane waves as basis set elements. We show that such methodology, once equipped with suitable numerical stabilization techniques necessary to deal with exponentials, products, and inversions of large matrices, gives access to the calculation of imaginary time correlation functions for medium-sized systems. We discuss the numerical stabilization techniques and the computational complexity of the methodology and we present the limitations related to the size of the systems on a quantitative basis. We perform the inverse Laplace transform of the obtained density-density correlation functions, assessing the ability of the phaseless auxiliary field quantum Monte Carlo method to evaluate dynamical properties of medium-sized homogeneous fermion systems.

Obtaining Hartree-Fock and density functional theory doubly excited states with Car-Parrinello density matrix search

Science.gov (United States)

Liang, Wenkel; Isborn, Christine M.; Li, Xiaosong

2009-11-01

The calculation of doubly excited states is one of the major problems plaguing the modern day excited state workhorse methodology of linear response time dependent Hartree-Fock (TDHF) and density function theory (TDDFT). We have previously shown that the use of a resonantly tuned field within real-time TDHF and TDDFT is able to simultaneously excite both the α and β electrons to achieve the two-electron excited states of minimal basis H2 and HeH+ [C. M. Isborn and X. Li, J. Chem. Phys. 129, 204107 (2008)]. We now extend this method to many electron systems with the use of our Car-Parrinello density matrix search (CP-DMS) with a first-principles fictitious mass method for wave function optimization [X. Li, C. L. Moss, W. Liang, and Y. Feng, J. Chem. Phys. 130, 234115 (2009)]. Real-time TDHF/TDDFT is used during the application of the laser field perturbation, driving the electron density toward the doubly excited state. The CP-DMS method then converges the density to the nearest stationary state. We present these stationary state doubly excited state energies and properties at the HF and DFT levels for H2, HeH+, lithium hydride, ethylene, and butadiene.

Two-dimensional electron density characterisation of arc interruption phenomenon in current-zero phase

Science.gov (United States)

Inada, Yuki; Kamiya, Tomoki; Matsuoka, Shigeyasu; Kumada, Akiko; Ikeda, Hisatoshi; Hidaka, Kunihiko

2018-01-01

Two-dimensional electron density imaging over free burning SF6 arcs and SF6 gas-blast arcs was conducted at current zero using highly sensitive Shack-Hartmann type laser wavefront sensors in order to experimentally characterise electron density distributions for the success and failure of arc interruption in the thermal reignition phase. The experimental results under an interruption probability of 50% showed that free burning SF6 arcs with axially asymmetric electron density profiles were interrupted with a success rate of 88%. On the other hand, the current interruption of SF6 gas-blast arcs was reproducibly achieved under locally reduced electron densities and the interruption success rate was 100%.

Nonorthogonal orbital based N-body reduced density matrices and their applications to valence bond theory. I. Hamiltonian matrix elements between internally contracted excited valence bond wave functions

Science.gov (United States)

Chen, Zhenhua; Chen, Xun; Wu, Wei

2013-04-01

In this series, the n-body reduced density matrix (n-RDM) approach for nonorthogonal orbitals and their applications to ab initio valence bond (VB) methods are presented. As the first paper of this series, Hamiltonian matrix elements between internally contracted VB wave functions are explicitly provided by means of nonorthogonal orbital based RDM approach. To this end, a more generalized Wick's theorem, called enhanced Wick's theorem, is presented both in arithmetical and in graphical forms, by which the deduction of expressions for the matrix elements between internally contracted VB wave functions is dramatically simplified, and the matrix elements are finally expressed in terms of tensor contractions of electronic integrals and n-RDMs of the reference VB self-consistent field wave function. A string-based algorithm is developed for the purpose of evaluating n-RDMs in an efficient way. Using the techniques presented in this paper, one is able to develop new methods and efficient algorithms for nonorthogonal orbital based many-electron theory much easier than by use of the first quantized formulism.

Simulation of electron beam from two strip electron guns and control of power density by rotation of gun

International Nuclear Information System (INIS)

Sahu, G K; Baruah, S; Thakur, K B

2012-01-01

Electron beam is preferably used for large scale evaporation of refractory materials. Material evaporation from a long and narrow source providing a well collimated wedge shaped atomic beam has applications in isotopic purification of metals relevant to nuclear industry. The electron beam from an electron gun with strip type filament provides a linear heating source. However, the high power density of the electron beam can lead to turbulence of the melt pool and undesirable splashing of molten metal. For obtaining quiet surface evaporation, the linear electron beam is generally scanned along its length. To further reduce the power density to maintain quiet evaporation the width of the vapour source can be controlled by rotating the electron gun on its plane, thereby scanning an inclined beam over the molten pool. The rotation of gun has further advantages. When multiple strip type electron guns are used for scaling up evaporation length, a dark zone appears between two beams due to physical separation of adjacent guns. This dark zone can be reduced by rotating the gun and thereby bringing two adjacent beams closer. The paper presented here provides the simulation results of the electron beam trajectory and incident power density originating from two strip electron guns by using in-house developed code. The effect of electron gun rotation on the electron beam trajectory and power density is studied. The simulation result is experimentally verified with the image of molten pool and heat affected zone taken after experiment. This technique can be gainfully utilized in controlling the time averaged power density of the electron beam and obtaining quiet evaporation from the metal molten pool.

Orbital functionals in density-matrix- and current-density-functional theory

Energy Technology Data Exchange (ETDEWEB)

Helbig, N

2006-05-15

Density-Functional Theory (DFT), although widely used and very successful in the calculation of several observables, fails to correctly describe strongly correlated materials. In the first part of this work we, therefore, introduce reduced-densitymatrix- functional theory (RDMFT) which is one possible way to treat electron correlation beyond DFT. Within this theory the one-body reduced density matrix (1- RDM) is used as the basic variable. Our main interest is the calculation of the fundamental gap which proves very problematic within DFT. In order to calculate the fundamental gap we generalize RDMFT to fractional particle numbers M by describing the system as an ensemble of an N and an N+1 particle system (with N{<=}M{<=}N+1). For each fixed particle number, M, the total energy is minimized with respect to the natural orbitals and their occupation numbers. This leads to the total energy as a function of M. The derivative of this function with respect to the particle number has a discontinuity at integer particle number which is identical to the gap. In addition, we investigate the necessary and sufficient conditions for the 1- RDM of a system with fractional particle number to be N-representable. Numerical results are presented for alkali atoms, small molecules, and periodic systems. Another problem within DFT is the description of non-relativistic many-electron systems in the presence of magnetic fields. It requires the paramagnetic current density and the spin magnetization to be used as basic variables besides the electron density. However, electron-gas-based functionals of current-spin-density-functional Theory (CSDFT) exhibit derivative discontinuities as a function of the magnetic field whenever a new Landau level is occupied, which makes them difficult to use in practice. Since the appearance of Landau levels is, intrinsically, an orbital effect it is appealing to use orbital-dependent functionals. We have developed a CSDFT version of the optimized

Local energy equation for two-electron atoms and relation between kinetic energy and electron densities

International Nuclear Information System (INIS)

March, N.H.

2002-08-01

In early work, Dawson and March [J. Chem. Phys. 81, 5850 (1984)] proposed a local energy method for treating both Hartree-Fock and correlated electron theory. Here, an exactly solvable model two-electron atom with pure harmonic interactions is treated in its ground state in the above context. A functional relation between the kinetic energy density t(r) at the origin r=0 and the electron density p(r) at the same point then emerges. The same approach is applied to the Hookean atom; in which the two electrons repel with Coulombic energy e 2 /r 12 , with r 12 the interelectronic separation, but are still harmonically confined. Again the kinetic energy density t(r) is the focal point, but now generalization away from r=0 is also effected. Finally, brief comments are added about He-like atomic ions in the limit of large atomic number. (author)

Theory of open quantum systems with bath of electrons and phonons and spins: many-dissipaton density matrixes approach.

Science.gov (United States)

Yan, YiJing

2014-02-07

This work establishes a strongly correlated system-and-bath dynamics theory, the many-dissipaton density operators formalism. It puts forward a quasi-particle picture for environmental influences. This picture unifies the physical descriptions and algebraic treatments on three distinct classes of quantum environments, electron bath, phonon bath, and two-level spin or exciton bath, as their participating in quantum dissipation processes. Dynamical variables for theoretical description are no longer just the reduced density matrix for system, but remarkably also those for quasi-particles of bath. The present theoretical formalism offers efficient and accurate means for the study of steady-state (nonequilibrium and equilibrium) and real-time dynamical properties of both systems and hybridizing environments. It further provides universal evaluations, exact in principle, on various correlation functions, including even those of environmental degrees of freedom in coupling with systems. Induced environmental dynamics could be reflected directly in experimentally measurable quantities, such as Fano resonances and quantum transport current shot noise statistics.

Electron paramagnetic resonance g-tensors from state interaction spin-orbit coupling density matrix renormalization group

Science.gov (United States)

Sayfutyarova, Elvira R.; Chan, Garnet Kin-Lic

2018-05-01

We present a state interaction spin-orbit coupling method to calculate electron paramagnetic resonance g-tensors from density matrix renormalization group wavefunctions. We apply the technique to compute g-tensors for the TiF3 and CuCl42 - complexes, a [2Fe-2S] model of the active center of ferredoxins, and a Mn4CaO5 model of the S2 state of the oxygen evolving complex. These calculations raise the prospects of determining g-tensors in multireference calculations with a large number of open shells.

Full Quantum Dynamics Simulation of a Realistic Molecular System Using the Adaptive Time-Dependent Density Matrix Renormalization Group Method.

Science.gov (United States)

Yao, Yao; Sun, Ke-Wei; Luo, Zhen; Ma, Haibo

2018-01-18

The accurate theoretical interpretation of ultrafast time-resolved spectroscopy experiments relies on full quantum dynamics simulations for the investigated system, which is nevertheless computationally prohibitive for realistic molecular systems with a large number of electronic and/or vibrational degrees of freedom. In this work, we propose a unitary transformation approach for realistic vibronic Hamiltonians, which can be coped with using the adaptive time-dependent density matrix renormalization group (t-DMRG) method to efficiently evolve the nonadiabatic dynamics of a large molecular system. We demonstrate the accuracy and efficiency of this approach with an example of simulating the exciton dissociation process within an oligothiophene/fullerene heterojunction, indicating that t-DMRG can be a promising method for full quantum dynamics simulation in large chemical systems. Moreover, it is also shown that the proper vibronic features in the ultrafast electronic process can be obtained by simulating the two-dimensional (2D) electronic spectrum by virtue of the high computational efficiency of the t-DMRG method.

Time-dependent density-functional tight-binding method with the third-order expansion of electron density.

Science.gov (United States)

Nishimoto, Yoshio

2015-09-07

We develop a formalism for the calculation of excitation energies and excited state gradients for the self-consistent-charge density-functional tight-binding method with the third-order contributions of a Taylor series of the density functional theory energy with respect to the fluctuation of electron density (time-dependent density-functional tight-binding (TD-DFTB3)). The formulation of the excitation energy is based on the existing time-dependent density functional theory and the older TD-DFTB2 formulae. The analytical gradient is computed by solving Z-vector equations, and it requires one to calculate the third-order derivative of the total energy with respect to density matrix elements due to the inclusion of the third-order contributions. The comparison of adiabatic excitation energies for selected small and medium-size molecules using the TD-DFTB2 and TD-DFTB3 methods shows that the inclusion of the third-order contributions does not affect excitation energies significantly. A different set of parameters, which are optimized for DFTB3, slightly improves the prediction of adiabatic excitation energies statistically. The application of TD-DFTB for the prediction of absorption and fluorescence energies of cresyl violet demonstrates that TD-DFTB3 reproduced the experimental fluorescence energy quite well.

The finite temperature density matrix and two-point correlations in the antiferromagnetic XXZ chain

Science.gov (United States)

Göhmann, Frank; Hasenclever, Nils P.; Seel, Alexander

2005-10-01

We derive finite temperature versions of integral formulae for the two-point correlation functions in the antiferromagnetic XXZ chain. The derivation is based on the summation of density matrix elements characterizing a finite chain segment of length m. On this occasion we also supply a proof of the basic integral formula for the density matrix presented in an earlier publication.

Truncation scheme of time-dependent density-matrix approach II

Energy Technology Data Exchange (ETDEWEB)

Tohyama, Mitsuru [Kyorin University School of Medicine, Mitaka, Tokyo (Japan); Schuck, Peter [Institut de Physique Nucleaire, IN2P3-CNRS, Universite Paris-Sud, Orsay (France); Laboratoire de Physique et de Modelisation des Milieux Condenses, CNRS et Universite Joseph Fourier, Grenoble (France)

2017-09-15

A truncation scheme of the Bogoliubov-Born-Green-Kirkwood-Yvon hierarchy for reduced density matrices, where a three-body density matrix is approximated by two-body density matrices, is improved to take into account a normalization effect. The truncation scheme is tested for the Lipkin model. It is shown that the obtained results are in good agreement with the exact solutions. (orig.)

Single-particle density matrix and superfluidity in the two-dimensional Bose Coulomb fluid

International Nuclear Information System (INIS)

Minguzzi, A.; Tosi, M.P.; Davoudi, B.

2002-01-01

A study by Magro and Ceperley [Phys. Rev. Lett. 73, 826 (1994)] has shown that the ground state of the two-dimensional fluid of charged bosons with logarithmic interactions is not Bose condensed, but exhibits algebraic off-diagonal order in the single-particle density matrix ρ(r). We use a hydrodynamic Hamiltonian expressed in terms of density and phase operators, in combination with an f-sum rule on the superfluid fraction, to reproduce these results and to extend the evaluation of the density matrix to finite temperature T. This approach allows us to treat the liquid as a superfluid in the absence of a condensate. The algebraic decay of the one-body density matrix is due to correlations between phase fluctuations, and we find that the exponent in the power law is determined by the superfluid density n s (T). We also find that the plasmon gap in the single-particle energy spectrum at long wavelengths decreases with increasing T and closes at the critical temperature for the onset of superfluidity

Multichannel quantum defect and reduced R-matrix

International Nuclear Information System (INIS)

Hategan, C.; Ionescu, R.A.; Cutoiu, D.; Gugiu, M.

2002-01-01

The collision of an electron with the atomic electronic core or the scattering of a nucleon on the atomic nucleus, usually, result into multiparticle excitations producing a resonance of a compound system, followed by its decay in reaction channels. Both in the electron-atom collisions and in nucleon-nucleus reactions, these multichannel resonances are described by poles of all R-Matrix elements. The resonances originating in single particle states, either in electron-atom collision or in nucleon-nucleus scattering, are approached in quite different descriptions. For example, the single-particle resonance in nuclear scattering is described, in R-Matrix Theory, by a perturbative method due to Bloch. The original single-nucleon state overlaps the actual states of the nucleus, resulting into a micro-giant description of the single particle resonance. The spectroscopic aspects of the single particle state, mixed with actual nuclear states, are subject of nucleon (or single particle) Strength Function. The electron, involving single particle Rydberg state in an atomic collision, 'avoids' its wave function mixing with that of inner multielectron core, because it is spatially far-away located from that core. This process is usually described by the Multichannel Quantum Defect Theory (MQDT). In the electron-atom scattering rather the effect of inner multielectron core on Rydberg electrons is studied by means of a global parameter, historically called 'Quantum Defect'. Both these types of resonances have in common the preserving of the single-particle wave function in a complex system with multiparticle excitations. In this work one approaches description of single-particle (electron or nucleon) resonance in a multichannel system. The single particle multichannel resonances are not longer described by a R-Matrix pole (specific for resonances originating in multiparticle excitations) but rather by a natural method for incorporating a single particle state in R-Matrix Theory

Unraveling multi-spin effects in rotational resonance nuclear magnetic resonance using effective reduced density matrix theory

International Nuclear Information System (INIS)

SivaRanjan, Uppala; Ramachandran, Ramesh

2014-01-01

A quantum-mechanical model integrating the concepts of reduced density matrix and effective Hamiltonians is proposed to explain the multi-spin effects observed in rotational resonance (R 2 ) nuclear magnetic resonance (NMR) experiments. Employing this approach, the spin system of interest is described in a reduced subspace inclusive of its coupling to the surroundings. Through suitable model systems, the utility of our theory is demonstrated and verified with simulations emerging from both analytic and numerical methods. The analytic results presented in this article provide an accurate description/interpretation of R 2 experimental results and could serve as a test-bed for distinguishing coherent/incoherent effects in solid-state NMR

Unraveling multi-spin effects in rotational resonance nuclear magnetic resonance using effective reduced density matrix theory

Energy Technology Data Exchange (ETDEWEB)

SivaRanjan, Uppala; Ramachandran, Ramesh, E-mail: rramesh@iisermohali.ac.in [Department of Chemical Sciences, Indian Institute of Science Education and Research (IISER) Mohali, Sector 81, Manauli, P.O. Box-140306, Mohali, Punjab (India)

2014-02-07

A quantum-mechanical model integrating the concepts of reduced density matrix and effective Hamiltonians is proposed to explain the multi-spin effects observed in rotational resonance (R{sup 2}) nuclear magnetic resonance (NMR) experiments. Employing this approach, the spin system of interest is described in a reduced subspace inclusive of its coupling to the surroundings. Through suitable model systems, the utility of our theory is demonstrated and verified with simulations emerging from both analytic and numerical methods. The analytic results presented in this article provide an accurate description/interpretation of R{sup 2} experimental results and could serve as a test-bed for distinguishing coherent/incoherent effects in solid-state NMR.

Design of Electronic Medical Record User Interfaces: A Matrix-Based Method for Improving Usability

Directory of Open Access Journals (Sweden)

Kushtrim Kuqi

2013-01-01

Full Text Available This study examines a new approach of using the Design Structure Matrix (DSM modeling technique to improve the design of Electronic Medical Record (EMR user interfaces. The usability of an EMR medication dosage calculator used for placing orders in an academic hospital setting was investigated. The proposed method captures and analyzes the interactions between user interface elements of the EMR system and groups elements based on information exchange, spatial adjacency, and similarity to improve screen density and time-on-task. Medication dose adjustment task time was recorded for the existing and new designs using a cognitive simulation model that predicts user performance. We estimate that the design improvement could reduce time-on-task by saving an average of 21 hours of hospital physicians’ time over the course of a month. The study suggests that the application of DSM can improve the usability of an EMR user interface.

Quantum Entanglement and Reduced Density Matrices

Science.gov (United States)

Purwanto, Agus; Sukamto, Heru; Yuwana, Lila

2018-05-01

We investigate entanglement and separability criteria of multipartite (n-partite) state by examining ranks of its reduced density matrices. Firstly, we construct the general formula to determine the criterion. A rank of origin density matrix always equals one, meanwhile ranks of reduced matrices have various ranks. Next, separability and entanglement criterion of multipartite is determined by calculating ranks of reduced density matrices. In this article we diversify multipartite state criteria into completely entangled state, completely separable state, and compound state, i.e. sub-entangled state and sub-entangledseparable state. Furthermore, we also shorten the calculation proposed by the previous research to determine separability of multipartite state and expand the methods to be able to differ multipartite state based on criteria above.

Implementing and evaluating a fictitious electron dynamics method for the calculation of electronic structure: Application to the Si(100) surface

International Nuclear Information System (INIS)

Hoffman, M J H; Claassens, C H

2006-01-01

A density matrix based fictitious electron dynamics method for calculating electronic structure has been implemented within a semi-empirical quantum chemistry environment. This method uses an equation of motion that implicitly ensures the idempotency constraint on the density matrix. Test calculations showed that this method has potential of being combined with simultaneous atomic dynamics, in analogy to the popular Car-Parrinello method. In addition, the sparsity of the density matrix and the sophisticated though flexible way of ensuring idempotency conservation while integrating the equation of motion creates the potential of developing a fast linear scaling method

The validity of the density scaling method in primary electron transport for photon and electron beams

International Nuclear Information System (INIS)

Woo, M.K.; Cunningham, J.R.

1990-01-01

In the convolution/superposition method of photon beam dose calculations, inhomogeneities are usually handled by using some form of scaling involving the relative electron densities of the inhomogeneities. In this paper the accuracy of density scaling as applied to primary electrons generated in photon interactions is examined. Monte Carlo calculations are compared with density scaling calculations for air and cork slab inhomogeneities. For individual primary photon kernels as well as for photon interactions restricted to a thin layer, the results can differ significantly, by up to 50%, between the two calculations. However, for realistic photon beams where interactions occur throughout the whole irradiated volume, the discrepancies are much less severe. The discrepancies for the kernel calculation are attributed to the scattering characteristics of the electrons and the consequent oversimplified modeling used in the density scaling method. A technique called the kernel integration technique is developed to analyze the general effects of air and cork inhomogeneities. It is shown that the discrepancies become significant only under rather extreme conditions, such as immediately beyond the surface after a large air gap. In electron beams all the primary electrons originate from the surface of the phantom and the errors caused by simple density scaling can be much more significant. Various aspects relating to the accuracy of density scaling for air and cork slab inhomogeneities are discussed

The localization-delocalization matrix and the electron-density-weighted connectivity matrix of a finite graphene nanoribbon reconstructed from kernel fragments.

Science.gov (United States)

Timm, Matthew J; Matta, Chérif F; Massa, Lou; Huang, Lulu

2014-11-26

Bader's quantum theory of atoms in molecules (QTAIM) and chemical graph theory, merged in the localization-delocalization matrices (LDMs) and the electron-density-weighted connectivity matrices (EDWCM), are shown to benefit in computational speed from the kernel energy method (KEM). The LDM and EDWCM quantum chemical graph matrices of a 66-atom C46H20 hydrogen-terminated armchair graphene nanoribbon, in 14 (2×7) rings of C2v symmetry, are accurately reconstructed from kernel fragments. (This includes the full sets of electron densities at 84 bond critical points and 19 ring critical points, and the full sets of 66 localization and 4290 delocalization indices (LIs and DIs).) The average absolute deviations between KEM and directly calculated atomic electron populations, obtained from the sum of the LIs and half of the DIs of an atom, are 0.0012 ± 0.0018 e(-) (∼0.02 ± 0.03%) for carbon atoms and 0.0007 ± 0.0003 e(-) (∼0.01 ± 0.01%) for hydrogen atoms. The integration errors in the total electron population (296 electrons) are +0.0003 e(-) for the direct calculation (+0.0001%) and +0.0022 e(-) for KEM (+0.0007%). The accuracy of the KEM matrix elements is, thus, probably of the order of magnitude of the combined precision of the electronic structure calculation and the atomic integrations. KEM appears capable of delivering not only the total energies with chemical accuracy (which is well documented) but also local and nonlocal properties accurately, including the DIs between the fragments (crossing fragmentation lines). Matrices of the intact ribbon, the kernels, the KEM-reconstructed ribbon, and errors are available as Supporting Information .

Density matrix in quantum electrodynamics, equivalence principle and Hawking effect

International Nuclear Information System (INIS)

Frolov, V.P.; Gitman, D.M.

1978-01-01

The expression for the density matrix describing particles of one sort (electrons or positrons) created by an external electromagnetic field from the vacuum is obtained. The explicit form of the density matrix is found for the case of constant and uniform electric field. Arguments are given for the presence of a connection between the thermal nature of the density matrix describing particles created by the gravitational field of a black hole and the equivalence principle. (author)

Electron-electron scattering in linear transport in two-dimensional systems

DEFF Research Database (Denmark)

Hu, Ben Yu-Kuang; Flensberg, Karsten

1996-01-01

We describe a method for numerically incorporating electron-electron scattering in quantum wells for small deviations of the distribution function from equilibrium, within the framework of the Boltzmann equation. For a given temperature T and density n, a symmetric matrix needs to be evaluated only...... once, and henceforth it can be used to describe electron-electron scattering in any Boltzmann equation linear-response calculation for that particular T and n. Using this method, we calculate the distribution function and mobility for electrons in a quantum well, including full finite...

Quantum entanglement and thermal reduced density matrices in fermion and spin systems on ladders

International Nuclear Information System (INIS)

Chen, Xiao; Fradkin, Eduardo

2013-01-01

Numerical studies of the reduced density matrix of a gapped spin-1/2 Heisenberg antiferromagnet on a two-leg ladder find that it has the same form as the Gibbs density matrix of a gapless spin-1/2 Heisenberg antiferromagnetic chain at a finite temperature determined by the spin gap of the ladder. We investigate this interesting result by considering a model of free fermions on a two-leg ladder (gapped by the inter-chain tunneling operator) and in spin systems on a ladder with a gapped ground state using exact solutions and several controlled approximations. We calculate the reduced density matrix and the entanglement entropy for a leg of the ladder (i.e. a cut made between the chains). In the fermionic system we find the exact form of the reduced density matrix for one of the chains and determine the entanglement spectrum explicitly. Here we find that in the weak tunneling limit of the ladder the entanglement entropy of one chain of the gapped ladder has a simple and universal form dictated by conformal invariance. In the case of the spin system, we consider the strong coupling limit by using perturbation theory and get the reduced density matrix by the Schmidt decomposition. The entanglement entropies of a general gapped system of two coupled conformal field theories (in 1 + 1 dimensions) are discussed using the replica trick and scaling arguments. We show that (1) for a system with a bulk gap the reduced density matrix has the form of a thermal density matrix, (2) the long-wavelength modes of one subsystem (a chain) of a gapped coupled system are always thermal, (3) the von Neumann entropy equals the thermodynamic entropy of one chain, and (4) the bulk gap plays the role of effective temperature. (paper)

Watching excitons move: the time-dependent transition density matrix

Science.gov (United States)

Ullrich, Carsten

2012-02-01

Time-dependent density-functional theory allows one to calculate excitation energies and the associated transition densities in principle exactly. The transition density matrix (TDM) provides additional information on electron-hole localization and coherence of specific excitations of the many-body system. We have extended the TDM concept into the real-time domain in order to visualize the excited-state dynamics in conjugated molecules. The time-dependent TDM is defined as an implicit density functional, and can be approximately obtained from the time-dependent Kohn-Sham orbitals. The quality of this approximation is assessed in simple model systems. A computational scheme for real molecular systems is presented: the time-dependent Kohn-Sham equations are solved with the OCTOPUS code and the time-dependent Kohn-Sham TDM is calculated using a spatial partitioning scheme. The method is applied to show in real time how locally created electron-hole pairs spread out over neighboring conjugated molecular chains. The coupling mechanism, electron-hole coherence, and the possibility of charge separation are discussed.

Fingerprint-based structure retrieval using electron density.

Science.gov (United States)

Yin, Shuangye; Dokholyan, Nikolay V

2011-03-01

We present a computational approach that can quickly search a large protein structural database to identify structures that fit a given electron density, such as determined by cryo-electron microscopy. We use geometric invariants (fingerprints) constructed using 3D Zernike moments to describe the electron density, and reduce the problem of fitting of the structure to the electron density to simple fingerprint comparison. Using this approach, we are able to screen the entire Protein Data Bank and identify structures that fit two experimental electron densities determined by cryo-electron microscopy. Copyright © 2010 Wiley-Liss, Inc.

Method for controlling low-energy high current density electron beams

International Nuclear Information System (INIS)

Lee, J.N.; Oswald, R.B. Jr.

1977-01-01

A method and an apparatus for controlling the angle of incidence of low-energy, high current density electron beams are disclosed. The apparatus includes a current generating diode arrangement with a mesh anode for producing a drifting electron beam. An auxiliary grounded screen electrode is placed between the anode and a target for controlling the average angle of incidence of electrons in the drifting electron beam. According to the method of the present invention, movement of the auxiliary screen electrode relative to the target and the anode permits reliable and reproducible adjustment of the average angle of incidence of the electrons in low energy, high current density relativistic electron beams

Correlation density matrices for one-dimensional quantum chains based on the density matrix renormalization group

International Nuclear Information System (INIS)

Muender, W; Weichselbaum, A; Holzner, A; Delft, Jan von; Henley, C L

2010-01-01

A useful concept for finding numerically the dominant correlations of a given ground state in an interacting quantum lattice system in an unbiased way is the correlation density matrix (CDM). For two disjoint, separated clusters, it is defined to be the density matrix of their union minus the direct product of their individual density matrices and contains all the correlations between the two clusters. We show how to extract from the CDM a survey of the relative strengths of the system's correlations in different symmetry sectors and the nature of their decay with distance (power law or exponential), as well as detailed information on the operators carrying long-range correlations and the spatial dependence of their correlation functions. To achieve this goal, we introduce a new method of analysing the CDM, termed the dominant operator basis (DOB) method, which identifies in an unbiased fashion a small set of operators for each cluster that serve as a basis for the dominant correlations of the system. We illustrate this method by analysing the CDM for a spinless extended Hubbard model that features a competition between charge density correlations and pairing correlations, and show that the DOB method successfully identifies their relative strengths and dominant correlators. To calculate the ground state of this model, we use the density matrix renormalization group, formulated in terms of a variational matrix product state (MPS) approach within which subsequent determination of the CDM is very straightforward. In an extended appendix, we give a detailed tutorial introduction to our variational MPS approach for ground state calculations for one-dimensional quantum chain models. We present in detail how MPSs overcome the problem of large Hilbert space dimensions in these models and describe all the techniques needed for handling them in practice.

Gamow-Jordan vectors and non-reducible density operators from higher-order S-matrix poles

International Nuclear Information System (INIS)

Bohm, A.; Loewe, M.; Maxson, S.; Patuleanu, P.; Puentmann, C.; Gadella, M.

1997-01-01

In analogy to Gamow vectors that are obtained from first-order resonance poles of the S-matrix, one can also define higher-order Gamow vectors which are derived from higher-order poles of the S-matrix. An S-matrix pole of r-th order at z R =E R -iΓ/2 leads to r generalized eigenvectors of order k=0,1,hor-ellipsis,r-1, which are also Jordan vectors of degree (k+1) with generalized eigenvalue (E R -iΓ/2). The Gamow-Jordan vectors are elements of a generalized complex eigenvector expansion, whose form suggests the definition of a state operator (density matrix) for the microphysical decaying state of this higher-order pole. This microphysical state is a mixture of non-reducible components. In spite of the fact that the k-th order Gamow-Jordan vectors has the polynomial time-dependence which one always associates with higher-order poles, the microphysical state obeys a purely exponential decay law. copyright 1997 American Institute of Physics

Gradient-based stochastic estimation of the density matrix

Science.gov (United States)

Wang, Zhentao; Chern, Gia-Wei; Batista, Cristian D.; Barros, Kipton

2018-03-01

Fast estimation of the single-particle density matrix is key to many applications in quantum chemistry and condensed matter physics. The best numerical methods leverage the fact that the density matrix elements f(H)ij decay rapidly with distance rij between orbitals. This decay is usually exponential. However, for the special case of metals at zero temperature, algebraic decay of the density matrix appears and poses a significant numerical challenge. We introduce a gradient-based probing method to estimate all local density matrix elements at a computational cost that scales linearly with system size. For zero-temperature metals, the stochastic error scales like S-(d+2)/2d, where d is the dimension and S is a prefactor to the computational cost. The convergence becomes exponential if the system is at finite temperature or is insulating.

Explicit treatment of N-body correlations within a density-matrix formalism

International Nuclear Information System (INIS)

Shun-Jin, W.; Cassing, W.

1985-01-01

The nuclear many-body problem is reformulated in the density-matrix approach such that n-body correlations are separated out from the reduced density matrix rho/sub n/. A set of equations for the time evolution of the n-body correlations c/sub n/ is derived which allows for physically transparent truncations with respect to the order of correlations. In the stationary limit (c/sub n/ = 0) a restriction to two-body correlations yields a generalized Bethe-Goldstone equation a restriction to body correlations yields generalized Faddeev equations in the density-matrix formulation. Furthermore it can be shown that any truncation of the set of equations (c/sub n/ = 0, n>m) is compatible with conservation laws, a quality which in general is not fulfilled if higher order correlations are treated perturbatively

Correlated density matrix theory of spatially inhomogeneous Bose fluids

International Nuclear Information System (INIS)

Gernoth, K.A.; Clark, J.W.; Ristig, M.L.

1994-06-01

In this paper, the variational Hartree-Jastrow theory of the ground state of spatially inhomogeneous Bose systems is extended to finite temperatures. The theory presented here is a generalization also in the sense that it extends the correlated density matrix approach, formulated previously for uniform Bose fluids, to systems with nonuniform density profiles. The method provides a framework in which the effects of thermal excitations on the spatial structure of a Bose fluid, as represented by the density profile and the two-body distribution functions, may be discussed on the basis on an ab initio microscopic description of the system. Thermal excitations make their appearance through self-consistently determined one-body and two-body potentials which enter the nonlinear, coupled Euler-Lagrange equations for the one-body density and for the pair distribution function. Since back-flow correlations are neglected, the excitations are described by a Feynman eigenvalue equation, suitably generalized to nonzero temperatures. The only external quantities entering the correlated density matrix theory elaborated here are the bare two-body interaction potential and, in actual applications, the boundary conditions to be imposed on the one-body density. 30 refs

Computing the Density Matrix in Electronic Structure Theory on Graphics Processing Units.

Science.gov (United States)

Cawkwell, M J; Sanville, E J; Mniszewski, S M; Niklasson, Anders M N

2012-11-13

The self-consistent solution of a Schrödinger-like equation for the density matrix is a critical and computationally demanding step in quantum-based models of interatomic bonding. This step was tackled historically via the diagonalization of the Hamiltonian. We have investigated the performance and accuracy of the second-order spectral projection (SP2) algorithm for the computation of the density matrix via a recursive expansion of the Fermi operator in a series of generalized matrix-matrix multiplications. We demonstrate that owing to its simplicity, the SP2 algorithm [Niklasson, A. M. N. Phys. Rev. B2002, 66, 155115] is exceptionally well suited to implementation on graphics processing units (GPUs). The performance in double and single precision arithmetic of a hybrid GPU/central processing unit (CPU) and full GPU implementation of the SP2 algorithm exceed those of a CPU-only implementation of the SP2 algorithm and traditional matrix diagonalization when the dimensions of the matrices exceed about 2000 × 2000. Padding schemes for arrays allocated in the GPU memory that optimize the performance of the CUBLAS implementations of the level 3 BLAS DGEMM and SGEMM subroutines for generalized matrix-matrix multiplications are described in detail. The analysis of the relative performance of the hybrid CPU/GPU and full GPU implementations indicate that the transfer of arrays between the GPU and CPU constitutes only a small fraction of the total computation time. The errors measured in the self-consistent density matrices computed using the SP2 algorithm are generally smaller than those measured in matrices computed via diagonalization. Furthermore, the errors in the density matrices computed using the SP2 algorithm do not exhibit any dependence of system size, whereas the errors increase linearly with the number of orbitals when diagonalization is employed.

Collective excitations in semiconductor superlattices and plasma modes of a two-dimensional electron gas with spatially modulated charge density

International Nuclear Information System (INIS)

Eliasson, G.L.

1987-01-01

The theory of collective excitations in semiconductor superlattices is formulated by using linear response theory. Different kinds of collective excitations in type I (GaAs-GaAlAs) and type II (GaSb-InAs) superlattices are surveyed. Special attention is paid to the presence of surface and finite-size effects. In calculating the dielectric matrix, the effect of different approximations of the system is discussed. The theory for inelastic length scattering (Raman scattering), and for Electron Energy Loss (EEL) due to collective excitations, is formulated. Calculations for several model systems are presented and the main features of the spectra are discussed. In part II the theory of collective excitations of a two-dimensional electron gas with a spatially periodic equilibrium density is formulated. As a first example a periodic array of two-dimensional electron gas strips with constant equilibrium density is studied. The integral equation that describes the charge fluctuations on the strips is derived and solved numerically. The spatial dependence of the density fluctuation across a single strip can be in the form of either propagating or evanescent waves

Variational Optimization of the Second-Order Density Matrix Corresponding to a Seniority-Zero Configuration Interaction Wave Function.

Science.gov (United States)

Poelmans, Ward; Van Raemdonck, Mario; Verstichel, Brecht; De Baerdemacker, Stijn; Torre, Alicia; Lain, Luis; Massaccesi, Gustavo E; Alcoba, Diego R; Bultinck, Patrick; Van Neck, Dimitri

2015-09-08

We perform a direct variational determination of the second-order (two-particle) density matrix corresponding to a many-electron system, under a restricted set of the two-index N-representability P-, Q-, and G-conditions. In addition, we impose a set of necessary constraints that the two-particle density matrix must be derivable from a doubly occupied many-electron wave function, i.e., a singlet wave function for which the Slater determinant decomposition only contains determinants in which spatial orbitals are doubly occupied. We rederive the two-index N-representability conditions first found by Weinhold and Wilson and apply them to various benchmark systems (linear hydrogen chains, He, N2, and CN(-)). This work is motivated by the fact that a doubly occupied many-electron wave function captures in many cases the bulk of the static correlation. Compared to the general case, the structure of doubly occupied two-particle density matrices causes the associate semidefinite program to have a very favorable scaling as L(3), where L is the number of spatial orbitals. Since the doubly occupied Hilbert space depends on the choice of the orbitals, variational calculation steps of the two-particle density matrix are interspersed with orbital-optimization steps (based on Jacobi rotations in the space of the spatial orbitals). We also point to the importance of symmetry breaking of the orbitals when performing calculations in a doubly occupied framework.

Matrix-product-state method with local basis optimization for nonequilibrium electron-phonon systems

Science.gov (United States)

Heidrich-Meisner, Fabian; Brockt, Christoph; Dorfner, Florian; Vidmar, Lev; Jeckelmann, Eric

We present a method for simulating the time evolution of quasi-one-dimensional correlated systems with strongly fluctuating bosonic degrees of freedom (e.g., phonons) using matrix product states. For this purpose we combine the time-evolving block decimation (TEBD) algorithm with a local basis optimization (LBO) approach. We discuss the performance of our approach in comparison to TEBD with a bare boson basis, exact diagonalization, and diagonalization in a limited functional space. TEBD with LBO can reduce the computational cost by orders of magnitude when boson fluctuations are large and thus it allows one to investigate problems that are out of reach of other approaches. First, we test our method on the non-equilibrium dynamics of a Holstein polaron and show that it allows us to study the regime of strong electron-phonon coupling. Second, the method is applied to the scattering of an electronic wave packet off a region with electron-phonon coupling. Our study reveals a rich physics including transient self-trapping and dissipation. Supported by Deutsche Forschungsgemeinschaft (DFG) via FOR 1807.

Three-body forces for electrons by the S-matrix method

International Nuclear Information System (INIS)

Margaritelli, R.

1989-01-01

A electromagnetic three-body potential between eletrons is derived by the S-matrix method. This potential can be compared up to a certain point with other electromagnetic potentials (obtained by other methods) encountered in the literature. However, since the potential derived here is far more complete than others, this turns direct comparison with the potentials found in the literature somewhat difficult. These calculations allow a better understanding of the S-matrix method as applied to problems which involve the calculations of three-body nuclear forces (these calculations are performed in order to understand the 3 He form factor). Furthermore, these results enable us to decide between two discrepant works which derive the two-pion exchange three-body potential, both by the S-matrix method. (author) [pt

Relative Contribution of Matrix Structure, Patch Resources and Management to the Local Densities of Two Large Blue Butterfly Species.

Science.gov (United States)

Kajzer-Bonk, Joanna; Skórka, Piotr; Nowicki, Piotr; Bonk, Maciej; Król, Wiesław; Szpiłyk, Damian; Woyciechowski, Michal

2016-01-01

The type of matrix, the landscape surrounding habitat patches, may determine the distribution and function of local populations. However, the matrix is often heterogeneous, and its various components may differentially contribute to metapopulation processes at different spatial scales, a phenomenon that has rarely been investigated. The aim of this study was to estimate the relative importance of matrix composition and spatial scale, habitat quality, and management intensity on the occurrence and density of local populations of two endangered large blue butterflies: Phengaris teleius and P. nausithous. Presence and abundance data were assessed over two years, 2011-12, in 100 local patches within two heterogeneous regions (near Kraków and Tarnów, southern Poland). The matrix composition was analyzed at eight spatial scales. We observed high occupancy rates in both species, regions and years. With the exception of area and isolation, almost all of the matrix components contributed to Phengaris sp. densities. The different matrix components acted at different spatial scales (grassland cover within 4 and 3 km, field cover within 0.4 and 0.3 km and water cover within 4 km radii for P. teleius and P. nausithous, respectively) and provided the highest independent contribution to the butterfly densities. Additionally, the effects of a 0.4 km radius of forest cover and a food plant cover on P. teleius, and a 1 km radius of settlement cover and management intensity on P. nausithous densities were observed. Contrary to former studies we conclude that the matrix heterogeneity and spatial scale rather than general matrix type are of relevance for densities of butterflies. Conservation strategies for these umbrella species should concentrate on maintaining habitat quality and managing matrix composition at the most appropriate spatial scales.

Relative Contribution of Matrix Structure, Patch Resources and Management to the Local Densities of Two Large Blue Butterfly Species

Science.gov (United States)

Skórka, Piotr; Nowicki, Piotr; Bonk, Maciej; Król, Wiesław; Szpiłyk, Damian; Woyciechowski, Michal

2016-01-01

The type of matrix, the landscape surrounding habitat patches, may determine the distribution and function of local populations. However, the matrix is often heterogeneous, and its various components may differentially contribute to metapopulation processes at different spatial scales, a phenomenon that has rarely been investigated. The aim of this study was to estimate the relative importance of matrix composition and spatial scale, habitat quality, and management intensity on the occurrence and density of local populations of two endangered large blue butterflies: Phengaris teleius and P. nausithous. Presence and abundance data were assessed over two years, 2011–12, in 100 local patches within two heterogeneous regions (near Kraków and Tarnów, southern Poland). The matrix composition was analyzed at eight spatial scales. We observed high occupancy rates in both species, regions and years. With the exception of area and isolation, almost all of the matrix components contributed to Phengaris sp. densities. The different matrix components acted at different spatial scales (grassland cover within 4 and 3 km, field cover within 0.4 and 0.3 km and water cover within 4 km radii for P. teleius and P. nausithous, respectively) and provided the highest independent contribution to the butterfly densities. Additionally, the effects of a 0.4 km radius of forest cover and a food plant cover on P. teleius, and a 1 km radius of settlement cover and management intensity on P. nausithous densities were observed. Contrary to former studies we conclude that the matrix heterogeneity and spatial scale rather than general matrix type are of relevance for densities of butterflies. Conservation strategies for these umbrella species should concentrate on maintaining habitat quality and managing matrix composition at the most appropriate spatial scales. PMID:28005942

Extending the range of real time density matrix renormalization group simulations

Science.gov (United States)

Kennes, D. M.; Karrasch, C.

2016-03-01

We discuss a few simple modifications to time-dependent density matrix renormalization group (DMRG) algorithms which allow to access larger time scales. We specifically aim at beginners and present practical aspects of how to implement these modifications within any standard matrix product state (MPS) based formulation of the method. Most importantly, we show how to 'combine' the Schrödinger and Heisenberg time evolutions of arbitrary pure states | ψ 〉 and operators A in the evaluation of 〈A〉ψ(t) = 〈 ψ | A(t) | ψ 〉 . This includes quantum quenches. The generalization to (non-)thermal mixed state dynamics 〈A〉ρ(t) =Tr [ ρA(t) ] induced by an initial density matrix ρ is straightforward. In the context of linear response (ground state or finite temperature T > 0) correlation functions, one can extend the simulation time by a factor of two by 'exploiting time translation invariance', which is efficiently implementable within MPS DMRG. We present a simple analytic argument for why a recently-introduced disentangler succeeds in reducing the effort of time-dependent simulations at T > 0. Finally, we advocate the python programming language as an elegant option for beginners to set up a DMRG code.

Relativistic density functional theory with picture-change corrected electron density based on infinite-order Douglas-Kroll-Hess method

Science.gov (United States)

Oyama, Takuro; Ikabata, Yasuhiro; Seino, Junji; Nakai, Hiromi

2017-07-01

This Letter proposes a density functional treatment based on the two-component relativistic scheme at the infinite-order Douglas-Kroll-Hess (IODKH) level. The exchange-correlation energy and potential are calculated using the electron density based on the picture-change corrected density operator transformed by the IODKH method. Numerical assessments indicated that the picture-change uncorrected density functional terms generate significant errors, on the order of hartree for heavy atoms. The present scheme was found to reproduce the energetics in the four-component treatment with high accuracy.

Time-Dependent Density Functional Theory for Open Systems and Its Applications.

Science.gov (United States)

Chen, Shuguang; Kwok, YanHo; Chen, GuanHua

2018-02-20

Photovoltaic devices, electrochemical cells, catalysis processes, light emitting diodes, scanning tunneling microscopes, molecular electronics, and related devices have one thing in common: open quantum systems where energy and matter are not conserved. Traditionally quantum chemistry is confined to isolated and closed systems, while quantum dissipation theory studies open quantum systems. The key quantity in quantum dissipation theory is the reduced system density matrix. As the reduced system density matrix is an O(M! × M!) matrix, where M is the number of the particles of the system of interest, quantum dissipation theory can only be employed to simulate systems of a few particles or degrees of freedom. It is thus important to combine quantum chemistry and quantum dissipation theory so that realistic open quantum systems can be simulated from first-principles. We have developed a first-principles method to simulate the dynamics of open electronic systems, the time-dependent density functional theory for open systems (TDDFT-OS). Instead of the reduced system density matrix, the key quantity is the reduced single-electron density matrix, which is an N × N matrix where N is the number of the atomic bases of the system of interest. As the dimension of the key quantity is drastically reduced, the TDDFT-OS can thus be used to simulate the dynamics of realistic open electronic systems and efficient numerical algorithms have been developed. As an application, we apply the method to study how quantum interference develops in a molecular transistor in time domain. We include electron-phonon interaction in our simulation and show that quantum interference in the given system is robust against nuclear vibration not only in the steady state but also in the transient dynamics. As another application, by combining TDDFT-OS with Ehrenfest dynamics, we study current-induced dissociation of water molecules under scanning tunneling microscopy and follow its time dependent

Path integral density matrix dynamics: a method for calculating time-dependent properties in thermal adiabatic and non-adiabatic systems.

Science.gov (United States)

Habershon, Scott

2013-09-14

We introduce a new approach for calculating quantum time-correlation functions and time-dependent expectation values in many-body thermal systems; both electronically adiabatic and non-adiabatic cases can be treated. Our approach uses a path integral simulation to sample an initial thermal density matrix; subsequent evolution of this density matrix is equivalent to solution of the time-dependent Schrödinger equation, which we perform using a linear expansion of Gaussian wavepacket basis functions which evolve according to simple classical-like trajectories. Overall, this methodology represents a formally exact approach for calculating time-dependent quantum properties; by introducing approximations into both the imaginary-time and real-time propagations, this approach can be adapted for complex many-particle systems interacting through arbitrary potentials. We demonstrate this method for the spin Boson model, where we find good agreement with numerically exact calculations. We also discuss future directions of improvement for our approach with a view to improving accuracy and efficiency.

Path integral density matrix dynamics: A method for calculating time-dependent properties in thermal adiabatic and non-adiabatic systems

International Nuclear Information System (INIS)

Habershon, Scott

2013-01-01

We introduce a new approach for calculating quantum time-correlation functions and time-dependent expectation values in many-body thermal systems; both electronically adiabatic and non-adiabatic cases can be treated. Our approach uses a path integral simulation to sample an initial thermal density matrix; subsequent evolution of this density matrix is equivalent to solution of the time-dependent Schrödinger equation, which we perform using a linear expansion of Gaussian wavepacket basis functions which evolve according to simple classical-like trajectories. Overall, this methodology represents a formally exact approach for calculating time-dependent quantum properties; by introducing approximations into both the imaginary-time and real-time propagations, this approach can be adapted for complex many-particle systems interacting through arbitrary potentials. We demonstrate this method for the spin Boson model, where we find good agreement with numerically exact calculations. We also discuss future directions of improvement for our approach with a view to improving accuracy and efficiency

The ab-initio density matrix renormalization group in practice

Energy Technology Data Exchange (ETDEWEB)

Olivares-Amaya, Roberto; Hu, Weifeng; Sharma, Sandeep; Yang, Jun; Chan, Garnet Kin-Lic [Department of Chemistry, Princeton University, Princeton, New Jersey 08544 (United States); Nakatani, Naoki [Department of Chemistry, Princeton University, Princeton, New Jersey 08544 (United States); Catalysis Research Center, Hokkaido University, Kita 21 Nishi 10, Sapporo, Hokkaido 001-0021 (Japan)

2015-01-21

The ab-initio density matrix renormalization group (DMRG) is a tool that can be applied to a wide variety of interesting problems in quantum chemistry. Here, we examine the density matrix renormalization group from the vantage point of the quantum chemistry user. What kinds of problems is the DMRG well-suited to? What are the largest systems that can be treated at practical cost? What sort of accuracies can be obtained, and how do we reason about the computational difficulty in different molecules? By examining a diverse benchmark set of molecules: π-electron systems, benchmark main-group and transition metal dimers, and the Mn-oxo-salen and Fe-porphine organometallic compounds, we provide some answers to these questions, and show how the density matrix renormalization group is used in practice.

The ab-initio density matrix renormalization group in practice.

Science.gov (United States)

Olivares-Amaya, Roberto; Hu, Weifeng; Nakatani, Naoki; Sharma, Sandeep; Yang, Jun; Chan, Garnet Kin-Lic

2015-01-21

The ab-initio density matrix renormalization group (DMRG) is a tool that can be applied to a wide variety of interesting problems in quantum chemistry. Here, we examine the density matrix renormalization group from the vantage point of the quantum chemistry user. What kinds of problems is the DMRG well-suited to? What are the largest systems that can be treated at practical cost? What sort of accuracies can be obtained, and how do we reason about the computational difficulty in different molecules? By examining a diverse benchmark set of molecules: π-electron systems, benchmark main-group and transition metal dimers, and the Mn-oxo-salen and Fe-porphine organometallic compounds, we provide some answers to these questions, and show how the density matrix renormalization group is used in practice.

Method of measuring the current density distribution and emittance of pulsed electron beams

International Nuclear Information System (INIS)

Schilling, H.B.

1979-07-01

This method of current density measurement employs an array of many Faraday cups, each cup being terminated by an integrating capacitor. The voltages of the capacitors are subsequently displayed on a scope, thus giving the complete current density distribution with one shot. In the case of emittance measurements, a moveable small-diameter aperture is inserted at some distance in front of the cup array. Typical results with a two-cathode, two-energy electron source are presented. (orig.)

Quantum Stochastic Trajectories: The Fokker-Planck-Bohm Equation Driven by the Reduced Density Matrix.

Science.gov (United States)

Avanzini, Francesco; Moro, Giorgio J

2018-03-15

The quantum molecular trajectory is the deterministic trajectory, arising from the Bohm theory, that describes the instantaneous positions of the nuclei of molecules by assuring the agreement with the predictions of quantum mechanics. Therefore, it provides the suitable framework for representing the geometry and the motions of molecules without neglecting their quantum nature. However, the quantum molecular trajectory is extremely demanding from the computational point of view, and this strongly limits its applications. To overcome such a drawback, we derive a stochastic representation of the quantum molecular trajectory, through projection operator techniques, for the degrees of freedom of an open quantum system. The resulting Fokker-Planck operator is parametrically dependent upon the reduced density matrix of the open system. Because of the pilot role played by the reduced density matrix, this stochastic approach is able to represent accurately the main features of the open system motions both at equilibrium and out of equilibrium with the environment. To verify this procedure, the predictions of the stochastic and deterministic representation are compared for a model system of six interacting harmonic oscillators, where one oscillator is taken as the open quantum system of interest. The undeniable advantage of the stochastic approach is that of providing a simplified and self-contained representation of the dynamics of the open system coordinates. Furthermore, it can be employed to study the out of equilibrium dynamics and the relaxation of quantum molecular motions during photoinduced processes, like photoinduced conformational changes and proton transfers.

Evaluation of the thermodynamics of a four level system using canonical density matrix method

Directory of Open Access Journals (Sweden)

Awoga Oladunjoye A.

2013-02-01

Full Text Available We consider a four-level system with two subsystems coupled by weak interaction. The system is in thermal equilibrium. The thermodynamics of the system, namely internal energy, free energy, entropy and heat capacity, are evaluated using the canonical density matrix by two methods. First by Kronecker product method and later by treating the subsystems separately and then adding the evaluated thermodynamic properties of each subsystem. It is discovered that both methods yield the same result, the results obey the laws of thermodynamics and are the same as earlier obtained results. The results also show that each level of the subsystems introduces a new degree of freedom and increases the entropy of the entire system. We also found that the four-level system predicts a linear relationship between heat capacity and temperature at very low temperatures just as in metals. Our numerical results show the same trend.

Matrix method for two-dimensional waveguide mode solution

Science.gov (United States)

Sun, Baoguang; Cai, Congzhong; Venkatesh, Balajee Seshasayee

2018-05-01

In this paper, we show that the transfer matrix theory of multilayer optics can be used to solve the modes of any two-dimensional (2D) waveguide for their effective indices and field distributions. A 2D waveguide, even composed of numerous layers, is essentially a multilayer stack and the transmission through the stack can be analysed using the transfer matrix theory. The result is a transfer matrix with four complex value elements, namely A, B, C and D. The effective index of a guided mode satisfies two conditions: (1) evanescent waves exist simultaneously in the first (cladding) layer and last (substrate) layer, and (2) the complex element D vanishes. For a given mode, the field distribution in the waveguide is the result of a 'folded' plane wave. In each layer, there is only propagation and absorption; at each boundary, only reflection and refraction occur, which can be calculated according to the Fresnel equations. As examples, we show that this method can be used to solve modes supported by the multilayer step-index dielectric waveguide, slot waveguide, gradient-index waveguide and various plasmonic waveguides. The results indicate the transfer matrix method is effective for 2D waveguide mode solution in general.

Metal-insulator transition in disordered systems from the one-body density matrix

DEFF Research Database (Denmark)

Olsen, Thomas; Resta, Raffaele; Souza, Ivo

2017-01-01

The insulating state of matter can be probed by means of a ground state geometrical marker, which is closely related to the modern theory of polarization (based on a Berry phase). In the present work we show that this marker can be applied to determine the metal-insulator transition in disordered...... the one-body density matrix. The approach has a general ab initio formulation and could in principle be applied to realistic disordered materials by standard electronic structure methods....... systems. In particular, for noninteracting systems the geometrical marker can be obtained from the configurational average of the norm-squared one-body density matrix, which can be calculated within open as well as periodic boundary conditions. This is in sharp contrast to a classification based...

The density matrix - The story of a failed transfer

Energy Technology Data Exchange (ETDEWEB)

Blum, Alexander [MPI fuer Wissenschaftsgeschichte, Berlin (Germany)

2013-07-01

With the discovery of the positron in 1933, Paul Dirac (along with most other physicists) was forced to really take seriously his earlier suggestion that in the world as we know it all negative energy states are occupied and we are thus surrounded by an infinite sea of electrons. What was needed was a way to treat this large number of electrons in a manageable fashion. Dirac resorted to the use of the density matrix, a technique he had earlier used to describe the large number of electrons in complex atoms. Initially, this transfer from atomic physics to what we would nowadays call particle physics was quite successful, and for a few years the density matrix was the state of the art in describing the Dirac electron sea, but then rapidly fell out of favor. I investigate the causes of this ultimately failed transfer and how it relates to changes in the physical notion of the vacuum, changes which eventually eliminated the analogy on which the transfer had been based in the first place.

Large-distance and long-time asymptotic behavior of the reduced density matrix in the non-linear Schroedinger model

Energy Technology Data Exchange (ETDEWEB)

Kozlowski, K.K.

2010-12-15

Starting from the form factor expansion in finite volume, we derive the multidimensional generalization of the so-called Natte series for the zero-temperature, time and distance dependent reduced density matrix in the non-linear Schroedinger model. This representation allows one to read-off straightforwardly the long-time/large-distance asymptotic behavior of this correlator. Our method of analysis reduces the complexity of the computation of the asymptotic behavior of correlation functions in the so-called interacting integrable models, to the one appearing in free fermion equivalent models. We compute explicitly the first few terms appearing in the asymptotic expansion. Part of these terms stems from excitations lying away from the Fermi boundary, and hence go beyond what can be obtained by using the CFT/Luttinger liquid based predictions. (orig.)

Variational and robust density fitting of four-center two-electron integrals in local metrics

Science.gov (United States)

Reine, Simen; Tellgren, Erik; Krapp, Andreas; Kjærgaard, Thomas; Helgaker, Trygve; Jansik, Branislav; Høst, Stinne; Salek, Paweł

2008-09-01

Density fitting is an important method for speeding up quantum-chemical calculations. Linear-scaling developments in Hartree-Fock and density-functional theories have highlighted the need for linear-scaling density-fitting schemes. In this paper, we present a robust variational density-fitting scheme that allows for solving the fitting equations in local metrics instead of the traditional Coulomb metric, as required for linear scaling. Results of fitting four-center two-electron integrals in the overlap and the attenuated Gaussian damped Coulomb metric are presented, and we conclude that density fitting can be performed in local metrics at little loss of chemical accuracy. We further propose to use this theory in linear-scaling density-fitting developments.

Hydride vapor phase GaN films with reduced density of residual electrons and deep traps

International Nuclear Information System (INIS)

Polyakov, A. Y.; Smirnov, N. B.; Govorkov, A. V.; Yugova, T. G.; Cox, H.; Helava, H.; Makarov, Yu.; Usikov, A. S.

2014-01-01

Electrical properties and deep electron and hole traps spectra are compared for undoped n-GaN films grown by hydride vapor phase epitaxy (HVPE) in the regular process (standard HVPE samples) and in HVPE process optimized for decreasing the concentration of residual donor impurities (improved HVPE samples). It is shown that the residual donor density can be reduced by optimization from ∼10 17  cm −3 to (2–5) × 10 14  cm −3 . The density of deep hole traps and deep electron traps decreases with decreased donor density, so that the concentration of deep hole traps in the improved samples is reduced to ∼5 × 10 13  cm −3 versus 2.9 × 10 16  cm −3 in the standard samples, with a similar decrease in the electron traps concentration

Matrix product operators, matrix product states, and ab initio density matrix renormalization group algorithms

Science.gov (United States)

Chan, Garnet Kin-Lic; Keselman, Anna; Nakatani, Naoki; Li, Zhendong; White, Steven R.

2016-07-01

Current descriptions of the ab initio density matrix renormalization group (DMRG) algorithm use two superficially different languages: an older language of the renormalization group and renormalized operators, and a more recent language of matrix product states and matrix product operators. The same algorithm can appear dramatically different when written in the two different vocabularies. In this work, we carefully describe the translation between the two languages in several contexts. First, we describe how to efficiently implement the ab initio DMRG sweep using a matrix product operator based code, and the equivalence to the original renormalized operator implementation. Next we describe how to implement the general matrix product operator/matrix product state algebra within a pure renormalized operator-based DMRG code. Finally, we discuss two improvements of the ab initio DMRG sweep algorithm motivated by matrix product operator language: Hamiltonian compression, and a sum over operators representation that allows for perfect computational parallelism. The connections and correspondences described here serve to link the future developments with the past and are important in the efficient implementation of continuing advances in ab initio DMRG and related algorithms.

Measuring the electron density in plasmas from the difference of Lorentzian part of the widths of two Balmer series hydrogen lines

Energy Technology Data Exchange (ETDEWEB)

Yubero, C. [Grupo de Física de Plasmas: Diagnosis, Modelos y Aplicaciones (FQM-136), Edificio A. Einstein (C-2), Campus de Rabanales, Universidad de Córdoba, 14071 Córdoba (Spain); García, M.C., E-mail: fa1gamam@uco.es [Grupo de Física de Plasmas: Diagnosis, Modelos y Aplicaciones (FQM-136), Edificio A. Einstein (C-2), Campus de Rabanales, Universidad de Córdoba, 14071 Córdoba (Spain); Dimitrijevic, M.S. [Astronomical Observatory, Volgina 7, 11060 Belgrade (Serbia); Sola, A.; Gamero, A. [Grupo de Física de Plasmas: Diagnosis, Modelos y Aplicaciones (FQM-136), Edificio A. Einstein (C-2), Campus de Rabanales, Universidad de Córdoba, 14071 Córdoba (Spain)

2015-05-01

We present an alternative optical emission spectroscopy method to measure the plasma electron density from the difference of widths of two Balmer series hydrogen lines (H{sub α} and H{sub β}), especially convenient for non-thermal plasmas since with this method, there is no need to know either the gas temperature or the van der Waals contribution to the Lorentzian part of the line. In this paper it has been assumed that the part of full width at half maximum due to Stark broadening can be determined with the approximation of Lorentzian line shape. The method has been applied to the determination of the electron density in an argon microwave-induced plasma maintained at atmospheric pressure, and comparison with the results obtained using other diagnostic methods has been done. - Highlights: • An alternative method to measure the electron density in plasmas from two Balmer series hydrogen lines (H{sub α} and H{sub β}) is presented. • The method is very convenient for plasmas with electron densities of the order of 10{sup 14} cm{sup −3} and above, at low gas temperatures. • It has been applied to the determination of the electron density of an argon microwave plasma at atmospheric pressure. • Results from it are in good agreement with previous ones obtained using other diagnostic methods.

Apparatus and method for generating high density pulses of electrons

International Nuclear Information System (INIS)

Lee, C.; Oettinger, P.E.

1981-01-01

An apparatus and method are described for the production of high density pulses of electrons using a laser energized emitter. Caesium atoms from a low pressure vapour atmosphere are absorbed on and migrate from a metallic target rapidly heated by a laser to a high temperature. Due to this heating time being short compared with the residence time of the caesium atoms adsorbed on the target surface, copious electrons are emitted which form a high current density pulse. (U.K.)

Study of electron-molecule collision via finite-element method and r-matrix propagation technique: Exact exchange

International Nuclear Information System (INIS)

Abdolsalami, F.; Abdolsalami, M.; Perez, L.; Gomez, P.

1995-01-01

The authors have applied the finite-element method to electron-molecule collision with the exchange effect implemented rigorously. All the calculations are done in the body-frame within the fixed-nuclei approximation, where the exact treatment of exchange as a nonlocal effect results in a set of coupled integro-differential equations. The method is applied to e-H 2 and e-N 2 scatterings and the cross sections obtained are in very good agreement with the corresponding results the authors have generated from the linear-algebraic approach. This confirms the significant difference observed between their results generated by linear-algebraic method and the previously published e-N 2 cross sections. Their studies show that the finite-element method is clearly superior to the linear-algebraic approach in both memory usage and CPU time especially for large systems such as e-N 2 . The system coefficient matrix obtained from the finite-element method is often sparse and smaller in size by a factor of 12 to 16, compared to the linear-algebraic technique. Moreover, the CPU time required to obtain stable results with the finite-element method is significantly smaller than the linear-algebraic approach for one incident electron energy. The usage of computer resources in the finite-element method can even be reduced much further when (1) scattering calculations involving multiple electron energies are performed in one computer run and (2) exchange, which is a short range effect, is approximated by a sparse matrix. 17 refs., 7 figs., 5 tabs

Time-dependent quantum many-body systems. Linear response, electronic transport, and reduced density matrices

International Nuclear Information System (INIS)

Appel, H.

2007-05-01

In part I of this work we present a double-pole approximation (DPA) to the response equations of time-dependent density functional theory (TDDFT). The double-pole approximation provides an exact description of systems with two strongly coupled excitations which are isolated from the rest of the spectrum. In contrast to the traditional single-pole approximation of TDDFT the DPA also yields corrections to the Kohn-Sham oscillator strengths. We also demonstrate how to invert the double-pole solution which allows us to predict matrix elements of the exchange-correlation kernel f xc from experimental input. We attempt some first steps towards a time-dependent generalization of reduced density matrix functional theory (RDMFT). In part II we derive equations of motion for natural orbitals and occupation numbers. Using the equation of motion for the occupation numbers we show that an adiabatic extension of presently known ground-state functionals of static RDMFT always leads to occupation numbers which are constant in time. From the stationary conditions of the equations of motion for the N-body correlations (correlated parts of the N-body matrices) we derive a new class of ground-state functionals which can be used in static RDMFT. Applications are presented for a one-dimensional model system where the time-dependent many-body Schroedinger equation can be propagated numerically. We use optimal control theory to find optimized laser pulses for transitions in a model for atomic Helium. From the numerically exact correlated wavefunction we extract the exact time evolution of natural orbitals and occupation numbers for (i) laser-driven Helium and (ii) electron-ion scattering. Part III of this work considers time-dependent quantum transport within TDDFT. We present an algorithm for the calculation of extended eigenstates of single-particle Hamiltonians which is especially tailored to a finite-difference discretization of the Schroedinger equation. We consider the propagation

Time-dependent quantum many-body systems. Linear response, electronic transport, and reduced density matrices

Energy Technology Data Exchange (ETDEWEB)

Appel, H.

2007-05-15

In part I of this work we present a double-pole approximation (DPA) to the response equations of time-dependent density functional theory (TDDFT). The double-pole approximation provides an exact description of systems with two strongly coupled excitations which are isolated from the rest of the spectrum. In contrast to the traditional single-pole approximation of TDDFT the DPA also yields corrections to the Kohn-Sham oscillator strengths. We also demonstrate how to invert the double-pole solution which allows us to predict matrix elements of the exchange-correlation kernel f{sub xc} from experimental input. We attempt some first steps towards a time-dependent generalization of reduced density matrix functional theory (RDMFT). In part II we derive equations of motion for natural orbitals and occupation numbers. Using the equation of motion for the occupation numbers we show that an adiabatic extension of presently known ground-state functionals of static RDMFT always leads to occupation numbers which are constant in time. From the stationary conditions of the equations of motion for the N-body correlations (correlated parts of the N-body matrices) we derive a new class of ground-state functionals which can be used in static RDMFT. Applications are presented for a one-dimensional model system where the time-dependent many-body Schroedinger equation can be propagated numerically. We use optimal control theory to find optimized laser pulses for transitions in a model for atomic Helium. From the numerically exact correlated wavefunction we extract the exact time evolution of natural orbitals and occupation numbers for (i) laser-driven Helium and (ii) electron-ion scattering. Part III of this work considers time-dependent quantum transport within TDDFT. We present an algorithm for the calculation of extended eigenstates of single-particle Hamiltonians which is especially tailored to a finite-difference discretization of the Schroedinger equation. We consider the

Critical density for Landau damping in a two-electron-component plasma

Energy Technology Data Exchange (ETDEWEB)

Rupp, Constantin F.; López, Rodrigo A.; Araneda, Jaime A. [Departamento de Física, Facultad de Ciencias Físicas y Matemáticas, Universidad de Concepción, Concepción (Chile)

2015-10-15

The asymptotic evolution of an initial perturbation in a collisionless two-electron-component plasma with different temperatures is studied numerically. The transition between linear and nonlinear damping regimes is determined by slowly varying the density of the secondary electron-component using high-resolution Vlasov-Poisson simulations. It is shown that, for fixed amplitude perturbations, this transition behaves as a critical phenomenon with time scales and field amplitudes exhibiting power-law dependencies on the threshold density, similar to the critical amplitude behavior in a single-component plasma.

Potential and electron density calculated for freely expanding plasma by an electron beam

International Nuclear Information System (INIS)

Ho, C. Y.; Tsai, Y. H.; Ma, C.; Wen, M. Y.

2011-01-01

This paper investigates the radial distributions of potential and electron density in free expansion plasma induced by an electron beam irradiating on the plate. The region of plasma production is assumed to be cylindrical, and the plasma expansion is assumed to be from a cylindrical source. Therefore, the one-dimensional model in cylindrical coordinates is employed in order to analyze the radial distributions of the potential and electron density. The Runge-Kutta method and the perturbation method are utilized in order to obtain the numerical and approximate solutions, respectively. The results reveal that the decrease in the initial ion energy makes most of the ions gather near the plasma production region and reduces the distribution of the average positive potential, electron, and ion density along the radial direction. The oscillation of steady-state plasma along the radial direction is also presented in this paper. The ions induce a larger amplitude of oscillation along the radial direction than do electrons because the electrons oscillate around slowly moving ions due to a far smaller electron mass than ion mass. The radial distributions of the positive potential and electron density predicted from this study are compared with the available experimental data.

Two color interferometric electron density measurement in an axially blown arc

Science.gov (United States)

Stoller, Patrick; Carstensen, Jan; Galletti, Bernardo; Doiron, Charles; Sokolov, Alexey; Salzmann, René; Simon, Sandor; Jabs, Philipp

2016-09-01

High voltage circuit breakers protect the power grid by interrupting the current in case of a short circuit. To do so an arc is ignited between two contacts as they separate; transonic gas flow is used to cool and ultimately extinguish the arc at a current-zero crossing of the alternating current. A detailed understanding of the arc interruption process is needed to improve circuit breaker design. The conductivity of the partially ionized gas remaining after the current-zero crossing, a key parameter in determining whether the arc will be interrupted or not, is a function of the electron density. The electron density, in turn, is a function of the detailed dynamics of the arc cooling process, which does not necessarily occur under local thermodynamic equilibrium (LTE) conditions. In this work, we measure the spatially resolved line-integrated index of refraction in a near-current-zero arc stabilized in an axial flow of synthetic air with two nanosecond pulsed lasers at wavelengths of 532 nm and 671 nm. Generating a stable, cylindrically symmetric arc enables us to determine the three-dimensional index of refraction distribution using Abel inversion. Due to the wavelength dependence of the component of the index of refraction related to the free electrons, the information at two different wavelengths can be used to determine the electron density. This information allows us to determine how important it is to take into account non-equilibrium effects for accurate modeling of the physics of decaying arcs.

Density matrix embedding in an antisymmetrized geminal power bath

International Nuclear Information System (INIS)

Tsuchimochi, Takashi; Welborn, Matthew; Van Voorhis, Troy

2015-01-01

Density matrix embedding theory (DMET) has emerged as a powerful tool for performing wave function-in-wave function embedding for strongly correlated systems. In traditional DMET, an accurate calculation is performed on a small impurity embedded in a mean field bath. Here, we extend the original DMET equations to account for correlation in the bath via an antisymmetrized geminal power (AGP) wave function. The resulting formalism has a number of advantages. First, it allows one to properly treat the weak correlation limit of independent pairs, which DMET is unable to do with a mean-field bath. Second, it associates a size extensive correlation energy with a given density matrix (for the models tested), which AGP by itself is incapable of providing. Third, it provides a reasonable description of charge redistribution in strongly correlated but non-periodic systems. Thus, AGP-DMET appears to be a good starting point for describing electron correlation in molecules, which are aperiodic and possess both strong and weak electron correlation

On the dilute gas two particle density matrices of p--H2 and He4

International Nuclear Information System (INIS)

Weres, O.

1976-01-01

In the preceding paper we demonstrated that the reduced two- particle density matrix of simple quantum liquids could profitably be re-expressed in terms of a Taylor expansion of its logarithm about the diagonal. In the present publication we examine the Taylor coefficients which arise when the dilute gas two particle density matrix is expanded in this way. In particular, we evaluate the leading coefficients of p-H 2 and He 4 exactly and extend the Wigner--Kirkwood approximation to provided approximate expressions for them. We demonstrate how these approximate expressions may be applied to yield results superior to those yielded by the ordinary Wigner--Kirkwood approximation. In an appendix we demonstrate how the Block equation for the dilute gas two particle density matrix may be reduced to an equivalent closed set of equations for the leading Taylor coefficients

Bond index: relation to second-order density matrix and charge fluctuations

International Nuclear Information System (INIS)

Giambiagi, M.S. de; Giambiagi, M.; Jorge, F.E.

1985-01-01

It is shown that, in the same way as the atomic charge is an invariant built from the first-order density matrix, the closed-shell generalized bond index is an invariant associated with the second-order reduced density matrix. The active charge of an atom (sum of bond indices) is shown to be the sum of all density correlation functions between it and the other atoms in the molecule; similarly, the self-charge is the fluctuation of its total charge. (Author) [pt

'Anomalous electron transport' with 'Giant Current Density' at room temperature observed with nanogranular materials

International Nuclear Information System (INIS)

Koops, Hans W.P.

2013-01-01

Focused electron beam induced deposition is a novel bottom up nano-structurization technology. An electron beam of high power density is used to generate nano- structures with dimensions > 20 nm, but being composed from amorphous or nanogranular materials with crystals of 2 to 5 nm diameter embedded in a Fullerene matrix. Those compounds are generated in general by secondary or low energy electrons in layers of inorganic, organic, organometallic compounds absorbed to the sample. Those are converted into nanogranular materials by the electron beam following chemical and physical laws, as given by 'Mother Nature'. Metals and amorphous mixtures of chemical compounds from metals are normal resistors, which can carry a current density J 2 . Nanogranular composites like Au/C or Pt/C with metal nanocrystals embedded in a Fullerene matrix have hopping conduction with 0-dimensional Eigen-value characteristics and show 'anomalous electron transport' and can carry 'Giant Current Densities' with values from > 1 MA/cm 2 to 0.1 GA/cm 2 without destruction of the materials. However the area connecting the nanogranular material with a metal with a 3-dimensional electron gas needs to be designed, that the flowing current is reduced to the current density values which the 3-D metal can support without segregation. The basis for a theoretical explanation of the phenomenon can be geometry quantization for Coulomb blockade, of electron surface orbitals around the nanocrystals, hopping conduction, and the limitation of the density of states for phonons in geometry confined non percolated granular materials with strong difference in mass and orientation. Several applications in electronics, signal generators, light sources, detectors, and solar energy harvesting are suggested. (author)

Method of computer algebraic calculation of the matrix elements in the second quantization language

International Nuclear Information System (INIS)

Gotoh, Masashi; Mori, Kazuhide; Itoh, Reikichi

1995-01-01

An automated method by the algebraic programming language REDUCE3 for specifying the matrix elements expressed in second quantization language is presented and then applied to the case of the matrix elements in the TDHF theory. This program works in a very straightforward way by commuting the electron creation and annihilation operator (a † and a) until these operators have completely vanished from the expression of the matrix element under the appropriate elimination conditions. An improved method using singlet generators of unitary transformations in the place of the electron creation and annihilation operators is also presented. This improvement reduces the time and memory required for the calculation. These methods will make programming in the field of quantum chemistry much easier. 11 refs., 1 tab

Force-balance and differential equation for the ground-state electron density in atoms and molecules

International Nuclear Information System (INIS)

Amovilli, C.; March, N.H.; Gal, T.; Nagy, A.

2000-01-01

Holas and March (1995) established a force-balance equation from the many-electron Schroedinger equation. Here, the authors propose this as a basis for the construction of a (usually approximate) differential equation for the ground-state electron density. By way of example they present the simple case of two-electron systems with different external potentials but with weak electron-electron Coulomb repulsion λe 2 /r 12 . In this case first-order Rayleigh-Schroedinger (RS) perturbation theory of the ground-state wave function is known to lead to a compact expression for the first-order density matrix γ(r,rprime) in terms of its diagonal density ρ(r) and the density corresponding to λ = 0. This result allows the force-balance equation to be written as a third-order linear, differential homogeneous equation for the ground-state electron density ρ(r). The example of the two-electron Hookean atom is treated: For this case one can also transcend the first-order RS perturbation theory and get exact results for discrete choices of force constants (external potential)

C library for topological study of the electronic charge density.

Science.gov (United States)

Vega, David; Aray, Yosslen; Rodríguez, Jesús

2012-12-05

The topological study of the electronic charge density is useful to obtain information about the kinds of bonds (ionic or covalent) and the atom charges on a molecule or crystal. For this study, it is necessary to calculate, at every space point, the electronic density and its electronic density derivatives values up to second order. In this work, a grid-based method for these calculations is described. The library, implemented for three dimensions, is based on a multidimensional Lagrange interpolation in a regular grid; by differentiating the resulting polynomial, the gradient vector, the Hessian matrix and the Laplacian formulas were obtained for every space point. More complex functions such as the Newton-Raphson method (to find the critical points, where the gradient is null) and the Cash-Karp Runge-Kutta method (used to make the gradient paths) were programmed. As in some crystals, the unit cell has angles different from 90°, the described library includes linear transformations to correct the gradient and Hessian when the grid is distorted (inclined). Functions were also developed to handle grid containing files (grd from DMol® program, CUBE from Gaussian® program and CHGCAR from VASP® program). Each one of these files contains the data for a molecular or crystal electronic property (such as charge density, spin density, electrostatic potential, and others) in a three-dimensional (3D) grid. The library can be adapted to make the topological study in any regular 3D grid by modifying the code of these functions. Copyright © 2012 Wiley Periodicals, Inc.

The rigorous stochastic matrix multiplication scheme for the calculations of reduced equilibrium density matrices of open multilevel quantum systems

International Nuclear Information System (INIS)

Chen, Xin

2014-01-01

Understanding the roles of the temporary and spatial structures of quantum functional noise in open multilevel quantum molecular systems attracts a lot of theoretical interests. I want to establish a rigorous and general framework for functional quantum noises from the constructive and computational perspectives, i.e., how to generate the random trajectories to reproduce the kernel and path ordering of the influence functional with effective Monte Carlo methods for arbitrary spectral densities. This construction approach aims to unify the existing stochastic models to rigorously describe the temporary and spatial structure of Gaussian quantum noises. In this paper, I review the Euclidean imaginary time influence functional and propose the stochastic matrix multiplication scheme to calculate reduced equilibrium density matrices (REDM). In addition, I review and discuss the Feynman-Vernon influence functional according to the Gaussian quadratic integral, particularly its imaginary part which is critical to the rigorous description of the quantum detailed balance. As a result, I establish the conditions under which the influence functional can be interpreted as the average of exponential functional operator over real-valued Gaussian processes for open multilevel quantum systems. I also show the difference between the local and nonlocal phonons within this framework. With the stochastic matrix multiplication scheme, I compare the normalized REDM with the Boltzmann equilibrium distribution for open multilevel quantum systems

A simplified density matrix minimization for linear scaling self-consistent field theory

International Nuclear Information System (INIS)

Challacombe, M.

1999-01-01

A simplified version of the Li, Nunes and Vanderbilt [Phys. Rev. B 47, 10891 (1993)] and Daw [Phys. Rev. B 47, 10895 (1993)] density matrix minimization is introduced that requires four fewer matrix multiplies per minimization step relative to previous formulations. The simplified method also exhibits superior convergence properties, such that the bulk of the work may be shifted to the quadratically convergent McWeeny purification, which brings the density matrix to idempotency. Both orthogonal and nonorthogonal versions are derived. The AINV algorithm of Benzi, Meyer, and Tuma [SIAM J. Sci. Comp. 17, 1135 (1996)] is introduced to linear scaling electronic structure theory, and found to be essential in transformations between orthogonal and nonorthogonal representations. These methods have been developed with an atom-blocked sparse matrix algebra that achieves sustained megafloating point operations per second rates as high as 50% of theoretical, and implemented in the MondoSCF suite of linear scaling SCF programs. For the first time, linear scaling Hartree - Fock theory is demonstrated with three-dimensional systems, including water clusters and estane polymers. The nonorthogonal minimization is shown to be uncompetitive with minimization in an orthonormal representation. An early onset of linear scaling is found for both minimal and double zeta basis sets, and crossovers with a highly optimized eigensolver are achieved. Calculations with up to 6000 basis functions are reported. The scaling of errors with system size is investigated for various levels of approximation. copyright 1999 American Institute of Physics

Two-Dimensional Electron Density Measurement of Positive Streamer Discharge in Atmospheric-Pressure Air

Science.gov (United States)

Inada, Yuki; Ono, Ryo; Kumada, Akiko; Hidaka, Kunihiko; Maeyama, Mitsuaki

2016-09-01

The electron density of streamer discharges propagating in atmospheric-pressure air is crucially important for systematic understanding of the production mechanisms of reactive species utilized in wide ranging applications such as medical treatment, plasma-assisted ignition and combustion, ozone production and environmental pollutant processing. However, electron density measurement during the propagation of the atmospheric-pressure streamers is extremely difficult by using the conventional localized type measurement systems due to the streamer initiation jitters and the irreproducibility in the discharge paths. In order to overcome the difficulties, single-shot two-dimensional electron density measurement was conducted by using a Shack-Hartmann type laser wavefront sensor. The Shack-Hartmann sensor with a temporal resolution of 2 ns was applied to pulsed positive streamer discharges generated in an air gap between pin-to-plate electrodes. The electron density a few ns after the streamer initiation was 7*1021m-3 and uniformly distributed along the streamer channel. The electron density and its distribution profile were compared with a previous study simulating similar streamers, demonstrating good agreement. This work was supported in part by JKA and its promotion funds from KEIRIN RACE. The authors like to thank Mr. Kazuaki Ogura and Mr. Kaiho Aono of The University of Tokyo for their support during this work.

Comparison of three methods to reduce energy density. Effects on daily energy intake.

Science.gov (United States)

Williams, Rachel A; Roe, Liane S; Rolls, Barbara J

2013-07-01

Reductions in food energy density can decrease energy intake, but it is not known if the effects depend on the way that energy density is reduced. We investigated whether three methods of reducing energy density (decreasing fat, increasing fruit and vegetables, and adding water) differed in their effects on energy intake across the day. In a crossover design, 59 adults ate breakfast, lunch, and dinner in the laboratory once a week for 4 weeks. Across conditions, the entrées were either standard in energy density or were reduced in energy density by 20% using one of the three methods. Each meal included a manipulated entrée along with unmanipulated side dishes, and all foods were consumed ad libitum. Reducing the energy density of entrées significantly decreased daily energy intake compared to standard entrées (mean intake 2667 ± 77 kcal/day; 11,166 ± 322 kJ/day). The mean decrease was 396 ± 44 kcal/day (1658 ± 184 kJ/day) when fat was reduced, 308 ± 41 kcal/day (1290 ± 172 kJ/day) when fruit and vegetables were increased, and 230 ± 35 kcal/day (963 ± 147 kJ/day) when water was added. Daily energy intake was lower when fat was decreased compared to the other methods. These findings indicate that a variety of diet compositions can be recommended to reduce overall dietary energy density in order to moderate energy intake. Copyright © 2013 Elsevier Ltd. All rights reserved.

Exchange-correlation energies of atoms from efficient density functionals: influence of the electron density

Science.gov (United States)

Tao, Jianmin; Ye, Lin-Hui; Duan, Yuhua

2017-12-01

The primary goal of Kohn-Sham density functional theory is to evaluate the exchange-correlation contribution to electronic properties. However, the accuracy of a density functional can be affected by the electron density. Here we apply the nonempirical Tao-Mo (TM) semilocal functional to study the influence of the electron density on the exchange and correlation energies of atoms and ions, and compare the results with the commonly used nonempirical semilocal functionals local spin-density approximation (LSDA), Perdew-Burke-Ernzerhof (PBE), Tao-Perdew-Staroverov-Scuseria (TPSS), and hybrid functional PBE0. We find that the spin-restricted Hartree-Fock density yields the exchange and correlation energies in good agreement with the Optimized Effective Potential method, particularly for spherical atoms and ions. However, the errors of these semilocal and hybrid functionals become larger for self-consistent densities. We further find that the quality of the electron density have greater effect on the exchange-correlation energies of kinetic energy density-dependent meta-GGA functionals TPSS and TM than on those of the LSDA and GGA, and therefore, should have greater influence on the performance of meta-GGA functionals. Finally, we show that the influence of the density quality on PBE0 is slightly reduced, compared to that of PBE, due to the exact mixing.

Comment on "Nonuniqueness of algebraic first-order density-matrix functionals"

Science.gov (United States)

Gritsenko, O. V.

2018-02-01

Wang and Knowles (WK) [Phys. Rev. A 92, 012520 (2015), 10.1103/PhysRevA.92.012520] have given a counterexample to the conventional in reduced density-matrix functional theory representation of the second-order reduced density matrix (2RDM) Γi j ,k l in the basis of the natural orbitals as a function Γi j ,k l(n ) of the orbital occupation numbers (ONs) ni. The observed nonuniqueness of Γi j ,k l for prototype systems of different symmetry has been interpreted as the inherent inability of ON functions to reproduce the 2RDM, due to the insufficient information contained in the 1RDM spectrum. In this Comment, it is argued that, rather than totally invalidating Γi j ,k l(n ) , the WK example exposes its symmetry dependence which, as well as the previously established analogous dependence in density functional theory, is demonstrated with a general formulation based on the Levy constrained search.

Influence of electron correlation and degeneracy on the Fukui matrix and extension of frontier molecular orbital theory to correlated quantum chemical methods.

Science.gov (United States)

Bultinck, Patrick; Van Neck, Dimitri; Acke, Guillaume; Ayers, Paul W

2012-02-21

The Fukui function is considered as the diagonal element of the Fukui matrix in position space, where the Fukui matrix is the derivative of the one particle density matrix (1DM) with respect to the number of electrons. Diagonalization of the Fukui matrix, expressed in an orthogonal orbital basis, explains why regions in space with negative Fukui functions exist. Using a test set of molecules, electron correlation is found to have a remarkable effect on the eigenvalues of the Fukui matrix. The Fukui matrices at the independent electron model level are mathematically proven to always have an eigenvalue equal to exactly unity while the rest of the eigenvalues possibly differ from zero but sum to zero. The loss of idempotency of the 1DM at correlated levels of theory causes the loss of these properties. The influence of electron correlation is examined in detail and the frontier molecular orbital concept is extended to correlated levels of theory by defining it as the eigenvector of the Fukui matrix with the largest eigenvalue. The effect of degeneracy on the Fukui matrix is examined in detail, revealing that this is another way by which the unity eigenvalue and perfect pairing of eigenvalues can disappear.

Exact correlated kinetic energy related to the electron density for two-electron model atoms with harmonic confinement

International Nuclear Information System (INIS)

March, Norman H.; Akbari, Ali; Rubio, Angel

2007-01-01

For arbitrary interparticle interaction u(r 12 ), the model two-electron atom in the title is shown to be such that the ground-state electron density ρ(r) is determined uniquely by the correlated kinetic energy density t R (r) of the relative motion. Explicit results for t R (r) are presented for the Hookean atom with force constant k=1/4, and also for u(r 12 )=(λ)/(r 12 2 ) . Possible relevance of the Hookean atom treatment to the ground state of the helium atom itself is briefly discussed

DensToolKit: A comprehensive open-source package for analyzing the electron density and its derivative scalar and vector fields

Science.gov (United States)

Solano-Altamirano, J. M.; Hernández-Pérez, Julio M.

2015-11-01

DensToolKit is a suite of cross-platform, optionally parallelized, programs for analyzing the molecular electron density (ρ) and several fields derived from it. Scalar and vector fields, such as the gradient of the electron density (∇ρ), electron localization function (ELF) and its gradient, localized orbital locator (LOL), region of slow electrons (RoSE), reduced density gradient, localized electrons detector (LED), information entropy, molecular electrostatic potential, kinetic energy densities K and G, among others, can be evaluated on zero, one, two, and three dimensional grids. The suite includes a program for searching critical points and bond paths of the electron density, under the framework of Quantum Theory of Atoms in Molecules. DensToolKit also evaluates the momentum space electron density on spatial grids, and the reduced density matrix of order one along lines joining two arbitrary atoms of a molecule. The source code is distributed under the GNU-GPLv3 license, and we release the code with the intent of establishing an open-source collaborative project. The style of DensToolKit's code follows some of the guidelines of an object-oriented program. This allows us to supply the user with a simple manner for easily implement new scalar or vector fields, provided they are derived from any of the fields already implemented in the code. In this paper, we present some of the most salient features of the programs contained in the suite, some examples of how to run them, and the mathematical definitions of the implemented fields along with hints of how we optimized their evaluation. We benchmarked our suite against both a freely-available program and a commercial package. Speed-ups of ˜2×, and up to 12× were obtained using a non-parallel compilation of DensToolKit for the evaluation of fields. DensToolKit takes similar times for finding critical points, compared to a commercial package. Finally, we present some perspectives for the future development and

Teaching Chemistry with Electron Density Models

Science.gov (United States)

Shusterman, Gwendolyn P.; Shusterman, Alan J.

1997-07-01

Linus Pauling once said that a topic must satisfy two criteria before it can be taught to students. First, students must be able to assimilate the topic within a reasonable amount of time. Second, the topic must be relevant to the educational needs and interests of the students. Unfortunately, the standard general chemistry textbook presentation of "electronic structure theory", set as it is in the language of molecular orbitals, has a difficult time satisfying either criterion. Many of the quantum mechanical aspects of molecular orbitals are too difficult for most beginning students to appreciate, much less master, and the few applications that are presented in the typical textbook are too limited in scope to excite much student interest. This article describes a powerful new method for teaching students about electronic structure and its relevance to chemical phenomena. This method, which we have developed and used for several years in general chemistry (G.P.S.) and organic chemistry (A.J.S.) courses, relies on computer-generated three-dimensional models of electron density distributions, and largely satisfies Pauling's two criteria. Students find electron density models easy to understand and use, and because these models are easily applied to a broad range of topics, they successfully convey to students the importance of electronic structure. In addition, when students finally learn about orbital concepts they are better prepared because they already have a well-developed three-dimensional picture of electronic structure to fall back on. We note in this regard that the types of models we use have found widespread, rigorous application in chemical research (1, 2), so students who understand and use electron density models do not need to "unlearn" anything before progressing to more advanced theories.

Geometrical separation method for lipoproteins using bioformulated-fiber matrix electrophoresis: size of high-density lipoprotein does not reflect its density.

Science.gov (United States)

Tabuchi, Mari; Seo, Makoto; Inoue, Takayuki; Ikeda, Takeshi; Kogure, Akinori; Inoue, Ikuo; Katayama, Shigehiro; Matsunaga, Toshiyuki; Hara, Akira; Komoda, Tsugikazu

2011-02-01

The increasing number of patients with metabolic syndrome is a critical global problem. In this study, we describe a novel geometrical electrophoretic separation method using a bioformulated-fiber matrix to analyze high-density lipoprotein (HDL) particles. HDL particles are generally considered to be a beneficial component of the cholesterol fraction. Conventional electrophoresis is widely used but is not necessarily suitable for analyzing HDL particles. Furthermore, a higher HDL density is generally believed to correlate with a smaller particle size. Here, we use a novel geometrical separation technique incorporating recently developed nanotechnology (Nata de Coco) to contradict this belief. A dyslipidemia patient given a 1-month treatment of fenofibrate showed an inverse relationship between HDL density and size. Direct microscopic observation and morphological observation of fractionated HDL particles confirmed a lack of relationship between particle density and size. This new technique may improve diagnostic accuracy and medical treatment for lipid related diseases.

Derivation of general analytic gradient expressions for density-fitted post-Hartree-Fock methods: An efficient implementation for the density-fitted second-order Møller–Plesset perturbation theory

Energy Technology Data Exchange (ETDEWEB)

Bozkaya, Uğur, E-mail: ugur.bozkaya@atauni.edu.tr [Department of Chemistry, Atatürk University, Erzurum 25240, Turkey and Department of Chemistry, Middle East Technical University, Ankara 06800 (Turkey)

2014-09-28

General analytic gradient expressions (with the frozen-core approximation) are presented for density-fitted post-HF methods. An efficient implementation of frozen-core analytic gradients for the second-order Møller–Plesset perturbation theory (MP2) with the density-fitting (DF) approximation (applying to both reference and correlation energies), which is denoted as DF-MP2, is reported. The DF-MP2 method is applied to a set of alkanes, conjugated dienes, and noncovalent interaction complexes to compare the computational cost of single point analytic gradients with MP2 with the resolution of the identity approach (RI-MP2) [F. Weigend and M. Häser, Theor. Chem. Acc. 97, 331 (1997); R. A. Distasio, R. P. Steele, Y. M. Rhee, Y. Shao, and M. Head-Gordon, J. Comput. Chem. 28, 839 (2007)]. In the RI-MP2 method, the DF approach is used only for the correlation energy. Our results demonstrate that the DF-MP2 method substantially accelerate the RI-MP2 method for analytic gradient computations due to the reduced input/output (I/O) time. Because in the DF-MP2 method the DF approach is used for both reference and correlation energies, the storage of 4-index electron repulsion integrals (ERIs) are avoided, 3-index ERI tensors are employed instead. Further, as in case of integrals, our gradient equation is completely avoid construction or storage of the 4-index two-particle density matrix (TPDM), instead we use 2- and 3-index TPDMs. Hence, the I/O bottleneck of a gradient computation is significantly overcome. Therefore, the cost of the generalized-Fock matrix (GFM), TPDM, solution of Z-vector equations, the back transformation of TPDM, and integral derivatives are substantially reduced when the DF approach is used for the entire energy expression. Further application results show that the DF approach introduce negligible errors for closed-shell reaction energies and equilibrium bond lengths.

Applications of density matrix in the fractional quantum mechanics: Thomas-Fermi model and Hohenberg-Kohn theorems revisited

International Nuclear Information System (INIS)

Dong, Jianping

2011-01-01

The many-body space fractional quantum system is studied using the density matrix method. We give the new results of the Thomas-Fermi model, obtain the quantum pressure of the free electron gas. We also show the validity of the Hohenberg-Kohn theorems in the space fractional quantum mechanics and generalize the density functional theory to the fractional quantum mechanics. -- Highlights: → Thomas-Fermi model under the framework of fractional quantum mechanics is studied. → We show the validity of the HK theorems in the space fractional quantum mechanics. → The density functional theory is generalized to the fractional quantum mechanics.

Unbiased reduced density matrices and electronic properties from full configuration interaction quantum Monte Carlo

International Nuclear Information System (INIS)

Overy, Catherine; Blunt, N. S.; Shepherd, James J.; Booth, George H.; Cleland, Deidre; Alavi, Ali

2014-01-01

Properties that are necessarily formulated within pure (symmetric) expectation values are difficult to calculate for projector quantum Monte Carlo approaches, but are critical in order to compute many of the important observable properties of electronic systems. Here, we investigate an approach for the sampling of unbiased reduced density matrices within the full configuration interaction quantum Monte Carlo dynamic, which requires only small computational overheads. This is achieved via an independent replica population of walkers in the dynamic, sampled alongside the original population. The resulting reduced density matrices are free from systematic error (beyond those present via constraints on the dynamic itself) and can be used to compute a variety of expectation values and properties, with rapid convergence to an exact limit. A quasi-variational energy estimate derived from these density matrices is proposed as an accurate alternative to the projected estimator for multiconfigurational wavefunctions, while its variational property could potentially lend itself to accurate extrapolation approaches in larger systems

A second-order unconstrained optimization method for canonical-ensemble density-functional methods

Science.gov (United States)

Nygaard, Cecilie R.; Olsen, Jeppe

2013-03-01

A second order converging method of ensemble optimization (SOEO) in the framework of Kohn-Sham Density-Functional Theory is presented, where the energy is minimized with respect to an ensemble density matrix. It is general in the sense that the number of fractionally occupied orbitals is not predefined, but rather it is optimized by the algorithm. SOEO is a second order Newton-Raphson method of optimization, where both the form of the orbitals and the occupation numbers are optimized simultaneously. To keep the occupation numbers between zero and two, a set of occupation angles is defined, from which the occupation numbers are expressed as trigonometric functions. The total number of electrons is controlled by a built-in second order restriction of the Newton-Raphson equations, which can be deactivated in the case of a grand-canonical ensemble (where the total number of electrons is allowed to change). To test the optimization method, dissociation curves for diatomic carbon are produced using different functionals for the exchange-correlation energy. These curves show that SOEO favors symmetry broken pure-state solutions when using functionals with exact exchange such as Hartree-Fock and Becke three-parameter Lee-Yang-Parr. This is explained by an unphysical contribution to the exact exchange energy from interactions between fractional occupations. For functionals without exact exchange, such as local density approximation or Becke Lee-Yang-Parr, ensemble solutions are favored at interatomic distances larger than the equilibrium distance. Calculations on the chromium dimer are also discussed. They show that SOEO is able to converge to ensemble solutions for systems that are more complicated than diatomic carbon.

Density-functional expansion methods: Grand challenges.

Science.gov (United States)

Giese, Timothy J; York, Darrin M

2012-03-01

We discuss the source of errors in semiempirical density functional expansion (VE) methods. In particular, we show that VE methods are capable of well-reproducing their standard Kohn-Sham density functional method counterparts, but suffer from large errors upon using one or more of these approximations: the limited size of the atomic orbital basis, the Slater monopole auxiliary basis description of the response density, and the one- and two-body treatment of the core-Hamiltonian matrix elements. In the process of discussing these approximations and highlighting their symptoms, we introduce a new model that supplements the second-order density-functional tight-binding model with a self-consistent charge-dependent chemical potential equalization correction; we review our recently reported method for generalizing the auxiliary basis description of the atomic orbital response density; and we decompose the first-order potential into a summation of additive atomic components and many-body corrections, and from this examination, we provide new insights and preliminary results that motivate and inspire new approximate treatments of the core-Hamiltonian.

Hybrid matrix method for stable numerical analysis of the propagation of Dirac electrons in gapless bilayer graphene superlattices

Science.gov (United States)

Briones-Torres, J. A.; Pernas-Salomón, R.; Pérez-Álvarez, R.; Rodríguez-Vargas, I.

2016-05-01

Gapless bilayer graphene (GBG), like monolayer graphene, is a material system with unique properties, such as anti-Klein tunneling and intrinsic Fano resonances. These properties rely on the gapless parabolic dispersion relation and the chiral nature of bilayer graphene electrons. In addition, propagating and evanescent electron states coexist inherently in this material, giving rise to these exotic properties. In this sense, bilayer graphene is unique, since in most material systems in which Fano resonance phenomena are manifested an external source that provides extended states is required. However, from a numerical standpoint, the presence of evanescent-divergent states in the eigenfunctions linear superposition representing the Dirac spinors, leads to a numerical degradation (the so called Ωd problem) in the practical applications of the standard Coefficient Transfer Matrix (K) method used to study charge transport properties in Bilayer Graphene based multi-barrier systems. We present here a straightforward procedure based in the hybrid compliance-stiffness matrix method (H) that can overcome this numerical degradation. Our results show that in contrast to standard matrix method, the proposed H method is suitable to study the transmission and transport properties of electrons in GBG superlattice since it remains numerically stable regardless the size of the superlattice and the range of values taken by the input parameters: the energy and angle of the incident electrons, the barrier height and the thickness and number of barriers. We show that the matrix determinant can be used as a test of the numerical accuracy in real calculations.

A unifying probabilistic Bayesian approach to derive electron density from MRI for radiation therapy treatment planning

International Nuclear Information System (INIS)

Gudur, Madhu Sudhan Reddy; Hara, Wendy; Le, Quynh-Thu; Wang, Lei; Xing, Lei; Li, Ruijiang

2014-01-01

MRI significantly improves the accuracy and reliability of target delineation in radiation therapy for certain tumors due to its superior soft tissue contrast compared to CT. A treatment planning process with MRI as the sole imaging modality will eliminate systematic CT/MRI co-registration errors, reduce cost and radiation exposure, and simplify clinical workflow. However, MRI lacks the key electron density information necessary for accurate dose calculation and generating reference images for patient setup. The purpose of this work is to develop a unifying method to derive electron density from standard T1-weighted MRI. We propose to combine both intensity and geometry information into a unifying probabilistic Bayesian framework for electron density mapping. For each voxel, we compute two conditional probability density functions (PDFs) of electron density given its: (1) T1-weighted MRI intensity, and (2) geometry in a reference anatomy, obtained by deformable image registration between the MRI of the atlas and test patient. The two conditional PDFs containing intensity and geometry information are combined into a unifying posterior PDF, whose mean value corresponds to the optimal electron density value under the mean-square error criterion. We evaluated the algorithm’s accuracy of electron density mapping and its ability to detect bone in the head for eight patients, using an additional patient as the atlas or template. Mean absolute HU error between the estimated and true CT, as well as receiver operating characteristics for bone detection (HU > 200) were calculated. The performance was compared with a global intensity approach based on T1 and no density correction (set whole head to water). The proposed technique significantly reduced the errors in electron density estimation, with a mean absolute HU error of 126, compared with 139 for deformable registration (p = 2  ×  10 −4 ), 283 for the intensity approach (p = 2  ×  10 −6 ) and 282

Significance of matrix diagonalization in modelling inelastic electron scattering

Energy Technology Data Exchange (ETDEWEB)

Lee, Z. [University of Ulm, Ulm 89081 (Germany); Hambach, R. [University of Ulm, Ulm 89081 (Germany); University of Jena, Jena 07743 (Germany); Kaiser, U.; Rose, H. [University of Ulm, Ulm 89081 (Germany)

2017-04-15

Electron scattering is always applied as one of the routines to investigate nanostructures. Nowadays the development of hardware offers more and more prospect for this technique. For example imaging nanostructures with inelastic scattered electrons may allow to produce component-sensitive images with atomic resolution. Modelling inelastic electron scattering is therefore essential for interpreting these images. The main obstacle to study inelastic scattering problem is its complexity. During inelastic scattering, incident electrons entangle with objects, and the description of this process involves a multidimensional array. Since the simulation usually involves fourdimensional Fourier transforms, the computation is highly inefficient. In this work we have offered one solution to handle the multidimensional problem. By transforming a high dimensional array into twodimensional array, we are able to perform matrix diagonalization and approximate the original multidimensional array with its twodimensional eigenvectors. Our procedure reduces the complicated multidimensional problem to a twodimensional problem. In addition, it minimizes the number of twodimensional problems. This method is very useful for studying multiple inelastic scattering. - Highlights: • 4D problems are involved in modelling inelastic electron scattering. • By means of matrix diagonalization, the 4D problems can be simplified as 2D problems. • The number of 2D problems is minimized by using this approach.

Electronic structure of the Fe2 molecule in the local-spin-density approximation

International Nuclear Information System (INIS)

Dhar, S.; Kestner, N.R.

1988-01-01

Ab initio self-consistent all-electron spin-polarized calculations have been performed for the ground-state properties of the Fe 2 molecule using the local-spin-density approximation. A Gaussian orbital basis is employed and all the two-electron integrals are evaluated analytically. The matrix elements of the exchange-correlation potential are computed numerically. The total energy, the binding energy, the equilibrium distance, vibrational frequency, and the ground-state configurations are reported and compared with other calculations and experimental results

Application of texture analysis method for mammogram density classification

Science.gov (United States)

Nithya, R.; Santhi, B.

2017-07-01

Mammographic density is considered a major risk factor for developing breast cancer. This paper proposes an automated approach to classify breast tissue types in digital mammogram. The main objective of the proposed Computer-Aided Diagnosis (CAD) system is to investigate various feature extraction methods and classifiers to improve the diagnostic accuracy in mammogram density classification. Texture analysis methods are used to extract the features from the mammogram. Texture features are extracted by using histogram, Gray Level Co-Occurrence Matrix (GLCM), Gray Level Run Length Matrix (GLRLM), Gray Level Difference Matrix (GLDM), Local Binary Pattern (LBP), Entropy, Discrete Wavelet Transform (DWT), Wavelet Packet Transform (WPT), Gabor transform and trace transform. These extracted features are selected using Analysis of Variance (ANOVA). The features selected by ANOVA are fed into the classifiers to characterize the mammogram into two-class (fatty/dense) and three-class (fatty/glandular/dense) breast density classification. This work has been carried out by using the mini-Mammographic Image Analysis Society (MIAS) database. Five classifiers are employed namely, Artificial Neural Network (ANN), Linear Discriminant Analysis (LDA), Naive Bayes (NB), K-Nearest Neighbor (KNN), and Support Vector Machine (SVM). Experimental results show that ANN provides better performance than LDA, NB, KNN and SVM classifiers. The proposed methodology has achieved 97.5% accuracy for three-class and 99.37% for two-class density classification.

Effect of disorder on the density of states of a two-dimensional electron gas under magnetic field

International Nuclear Information System (INIS)

Bonifacie, S.; Meziani, Y.M.; Chaubet, C.; Jouault, B.; Raymond, A.

2004-01-01

We have calculated the density of states (DOS) of a two-dimensional electron gas in a perpendicular magnetic field, using a multiple scattering method, in the ultraquantum limit. We have considered doped and disordered 2D systems. The results of the scattering method are compared with direct simulations of disordered samples. Using the DOS, we have studied the metal-insulator transition and the magnetic freeze-out including a comparison with experimental results

Diagnosis of Unmagnetized Plasma Electron Number Density and Electron-neutral Collision Frequency by Using Microwave

International Nuclear Information System (INIS)

Yuan Zhongcai; Shi Jiaming; Xu Bo

2005-01-01

The plasma diagnostic method using the transmission attenuation of microwaves at double frequencies (PDMUTAMDF) indicates that the frequency and the electron-neutral collision frequency of the plasma can be deduced by utilizing the transmission attenuation of microwaves at two neighboring frequencies in a non-magnetized plasma. Then the electron density can be obtained from the plasma frequency. The PDMUTAMDF is a simple method to diagnose the plasma indirectly. In this paper, the interaction of electromagnetic waves and the plasma is analyzed. Then, based on the attenuation and the phase shift of a microwave in the plasma, the principle of the PDMUTAMDF is presented. With the diagnostic method, the spatially mean electron density and electron collision frequency of the plasma can be obtained. This method is suitable for the elementary diagnosis of the atmospheric-pressure plasma

Two-electron Rabi oscillations in real-time time-dependent density-functional theory

International Nuclear Information System (INIS)

Habenicht, Bradley F.; Tani, Noriyuki P.; Provorse, Makenzie R.; Isborn, Christine M.

2014-01-01

We investigate the Rabi oscillations of electrons excited by an applied electric field in several simple molecular systems using time-dependent configuration interaction (TDCI) and real-time time-dependent density-functional theory (RT-TDDFT) dynamics. While the TDCI simulations exhibit the expected single-electron Rabi oscillations at a single resonant electric field frequency, Rabi oscillations in the RT-TDDFT simulations are a two-electron process. The existence of two-electron Rabi oscillations is determined both by full population inversion between field-free molecular orbitals and the behavior of the instantaneous dipole moment during the simulations. Furthermore, the Rabi oscillations in RT-TDDFT are subject to an intensity threshold of the electric field, below which Rabi oscillations do not occur and above which the two-electron Rabi oscillations occur at a broad range of frequencies. It is also shown that at field intensities near the threshold intensity, the field frequency predicted to induce Rabi oscillations by linear response TDDFT only produces detuned Rabi oscillations. Instead, the field frequency that yields the full two-electron population inversion and Rabi oscillation behavior is shown to be the average of single-electron transition frequencies from the ground S 0 state and the doubly-excited S 2 state. The behavior of the two-electron Rabi oscillations is rationalized via two possible models. The first model is a multi-photon process that results from the electric field interacting with the three level system such that three level Rabi oscillations may occur. The second model suggests that the mean-field nature of RT-TDDFT induces paired electron propagation

Comparison of three methods to reduce energy density: effects on daily energy intake

OpenAIRE

Williams, Rachel A.; Roe, Liane S.; Rolls, Barbara J.

2013-01-01

Reductions in food energy density can decrease energy intake, but it is not known if the effects depend on the way that energy density is reduced. We investigated whether three methods of reducing energy density (decreasing fat, increasing fruit and vegetables, and adding water) differed in their effects on energy intake across the day. In a crossover design, 59 adults ate breakfast, lunch, and dinner in the laboratory once a week for four weeks. Across conditions, the entrées were either sta...

High performance computing of density matrix renormalization group method for 2-dimensional model. Parallelization strategy toward peta computing

International Nuclear Information System (INIS)

Yamada, Susumu; Igarashi, Ryo; Machida, Masahiko; Imamura, Toshiyuki; Okumura, Masahiko; Onishi, Hiroaki

2010-01-01

We parallelize the density matrix renormalization group (DMRG) method, which is a ground-state solver for one-dimensional quantum lattice systems. The parallelization allows us to extend the applicable range of the DMRG to n-leg ladders i.e., quasi two-dimension cases. Such an extension is regarded to bring about several breakthroughs in e.g., quantum-physics, chemistry, and nano-engineering. However, the straightforward parallelization requires all-to-all communications between all processes which are unsuitable for multi-core systems, which is a mainstream of current parallel computers. Therefore, we optimize the all-to-all communications by the following two steps. The first one is the elimination of the communications between all processes by only rearranging data distribution with the communication data amount kept. The second one is the avoidance of the communication conflict by rescheduling the calculation and the communication. We evaluate the performance of the DMRG method on multi-core supercomputers and confirm that our two-steps tuning is quite effective. (author)

Stabilization of electron-scale turbulence by electron density gradient in national spherical torus experiment

Energy Technology Data Exchange (ETDEWEB)

Ruiz Ruiz, J.; White, A. E. [MIT-Plasma Science and Fusion Center, Cambridge, Massachusetts 02139 (United States); Ren, Y.; Guttenfelder, W.; Kaye, S. M.; Leblanc, B. P.; Mazzucato, E. [Princeton Plasma Physics Laboratory, Princeton, New Jersey 08543 (United States); Lee, K. C. [National Fusion Research Institute, Daejeon (Korea, Republic of); Domier, C. W. [University of California at Davis, Davis, California 95616 (United States); Smith, D. R. [University of Wisconsin-Madison, Madison, Wisconsin 53706 (United States); Yuh, H. [Nova Photonics, Inc., Princeton, New Jersey 08540 (United States)

2015-12-15

Theory and experiments have shown that electron temperature gradient (ETG) turbulence on the electron gyro-scale, k{sub ⊥}ρ{sub e} ≲ 1, can be responsible for anomalous electron thermal transport in NSTX. Electron scale (high-k) turbulence is diagnosed in NSTX with a high-k microwave scattering system [D. R. Smith et al., Rev. Sci. Instrum. 79, 123501 (2008)]. Here we report on stabilization effects of the electron density gradient on electron-scale density fluctuations in a set of neutral beam injection heated H-mode plasmas. We found that the absence of high-k density fluctuations from measurements is correlated with large equilibrium density gradient, which is shown to be consistent with linear stabilization of ETG modes due to the density gradient using the analytical ETG linear threshold in F. Jenko et al. [Phys. Plasmas 8, 4096 (2001)] and linear gyrokinetic simulations with GS2 [M. Kotschenreuther et al., Comput. Phys. Commun. 88, 128 (1995)]. We also found that the observed power of electron-scale turbulence (when it exists) is anti-correlated with the equilibrium density gradient, suggesting density gradient as a nonlinear stabilizing mechanism. Higher density gradients give rise to lower values of the plasma frame frequency, calculated based on the Doppler shift of the measured density fluctuations. Linear gyrokinetic simulations show that higher values of the electron density gradient reduce the value of the real frequency, in agreement with experimental observation. Nonlinear electron-scale gyrokinetic simulations show that high electron density gradient reduces electron heat flux and stiffness, and increases the ETG nonlinear threshold, consistent with experimental observations.

Single-particle density matrix of liquid 4He

International Nuclear Information System (INIS)

Vakarchuk, I.A.

2008-01-01

The density single-particle matrix in the coordinate notation was calculated based on the expression for the interacting Bose-particle N system density matrix. Under the low temperatures the mentioned matrix in the first approximation enables to reproduce the Bogoliubov theory results. In the classical terms the mentioned theory enables to reproduce the results of the theory of the classical fluids in the approximation of the chaotic phases. On the basis of the density single-particle matrix one managed to obtain the function of the pulse distribution of the particles, the Bose-liquid average kinetic energy, and to study the Bose-Einstein condensation phenomenon [ru

SU-G-JeP2-02: A Unifying Multi-Atlas Approach to Electron Density Mapping Using Multi-Parametric MRI for Radiation Treatment Planning

Energy Technology Data Exchange (ETDEWEB)

Ren, S [Stanford University, Stanford, CA (United States); Tianjin University, Tianjin (China); Hara, W; Le, Q; Wang, L; Xing, L; Li, R [Stanford University, Stanford, CA (United States)

2016-06-15

Purpose: MRI has a number of advantages over CT as a primary modality for radiation treatment planning (RTP). However, one key bottleneck problem still remains, which is the lack of electron density information in MRI. In the work, a reliable method to map electron density is developed by leveraging the differential contrast of multi-parametric MRI. Methods: We propose a probabilistic Bayesian approach for electron density mapping based on T1 and T2-weighted MRI, using multiple patients as atlases. For each voxel, we compute two conditional probabilities: (1) electron density given its image intensity on T1 and T2-weighted MR images, and (2) electron density given its geometric location in a reference anatomy. The two sources of information (image intensity and spatial location) are combined into a unifying posterior probability density function using the Bayesian formalism. The mean value of the posterior probability density function provides the estimated electron density. Results: We evaluated the method on 10 head and neck patients and performed leave-one-out cross validation (9 patients as atlases and remaining 1 as test). The proposed method significantly reduced the errors in electron density estimation, with a mean absolute HU error of 138, compared with 193 for the T1-weighted intensity approach and 261 without density correction. For bone detection (HU>200), the proposed method had an accuracy of 84% and a sensitivity of 73% at specificity of 90% (AUC = 87%). In comparison, the AUC for bone detection is 73% and 50% using the intensity approach and without density correction, respectively. Conclusion: The proposed unifying method provides accurate electron density estimation and bone detection based on multi-parametric MRI of the head with highly heterogeneous anatomy. This could allow for accurate dose calculation and reference image generation for patient setup in MRI-based radiation treatment planning.

SU-G-JeP2-02: A Unifying Multi-Atlas Approach to Electron Density Mapping Using Multi-Parametric MRI for Radiation Treatment Planning

International Nuclear Information System (INIS)

Ren, S; Hara, W; Le, Q; Wang, L; Xing, L; Li, R

2016-01-01

Purpose: MRI has a number of advantages over CT as a primary modality for radiation treatment planning (RTP). However, one key bottleneck problem still remains, which is the lack of electron density information in MRI. In the work, a reliable method to map electron density is developed by leveraging the differential contrast of multi-parametric MRI. Methods: We propose a probabilistic Bayesian approach for electron density mapping based on T1 and T2-weighted MRI, using multiple patients as atlases. For each voxel, we compute two conditional probabilities: (1) electron density given its image intensity on T1 and T2-weighted MR images, and (2) electron density given its geometric location in a reference anatomy. The two sources of information (image intensity and spatial location) are combined into a unifying posterior probability density function using the Bayesian formalism. The mean value of the posterior probability density function provides the estimated electron density. Results: We evaluated the method on 10 head and neck patients and performed leave-one-out cross validation (9 patients as atlases and remaining 1 as test). The proposed method significantly reduced the errors in electron density estimation, with a mean absolute HU error of 138, compared with 193 for the T1-weighted intensity approach and 261 without density correction. For bone detection (HU>200), the proposed method had an accuracy of 84% and a sensitivity of 73% at specificity of 90% (AUC = 87%). In comparison, the AUC for bone detection is 73% and 50% using the intensity approach and without density correction, respectively. Conclusion: The proposed unifying method provides accurate electron density estimation and bone detection based on multi-parametric MRI of the head with highly heterogeneous anatomy. This could allow for accurate dose calculation and reference image generation for patient setup in MRI-based radiation treatment planning.

A unit density method of grain analysis used to identify GABEergic neurons for electron microscopic autoradiographs

International Nuclear Information System (INIS)

Burry, R.W.

1982-01-01

The distribution of electron microscopic autoradiographic grains over neurons in cerebellar cultures incubated with [ 3 H]gamma-aminobutyric acid ([ 3 H]GABA) was examined. With the unit density method of grain analysis, the number of grains over each structure was tested against the total grain density for the entire section. If an individual structure has a grain density higher than the expected grain density, it is considered one of the group of heavily labeled structures. The expected grain density for each structure is calculated based on the area for that structure, the total grain density and the Poisson distribution. A different expected grain density can be calculated for any P value required. The method provides an adequate population of structures for morphological analysis but excludes weakly labeled structures and thus may underestimate the number of labeled structures. The unit density method of grain analysis showed, as expected, a group of cell bodies and synapses that was labeled heavily. Cultures incubated with other [ 3 H]amino acids did not have any heavily labeled synaptic elements. In addition, serial section analysis of sections showed that synapses heavily labeled with [ 3 H]GABA are seen in adjacent sections. The advantage of the unit density method of grain analysis is that it can be used to separate two groups of metabolically different neurons even when no morphological differences are present. (Auth.)

General filtering method for electronic speckle pattern interferometry fringe images with various densities based on variational image decomposition.

Science.gov (United States)

Li, Biyuan; Tang, Chen; Gao, Guannan; Chen, Mingming; Tang, Shuwei; Lei, Zhenkun

2017-06-01

Filtering off speckle noise from a fringe image is one of the key tasks in electronic speckle pattern interferometry (ESPI). In general, ESPI fringe images can be divided into three categories: low-density fringe images, high-density fringe images, and variable-density fringe images. In this paper, we first present a general filtering method based on variational image decomposition that can filter speckle noise for ESPI fringe images with various densities. In our method, a variable-density ESPI fringe image is decomposed into low-density fringes, high-density fringes, and noise. A low-density fringe image is decomposed into low-density fringes and noise. A high-density fringe image is decomposed into high-density fringes and noise. We give some suitable function spaces to describe low-density fringes, high-density fringes, and noise, respectively. Then we construct several models and numerical algorithms for ESPI fringe images with various densities. And we investigate the performance of these models via our extensive experiments. Finally, we compare our proposed models with the windowed Fourier transform method and coherence enhancing diffusion partial differential equation filter. These two methods may be the most effective filtering methods at present. Furthermore, we use the proposed method to filter a collection of the experimentally obtained ESPI fringe images with poor quality. The experimental results demonstrate the performance of our proposed method.

Performance Tuning of Fock Matrix and Two-Electron Integral Calculations for NWChem on Leading HPC Platforms

Energy Technology Data Exchange (ETDEWEB)

Shan, Hongzhan; Austin, Brian M.; De Jong, Wibe A.; Oliker, Leonid; Wright, Nicholas J.; Apra, Edoardo

2014-10-01

Attaining performance in the evaluation of two-electron repulsion integrals and constructing the Fock matrix is of considerable importance to the computational chemistry community. Due to its numerical complexity improving the performance behavior across a variety of leading supercomputing platforms is an increasing challenge due to the significant diversity in high-performance computing architectures. In this paper, we present our successful tuning methodology for these important numerical methods on the Cray XE6, the Cray XC30, the IBM BG/Q, as well as the Intel Xeon Phi. Our optimization schemes leverage key architectural features including vectorization and simultaneous multithreading, and results in speedups of up to 2.5x compared with the original implementation.

ABCD Matrix Method a Case Study

CERN Document Server

Seidov, Zakir F; Yahalom, Asher

2004-01-01

In the Israeli Electrostatic Accelerator FEL, the distance between the accelerator's end and the wiggler's entrance is about 2.1 m, and 1.4 MeV electron beam is transported through this space using four similar quadrupoles (FODO-channel). The transfer matrix method (ABCD matrix method) was used for simulating the beam transport, a set of programs is written in the several programming languages (MATHEMATICA, MATLAB, MATCAD, MAPLE) and reasonable agreement is demonstrated between experimental results and simulations. Comparison of ABCD matrix method with the direct "numerical experiments" using EGUN, ELOP, and GPT programs with and without taking into account the space-charge effects showed the agreement to be good enough as well. Also the inverse problem of finding emittance of the electron beam at the S1 screen position (before FODO-channel), by using the spot image at S2 screen position (after FODO-channel) as function of quad currents, is considered. Spot and beam at both screens are described as tilted eel...

Convergent j-matrix calculation of electron-helium resonances

International Nuclear Information System (INIS)

Konovalov, D.A.; McCarthy, I.E.

1994-12-01

Resonance structures in n=2 and n=3 electron-helium excitation cross sections are calculated using the J-matrix method. The number of close-coupled helium bound and continuum states is taken to convergence, e.g. about 100 channels are coupled for each total spin and angular momentum. It is found that the present J-matrix results are in good shape agreement with recent 29-state R-matrix calculations. However the J-matrix absolute cross sections are slightly lower due to the influence of continuum channels included in the present method. Experiment and theory agree on the positions of n=2 and n=3 resonances. 22 refs., 1 tab.; 3 figs

A novel matrix approach for controlling the invariant densities of chaotic maps

International Nuclear Information System (INIS)

Rogers, Alan; Shorten, Robert; Heffernan, Daniel M.

2008-01-01

Recent work on positive matrices has resulted in a new matrix method for generating chaotic maps with arbitrary piecewise constant invariant densities, sometimes known as the inverse Frobenius-Perron problem (IFPP). In this paper, we give an extensive introduction to the IFPP, describing existing methods for solving it, and we describe our new matrix approach for solving the IFPP

A Comparative Study of Collagen Matrix Density Effect on Endothelial Sprout Formation Using Experimental and Computational Approaches.

Science.gov (United States)

Shamloo, Amir; Mohammadaliha, Negar; Heilshorn, Sarah C; Bauer, Amy L

2016-04-01

A thorough understanding of determining factors in angiogenesis is a necessary step to control the development of new blood vessels. Extracellular matrix density is known to have a significant influence on cellular behaviors and consequently can regulate vessel formation. The utilization of experimental platforms in combination with numerical models can be a powerful method to explore the mechanisms of new capillary sprout formation. In this study, using an integrative method, the interplay between the matrix density and angiogenesis was investigated. Owing the fact that the extracellular matrix density is a global parameter that can affect other parameters such as pore size, stiffness, cell-matrix adhesion and cross-linking, deeper understanding of the most important biomechanical or biochemical properties of the ECM causing changes in sprout morphogenesis is crucial. Here, we implemented both computational and experimental methods to analyze the mechanisms responsible for the influence of ECM density on the sprout formation that is difficult to be investigated comprehensively using each of these single methods. For this purpose, we first utilized an innovative approach to quantify the correspondence of the simulated collagen fibril density to the collagen density in the experimental part. Comparing the results of the experimental study and computational model led to some considerable achievements. First, we verified the results of the computational model using the experimental results. Then, we reported parameters such as the ratio of proliferating cells to migrating cells that was difficult to obtain from experimental study. Finally, this integrative system led to gain an understanding of the possible mechanisms responsible for the effect of ECM density on angiogenesis. The results showed that stable and long sprouts were observed at an intermediate collagen matrix density of 1.2 and 1.9 mg/ml due to a balance between the number of migrating and proliferating

Angular correlation of autoionization electrons and photons emitted from collisionally aligned atomic states

International Nuclear Information System (INIS)

Eichler, J.; Fritsch, W.

1976-01-01

The angular correlation of autoionization electrons or of photons ejected from collisionally aligned excited atoms is calculated assuming unpolarized beam and target, and polarization-insensitive detectors. Starting from the two-step hypothesis for the formation and decay of the intermediate excited atoms, the angular correlation is expressed in terms of the density matrix describing the excited system. Using the symmetries of the density matrix, a minimal set of independent matrix elements is given and the conditions for which a complete determination of this set is experimentally possible are discussed. For the case of electron emission, simple examples are pointed out in which the angular correlation is independent of the reduced Coulomb matrix elements describing the decay. (author)

Local-scaling density-functional method: Intraorbit and interorbit density optimizations

International Nuclear Information System (INIS)

Koga, T.; Yamamoto, Y.; Ludena, E.V.

1991-01-01

The recently proposed local-scaling density-functional theory provides us with a practical method for the direct variational determination of the electron density function ρ(r). The structure of ''orbits,'' which ensures the one-to-one correspondence between the electron density ρ(r) and the N-electron wave function Ψ({r k }), is studied in detail. For the realization of the local-scaling density-functional calculations, procedures for intraorbit and interorbit optimizations of the electron density function are proposed. These procedures are numerically illustrated for the helium atom in its ground state at the beyond-Hartree-Fock level

Conditional density matrix: systems and subsystems in quantum mechanics

International Nuclear Information System (INIS)

Belokurov, V.V.; Khrustalev, O.A.; Sadovnichij, V.A.; Timofeevskaya, O.D.

2003-01-01

A new quantum mechanical notion - Conditional Density Matrix - is discussed and is applied to describe some physical processes. This notion is a natural generalization of von Neumann density matrix for such processes as divisions of quantum systems into subsystems and reunifications of subsystems into new joint systems. Conditional Density Matrix assigns a quantum state to a subsystem of a composite system on condition that another part of the composite system is in some pure state

Determination of electron temperature and electron density in ...

African Journals Online (AJOL)

It is seen that the electron temperature increases from 5.8 × 102 oK to 7.83 × 104 oK as the pd is reduced from 130mm Hg × mm to 60 mm Hg × mm for argon. The electron densities increases from 2.8 × 1011/cm3 to 3.2 × 1011 /cm3 for the same variation of pds. For air the electron temperature increases from 3.6 × 104 oK to ...

A real-time extension of density matrix embedding theory for non-equilibrium electron dynamics

Science.gov (United States)

Kretchmer, Joshua S.; Chan, Garnet Kin-Lic

2018-02-01

We introduce real-time density matrix embedding theory (DMET), a dynamical quantum embedding theory for computing non-equilibrium electron dynamics in strongly correlated systems. As in the previously developed static DMET, real-time DMET partitions the system into an impurity corresponding to the region of interest coupled to the surrounding environment, which is efficiently represented by a quantum bath of the same size as the impurity. In this work, we focus on a simplified single-impurity time-dependent formulation as a first step toward a multi-impurity theory. The equations of motion of the coupled impurity and bath embedding problem are derived using the time-dependent variational principle. The accuracy of real-time DMET is compared to that of time-dependent complete active space self-consistent field (TD-CASSCF) theory and time-dependent Hartree-Fock (TDHF) theory for a variety of quantum quenches in the single impurity Anderson model (SIAM), in which the Hamiltonian is suddenly changed (quenched) to induce a non-equilibrium state. Real-time DMET shows a marked improvement over the mean-field TDHF, converging to the exact answer even in the non-trivial Kondo regime of the SIAM. However, as expected from analogous behavior in static DMET, the constrained structure of the real-time DMET wavefunction leads to a slower convergence with respect to active space size, in the single-impurity formulation, relative to TD-CASSCF. Our initial results suggest that real-time DMET provides a promising framework to simulate non-equilibrium electron dynamics in which strong electron correlation plays an important role, and lays the groundwork for future multi-impurity formulations.

Path integrals for electronic densities, reactivity indices, and localization functions in quantum systems.

Science.gov (United States)

Putz, Mihai V

2009-11-10

The density matrix theory, the ancestor of density functional theory, provides the immediate framework for Path Integral (PI) development, allowing the canonical density be extended for the many-electronic systems through the density functional closure relationship. Yet, the use of path integral formalism for electronic density prescription presents several advantages: assures the inner quantum mechanical description of the system by parameterized paths; averages the quantum fluctuations; behaves as the propagator for time-space evolution of quantum information; resembles Schrödinger equation; allows quantum statistical description of the system through partition function computing. In this framework, four levels of path integral formalism were presented: the Feynman quantum mechanical, the semiclassical, the Feynman-Kleinert effective classical, and the Fokker-Planck non-equilibrium ones. In each case the density matrix or/and the canonical density were rigorously defined and presented. The practical specializations for quantum free and harmonic motions, for statistical high and low temperature limits, the smearing justification for the Bohr's quantum stability postulate with the paradigmatic Hydrogen atomic excursion, along the quantum chemical calculation of semiclassical electronegativity and hardness, of chemical action and Mulliken electronegativity, as well as by the Markovian generalizations of Becke-Edgecombe electronic focalization functions - all advocate for the reliability of assuming PI formalism of quantum mechanics as a versatile one, suited for analytically and/or computationally modeling of a variety of fundamental physical and chemical reactivity concepts characterizing the (density driving) many-electronic systems.

Path Integrals for Electronic Densities, Reactivity Indices, and Localization Functions in Quantum Systems

Directory of Open Access Journals (Sweden)

Mihai V. Putz

2009-11-01

Full Text Available The density matrix theory, the ancestor of density functional theory, provides the immediate framework for Path Integral (PI development, allowing the canonical density be extended for the many-electronic systems through the density functional closure relationship. Yet, the use of path integral formalism for electronic density prescription presents several advantages: assures the inner quantum mechanical description of the system by parameterized paths; averages the quantum fluctuations; behaves as the propagator for time-space evolution of quantum information; resembles Schrödinger equation; allows quantum statistical description of the system through partition function computing. In this framework, four levels of path integral formalism were presented: the Feynman quantum mechanical, the semiclassical, the Feynman-Kleinert effective classical, and the Fokker-Planck non-equilibrium ones. In each case the density matrix or/and the canonical density were rigorously defined and presented. The practical specializations for quantum free and harmonic motions, for statistical high and low temperature limits, the smearing justification for the Bohr’s quantum stability postulate with the paradigmatic Hydrogen atomic excursion, along the quantum chemical calculation of semiclassical electronegativity and hardness, of chemical action and Mulliken electronegativity, as well as by the Markovian generalizations of Becke-Edgecombe electronic focalization functions – all advocate for the reliability of assuming PI formalism of quantum mechanics as a versatile one, suited for analytically and/or computationally modeling of a variety of fundamental physical and chemical reactivity concepts characterizing the (density driving many-electronic systems.

Decoherence dynamics of two charge qubits in vertically coupled quantum dots

International Nuclear Information System (INIS)

Ben Chouikha, W.; Bennaceur, R.; Jaziri, S.

2007-01-01

The decoherence dynamics of two charge qubits in a double quantum dot is investigated theoretically. We consider the quantum dynamics of two interacting electrons in a vertically coupled quantum dot driven by an external electric field. We derive the equations of motion for the density matrix, in which the presence of an electron confined in the double dot represents one qubit. A Markovian approach to the dynamical evolution of the reduced density matrix is adopted. We evaluate the concurrence of two qubits in order to study the effect of acoustic phonons on the entanglement. We also show that the disentanglement effect depends on the double dot parameters and increases with the temperature

A method to determine exactly the effective atomic number, electron density and absorbtion coefficient of materials from two Computer-Tomography measurement

International Nuclear Information System (INIS)

Christ, G.

1981-01-01

By the method of computer tomography, which is in use since about 10 years, X-ray images of a layer of interest can be produced without interference from the material present above this layer. An integral measurement of the attenuation of continuous X-radiation is sufficient to record the different attenuation behaviour in a layer for the purpose of image formation. For more information, however, can be obtained by taking into account the spectral distribution of the X-ray source and the energy dependence of the attenuation, which varies for different materials. In the experimental part of this work the measurement of the spectral distribution is described together with the necessary corrections, and the possible application of the cross sections for the relevant interaction processes, which are known from the literature is studied. As shown in the theoretical part, the attenuation coefficient can be described by an effective atomic number and the electron density of the absorber in the case of an arbitrary mixture of absorbing materials and a continuous X-ray spectrum. These two unknown material parameters can be determined by a method based on the measurement of two spectra with different spectral distribution. This is demonstrated by a one-dimensional and a two-dimensional computer simulation. (orig./WU) [de

Iterative approach as alternative to S-matrix in modal methods

Science.gov (United States)

Semenikhin, Igor; Zanuccoli, Mauro

2014-12-01

The continuously increasing complexity of opto-electronic devices and the rising demands of simulation accuracy lead to the need of solving very large systems of linear equations making iterative methods promising and attractive from the computational point of view with respect to direct methods. In particular, iterative approach potentially enables the reduction of required computational time to solve Maxwell's equations by Eigenmode Expansion algorithms. Regardless of the particular eigenmodes finding method used, the expansion coefficients are computed as a rule by scattering matrix (S-matrix) approach or similar techniques requiring order of M3 operations. In this work we consider alternatives to the S-matrix technique which are based on pure iterative or mixed direct-iterative approaches. The possibility to diminish the impact of M3 -order calculations to overall time and in some cases even to reduce the number of arithmetic operations to M2 by applying iterative techniques are discussed. Numerical results are illustrated to discuss validity and potentiality of the proposed approaches.

Two-Level Chebyshev Filter Based Complementary Subspace Method: Pushing the Envelope of Large-Scale Electronic Structure Calculations.

Science.gov (United States)

Banerjee, Amartya S; Lin, Lin; Suryanarayana, Phanish; Yang, Chao; Pask, John E

2018-06-12

We describe a novel iterative strategy for Kohn-Sham density functional theory calculations aimed at large systems (>1,000 electrons), applicable to metals and insulators alike. In lieu of explicit diagonalization of the Kohn-Sham Hamiltonian on every self-consistent field (SCF) iteration, we employ a two-level Chebyshev polynomial filter based complementary subspace strategy to (1) compute a set of vectors that span the occupied subspace of the Hamiltonian; (2) reduce subspace diagonalization to just partially occupied states; and (3) obtain those states in an efficient, scalable manner via an inner Chebyshev filter iteration. By reducing the necessary computation to just partially occupied states and obtaining these through an inner Chebyshev iteration, our approach reduces the cost of large metallic calculations significantly, while eliminating subspace diagonalization for insulating systems altogether. We describe the implementation of the method within the framework of the discontinuous Galerkin (DG) electronic structure method and show that this results in a computational scheme that can effectively tackle bulk and nano systems containing tens of thousands of electrons, with chemical accuracy, within a few minutes or less of wall clock time per SCF iteration on large-scale computing platforms. We anticipate that our method will be instrumental in pushing the envelope of large-scale ab initio molecular dynamics. As a demonstration of this, we simulate a bulk silicon system containing 8,000 atoms at finite temperature, and obtain an average SCF step wall time of 51 s on 34,560 processors; thus allowing us to carry out 1.0 ps of ab initio molecular dynamics in approximately 28 h (of wall time).

Density-matrix-functional calculations for matter in strong magnetic fields: Ground states of heavy atoms

DEFF Research Database (Denmark)

Johnsen, Kristinn; Yngvason, Jakob

1996-01-01

We report on a numerical study of the density matrix functional introduced by Lieb, Solovej, and Yngvason for the investigation of heavy atoms in high magnetic fields. This functional describes exactly the quantum mechanical ground state of atoms and ions in the limit when the nuclear charge Z...... and the electron number N tend to infinity with N/Z fixed, and the magnetic field B tends to infinity in such a way that B/Z4/3→∞. We have calculated electronic density profiles and ground-state energies for values of the parameters that prevail on neutron star surfaces and compared them with results obtained...... by other methods. For iron at B=1012 G the ground-state energy differs by less than 2% from the Hartree-Fock value. We have also studied the maximal negative ionization of heavy atoms in this model at various field strengths. In contrast to Thomas-Fermi type theories atoms can bind excess negative charge...

Application of the R-matrix method to photoionization of molecules.

Science.gov (United States)

Tashiro, Motomichi

2010-04-07

The R-matrix method has been used for theoretical calculation of electron collision with atoms and molecules for long years. The method was also formulated to treat photoionization process, however, its application has been mostly limited to photoionization of atoms. In this work, we implement the R-matrix method to treat molecular photoionization problem based on the UK R-matrix codes. This method can be used for diatomic as well as polyatomic molecules, with multiconfigurational description for electronic states of both target neutral molecule and product molecular ion. Test calculations were performed for valence electron photoionization of nitrogen (N(2)) as well as nitric oxide (NO) molecules. Calculated photoionization cross sections and asymmetry parameters agree reasonably well with the available experimental results, suggesting usefulness of the method for molecular photoionization.

Analysis of self-consistency effects in range-separated density-functional theory with Møller-Plesset perturbation theory

DEFF Research Database (Denmark)

Fromager, Emmanuel; Jensen, Hans Jørgen Aagaard

2011-01-01

Range-separated density-functional theory combines wave function theory for the long-range part of the two-electron interaction with density-functional theory for the short-range part. When describing the long-range interaction with non-variational methods, such as perturbation or coupled......-cluster theories, self-consistency effects are introduced in the density functional part, which for an exact solution requires iterations. They are generally assumed to be small but no detailed study has been performed so far. Here, the authors analyze self-consistency when using Møller-Plesset-type (MP......) perturbation theory for the long range interaction. The lowest-order self-consistency corrections to the wave function and the energy, that enter the perturbation expansions at the second and fourth order, respectively, are both expressed in terms of the one-electron reduced density matrix. The computational...

Electronic annealing Fermi operator expansion for DFT calculations on metallic systems

Science.gov (United States)

Aarons, Jolyon; Skylaris, Chris-Kriton

2018-02-01

Density Functional Theory (DFT) calculations with computational effort which increases linearly with the number of atoms (linear-scaling DFT) have been successfully developed for insulators, taking advantage of the exponential decay of the one-particle density matrix. For metallic systems, the density matrix is also expected to decay exponentially at finite electronic temperature and linear-scaling DFT methods should be possible by taking advantage of this decay. Here we present a method for DFT calculations at finite electronic temperature for metallic systems which is effectively linear-scaling (O(N)). Our method generates the elements of the one-particle density matrix and also finds the required chemical potential and electronic entropy using polynomial expansions. A fixed expansion length is always employed to generate the density matrix, without any loss in accuracy by the application of a high electronic temperature followed by successive steps of temperature reduction until the desired (low) temperature density matrix is obtained. We have implemented this method in the ONETEP linear-scaling (for insulators) DFT code which employs local orbitals that are optimised in situ. By making use of the sparse matrix machinery of ONETEP, our method exploits the sparsity of Hamiltonian and density matrices to perform calculations on metallic systems with computational cost that increases asymptotically linearly with the number of atoms. We demonstrate the linear-scaling computational cost of our method with calculation times on palladium nanoparticles with up to ˜13 000 atoms.

Stationary solution of a time dependent density matrix formalism

International Nuclear Information System (INIS)

Tohyama, Mitsuru

1994-01-01

A stationary solution of a time-dependent density-matrix formalism, which is an extension of the time-dependent Hartree-Fock theory to include the effects of two-body correlations, is obtained for the Lipkin model hamiltonian, using an adiabatic treatment of the two-body interaction. It is found that the obtained result is a reasonable approximation for the exact solution of the model. (author)

A method for conversion of Hounsfield number to electron density and prediction of macroscopic pair production cross-sections

International Nuclear Information System (INIS)

Knoeoes, T.; Nilsson, M.; Ahlgren, L.

1986-01-01

A method for the determination of electron density using a narrow beam attenuation geometry is described. The method does not require that the elemental composition of the phantom materials is known. The Hounsfield numbers for the phantom materials used were determined using five different CT scanners. A relationship between Hounsfield number and electron density can thus be established, which is of considerable value in radiation therapy treatment planning procedures. Measurements of the ratio coherent/incoherent scattering of low energy photons in a certain geometry has proven valuable for determination of atomic number, which in its turn can be used for estimation of macroscopic pair production coefficients for high energy photons. The combination of knowledge of electron density with methods for determination of processes, dependent on atomic number, can form a base for adequate composition of phantom materials for purposes of testing dose calculation algorithms for photons and electrons. (orig.)

Device and method for relativistic electron beam heating of a high-density plasma to drive fast liners

International Nuclear Information System (INIS)

Thode, L.E.

1981-01-01

A device and method for relativistic electron beam heating of a high-density plasma in a small localized region are described. A relativistic electron beam generator or accelerator produces a high-voltage electron beam which propagates along a vacuum drift tube and is modulated to initiate electron bunching within the beam. The beam is then directed through a low-density gas chamber which provides isolation between the vacuum modulator and the relativistic electron beam target. The relativistic beam is then applied to a high-density target plasma which typically comprises dt, dd, hydrogen boron or similar thermonuclear gas at a density of 1017 to 1020 electrons per cubic centimeter. The target gas is ionized prior to application of the electron beam by means of a laser or other preionization source to form a plasma. Utilizing a relativistic electron beam with an individual particle energy exceeding 3 mev, classical scattering by relativistic electrons passing through isolation foils is negligible. As a result, relativistic streaming instabilities are initiated within the high-density target plasma causing the relativistic electron beam to efficiently deposit its energy and momentum into a small localized region of the high-density plasma target. Fast liners disposed in the high-density target plasma are explosively or ablatively driven to implosion by a heated annular plasma surrounding the fast liner which is generated by an annular relativistic electron beam. An azimuthal magnetic field produced by axial current flow in the annular plasma, causes the energy in the heated annular plasma to converge on the fast liner

Rapid model building of beta-sheets in electron-density maps.

Science.gov (United States)

Terwilliger, Thomas C

2010-03-01

A method for rapidly building beta-sheets into electron-density maps is presented. beta-Strands are identified as tubes of high density adjacent to and nearly parallel to other tubes of density. The alignment and direction of each strand are identified from the pattern of high density corresponding to carbonyl and C(beta) atoms along the strand averaged over all repeats present in the strand. The beta-strands obtained are then assembled into a single atomic model of the beta-sheet regions. The method was tested on a set of 42 experimental electron-density maps at resolutions ranging from 1.5 to 3.8 A. The beta-sheet regions were nearly completely built in all but two cases, the exceptions being one structure at 2.5 A resolution in which a third of the residues in beta-sheets were built and a structure at 3.8 A in which under 10% were built. The overall average r.m.s.d. of main-chain atoms in the residues built using this method compared with refined models of the structures was 1.5 A.

Electron density profile reconstruction by maximum entropy method with multichannel HCN laser interferometer system on SPAC VII

International Nuclear Information System (INIS)

Kubo, S.; Narihara, K.; Tomita, Y.; Hasegawa, M.; Tsuzuki, T.; Mohri, A.

1988-01-01

A multichannel HCN laser interferometer system has been developed to investigate the plasma electron confinement properties in SPAC VII device. Maximum entropy method is applied to reconstruct the electron density profile from measured line integrated data. Particle diffusion coefficient in the peripheral region of the REB ring core spherator was obtained from the evolution of the density profile. (author)

Spatial charge motion on an uniform density matrix-general equations in opened and closed circuits

International Nuclear Information System (INIS)

Aguiar Monsanto, S. de.

1983-01-01

The motion of a space charge cloud embedded in a matrix of constant immobile charge density is studied in open as well as in closed circuit. In the first case, open circuit, the solution is almost trivial as compared as the other one in which, after some work, the problem is reduced to an ordinary differential equation. The method of solution is parallel to that employed in the study of monopolar free space charge motion. The voltage and the current produced by a system with no net charge but with unbalanced local charge density were calculated using the general equations derived in the first part of the work. (Author) [pt

Regular perturbation theory for two-electron atoms

International Nuclear Information System (INIS)

Feranchuk, I.D.; Triguk, V.V.

2011-01-01

Regular perturbation theory (RPT) for the ground and excited states of two-electron atoms or ions is developed. It is shown for the first time that summation of the matrix elements from the electron-electron interaction operator over all intermediate states can be calculated in a closed form by means of the two-particle Coulomb Green's function constructed in the Letter. It is shown that the second order approximation of RPT includes the main part of the correlation energy both for the ground and excited states. This approach can be also useful for description of two-electron atoms in external fields. -- Highlights: → We develop regular perturbation theory for the two-electron atoms or ions. → We calculate the sum of the matrix elements over all intermediate states. → We construct the two-particle Coulomb Green's function.

Edge electron density profiles and fluctuations measured by two-dimensional beam emission spectroscopy in the KSTAR

Energy Technology Data Exchange (ETDEWEB)

Nam, Y. U., E-mail: yunam@nfri.re.kr; Wi, H. M. [National Fusion Research Institute, Daejeon (Korea, Republic of); Zoletnik, S.; Lampert, M. [Wigner RCP Institute for Particle and Nuclear Physics, Budapest (Hungary); Kovácsik, Ákos [Institute of Nuclear Techniques, Budapest Technical University, Budapest (Hungary)

2014-11-15

Beam emission spectroscopy (BES) system in Korea Superconducting Tokamak Advanced Research (KSTAR) has recently been upgraded. The background intensity was reduced from 30% to 2% by suppressing the stray lights. This allows acquisition of the relative electron density profiles on the plasma edge without background subtraction from the beam power modulation signals. The KSTAR BES system has its spatial resolution of 1 cm, the temporal resolution of 2 MHz, and a total 32 channel (8 radial × 4 poloidal) avalanche photo diode array. Most measurements were done on the plasma edge, r/a ∼ 0.9, with 8 cm radial measurement width that covers the pedestal range. High speed density profile measurements reveal temporal behaviors of fast transient events, such as the precursors of edge localized modes and the transitions between confinement modes. Low background level also allows analysis of the edge density fluctuation patterns with reduced background fluctuations. Propagation of the density structures can be investigated by comparing the phase delays between the spatially distributed channels.

Electron density and gas density measurements in a millimeter-wave discharge

Energy Technology Data Exchange (ETDEWEB)

Schaub, S. C., E-mail: sschaub@mit.edu; Hummelt, J. S.; Guss, W. C.; Shapiro, M. A.; Temkin, R. J. [Plasma Science and Fusion Center, Massachusetts Institute of Technology 167 Albany St., Bldg. NW16, Cambridge, Massachusetts 02139 (United States)

2016-08-15

Electron density and neutral gas density have been measured in a non-equilibrium air breakdown plasma using optical emission spectroscopy and two-dimensional laser interferometry, respectively. A plasma was created with a focused high frequency microwave beam in air. Experiments were run with 110 GHz and 124.5 GHz microwaves at powers up to 1.2 MW. Microwave pulses were 3 μs long at 110 GHz and 2.2 μs long at 124.5 GHz. Electron density was measured over a pressure range of 25 to 700 Torr as the input microwave power was varied. Electron density was found to be close to the critical density, where the collisional plasma frequency is equal to the microwave frequency, over the pressure range studied and to vary weakly with input power. Neutral gas density was measured over a pressure range from 150 to 750 Torr at power levels high above the threshold for initiating breakdown. The two-dimensional structure of the neutral gas density was resolved. Intense, localized heating was found to occur hundreds of nanoseconds after visible plasma formed. This heating led to neutral gas density reductions of greater than 80% where peak plasma densities occurred. Spatial structure and temporal dynamics of gas heating at atmospheric pressure were found to agree well with published numerical simulations.

Electron density and gas density measurements in a millimeter-wave discharge

International Nuclear Information System (INIS)

Schaub, S. C.; Hummelt, J. S.; Guss, W. C.; Shapiro, M. A.; Temkin, R. J.

2016-01-01

Electron density and neutral gas density have been measured in a non-equilibrium air breakdown plasma using optical emission spectroscopy and two-dimensional laser interferometry, respectively. A plasma was created with a focused high frequency microwave beam in air. Experiments were run with 110 GHz and 124.5 GHz microwaves at powers up to 1.2 MW. Microwave pulses were 3 μs long at 110 GHz and 2.2 μs long at 124.5 GHz. Electron density was measured over a pressure range of 25 to 700 Torr as the input microwave power was varied. Electron density was found to be close to the critical density, where the collisional plasma frequency is equal to the microwave frequency, over the pressure range studied and to vary weakly with input power. Neutral gas density was measured over a pressure range from 150 to 750 Torr at power levels high above the threshold for initiating breakdown. The two-dimensional structure of the neutral gas density was resolved. Intense, localized heating was found to occur hundreds of nanoseconds after visible plasma formed. This heating led to neutral gas density reductions of greater than 80% where peak plasma densities occurred. Spatial structure and temporal dynamics of gas heating at atmospheric pressure were found to agree well with published numerical simulations.

Multireference configuration interaction theory using cumulant reconstruction with internal contraction of density matrix renormalization group wave function.

Science.gov (United States)

Saitow, Masaaki; Kurashige, Yuki; Yanai, Takeshi

2013-07-28

We report development of the multireference configuration interaction (MRCI) method that can use active space scalable to much larger size references than has previously been possible. The recent development of the density matrix renormalization group (DMRG) method in multireference quantum chemistry offers the ability to describe static correlation in a large active space. The present MRCI method provides a critical correction to the DMRG reference by including high-level dynamic correlation through the CI treatment. When the DMRG and MRCI theories are combined (DMRG-MRCI), the full internal contraction of the reference in the MRCI ansatz, including contraction of semi-internal states, plays a central role. However, it is thought to involve formidable complexity because of the presence of the five-particle rank reduced-density matrix (RDM) in the Hamiltonian matrix elements. To address this complexity, we express the Hamiltonian matrix using commutators, which allows the five-particle rank RDM to be canceled out without any approximation. Then we introduce an approximation to the four-particle rank RDM by using a cumulant reconstruction from lower-particle rank RDMs. A computer-aided approach is employed to derive the exceedingly complex equations of the MRCI in tensor-contracted form and to implement them into an efficient parallel computer code. This approach extends to the size-consistency-corrected variants of MRCI, such as the MRCI+Q, MR-ACPF, and MR-AQCC methods. We demonstrate the capability of the DMRG-MRCI method in several benchmark applications, including the evaluation of single-triplet gap of free-base porphyrin using 24 active orbitals.

Intrinsic Density Matrices of the Nuclear Shell Model

International Nuclear Information System (INIS)

Deveikis, A.; Kamuntavichius, G.

1996-01-01

A new method for calculation of shell model intrinsic density matrices, defined as two-particle density matrices integrated over the centre-of-mass position vector of two last particles and complemented with isospin variables, has been developed. The intrinsic density matrices obtained are completely antisymmetric, translation-invariant, and do not employ a group-theoretical classification of antisymmetric states. They are used for exact realistic density matrix expansion within the framework of the reduced Hamiltonian method. The procedures based on precise arithmetic for calculation of the intrinsic density matrices that involve no numerical diagonalization or orthogonalization have been developed and implemented in the computer code. (author). 11 refs., 2 tabs

New method of ionization energy calculation for two-electron ions

International Nuclear Information System (INIS)

Ershov, D.K.

1997-01-01

A new method for calculation of the ionization energy of two-electron ions is proposed. The method is based on the calculation of the energy of second electron interaction with the field of an one-electron ion the potential of which is well known

How electronic dynamics with Pauli exclusion produces Fermi-Dirac statistics.

Science.gov (United States)

Nguyen, Triet S; Nanguneri, Ravindra; Parkhill, John

2015-04-07

It is important that any dynamics method approaches the correct population distribution at long times. In this paper, we derive a one-body reduced density matrix dynamics for electrons in energetic contact with a bath. We obtain a remarkable equation of motion which shows that in order to reach equilibrium properly, rates of electron transitions depend on the density matrix. Even though the bath drives the electrons towards a Boltzmann distribution, hole blocking factors in our equation of motion cause the electronic populations to relax to a Fermi-Dirac distribution. These factors are an old concept, but we show how they can be derived with a combination of time-dependent perturbation theory and the extended normal ordering of Mukherjee and Kutzelnigg for a general electronic state. The resulting non-equilibrium kinetic equations generalize the usual Redfield theory to many-electron systems, while ensuring that the orbital occupations remain between zero and one. In numerical applications of our equations, we show that relaxation rates of molecules are not constant because of the blocking effect. Other applications to model atomic chains are also presented which highlight the importance of treating both dephasing and relaxation. Finally, we show how the bath localizes the electron density matrix.

How electronic dynamics with Pauli exclusion produces Fermi-Dirac statistics

International Nuclear Information System (INIS)

Nguyen, Triet S.; Nanguneri, Ravindra; Parkhill, John

2015-01-01

It is important that any dynamics method approaches the correct population distribution at long times. In this paper, we derive a one-body reduced density matrix dynamics for electrons in energetic contact with a bath. We obtain a remarkable equation of motion which shows that in order to reach equilibrium properly, rates of electron transitions depend on the density matrix. Even though the bath drives the electrons towards a Boltzmann distribution, hole blocking factors in our equation of motion cause the electronic populations to relax to a Fermi-Dirac distribution. These factors are an old concept, but we show how they can be derived with a combination of time-dependent perturbation theory and the extended normal ordering of Mukherjee and Kutzelnigg for a general electronic state. The resulting non-equilibrium kinetic equations generalize the usual Redfield theory to many-electron systems, while ensuring that the orbital occupations remain between zero and one. In numerical applications of our equations, we show that relaxation rates of molecules are not constant because of the blocking effect. Other applications to model atomic chains are also presented which highlight the importance of treating both dephasing and relaxation. Finally, we show how the bath localizes the electron density matrix

Matrix Methods for Solving Hartree-Fock Equations in Atomic Structure Calculations and Line Broadening

Directory of Open Access Journals (Sweden)

Thomas Gomez

2018-04-01

Full Text Available Atomic structure of N-electron atoms is often determined by solving the Hartree-Fock equations, which are a set of integro-differential equations. The integral part of the Hartree-Fock equations treats electron exchange, but the Hartree-Fock equations are not often treated as an integro-differential equation. The exchange term is often approximated as an inhomogeneous or an effective potential so that the Hartree-Fock equations become a set of ordinary differential equations (which can be solved using the usual shooting methods. Because the Hartree-Fock equations are an iterative-refinement method, the inhomogeneous term relies on the previous guess of the wavefunction. In addition, there are numerical complications associated with solving inhomogeneous differential equations. This work uses matrix methods to solve the Hartree-Fock equations as an integro-differential equation. It is well known that a derivative operator can be expressed as a matrix made of finite-difference coefficients; energy eigenvalues and eigenvectors can be obtained by using linear-algebra packages. The integral (exchange part of the Hartree-Fock equation can be approximated as a sum and written as a matrix. The Hartree-Fock equations can be solved as a matrix that is the sum of the differential and integral matrices. We compare calculations using this method against experiment and standard atomic structure calculations. This matrix method can also be used to solve for free-electron wavefunctions, thus improving how the atoms and free electrons interact. This technique is important for spectral line broadening in two ways: it improves the atomic structure calculations, and it improves the motion of the plasma electrons that collide with the atom.

Efficient density matrix renormalization group algorithm to study Y junctions with integer and half-integer spin

KAUST Repository

Kumar, Manoranjan

2016-02-03

An efficient density matrix renormalization group (DMRG) algorithm is presented and applied to Y junctions, systems with three arms of n sites that meet at a central site. The accuracy is comparable to DMRG of chains. As in chains, new sites are always bonded to the most recently added sites and the superblock Hamiltonian contains only new or once renormalized operators. Junctions of up to N=3n+1≈500 sites are studied with antiferromagnetic (AF) Heisenberg exchange J between nearest-neighbor spins S or electron transfer t between nearest neighbors in half-filled Hubbard models. Exchange or electron transfer is exclusively between sites in two sublattices with NA≠NB. The ground state (GS) and spin densities ρr=⟨Szr⟩ at site r are quite different for junctions with S=1/2, 1, 3/2, and 2. The GS has finite total spin SG=2S(S) for even (odd) N and for MG=SG in the SG spin manifold, ρr>0(<0) at sites of the larger (smaller) sublattice. S=1/2 junctions have delocalized states and decreasing spin densities with increasing N. S=1 junctions have four localized Sz=1/2 states at the end of each arm and centered on the junction, consistent with localized states in S=1 chains with finite Haldane gap. The GS of S=3/2 or 2 junctions of up to 500 spins is a spin density wave with increased amplitude at the ends of arms or near the junction. Quantum fluctuations completely suppress AF order in S=1/2 or 1 junctions, as well as in half-filled Hubbard junctions, but reduce rather than suppress AF order in S=3/2 or 2 junctions.

Efficient density matrix renormalization group algorithm to study Y junctions with integer and half-integer spin

KAUST Repository

Kumar, Manoranjan; Parvej, Aslam; Thomas, Simil; Ramasesha, S.; Soos, Z. G.

2016-01-01

An efficient density matrix renormalization group (DMRG) algorithm is presented and applied to Y junctions, systems with three arms of n sites that meet at a central site. The accuracy is comparable to DMRG of chains. As in chains, new sites are always bonded to the most recently added sites and the superblock Hamiltonian contains only new or once renormalized operators. Junctions of up to N=3n+1≈500 sites are studied with antiferromagnetic (AF) Heisenberg exchange J between nearest-neighbor spins S or electron transfer t between nearest neighbors in half-filled Hubbard models. Exchange or electron transfer is exclusively between sites in two sublattices with NA≠NB. The ground state (GS) and spin densities ρr=⟨Szr⟩ at site r are quite different for junctions with S=1/2, 1, 3/2, and 2. The GS has finite total spin SG=2S(S) for even (odd) N and for MG=SG in the SG spin manifold, ρr>0(<0) at sites of the larger (smaller) sublattice. S=1/2 junctions have delocalized states and decreasing spin densities with increasing N. S=1 junctions have four localized Sz=1/2 states at the end of each arm and centered on the junction, consistent with localized states in S=1 chains with finite Haldane gap. The GS of S=3/2 or 2 junctions of up to 500 spins is a spin density wave with increased amplitude at the ends of arms or near the junction. Quantum fluctuations completely suppress AF order in S=1/2 or 1 junctions, as well as in half-filled Hubbard junctions, but reduce rather than suppress AF order in S=3/2 or 2 junctions.

Tight-binding approximations to time-dependent density functional theory — A fast approach for the calculation of electronically excited states

Energy Technology Data Exchange (ETDEWEB)

Rüger, Robert, E-mail: rueger@scm.com [Scientific Computing & Modelling NV, De Boelelaan 1083, 1081 HV Amsterdam (Netherlands); Department of Theoretical Chemistry, Vrije Universiteit Amsterdam, De Boelelaan 1083, 1081 HV Amsterdam (Netherlands); Wilhelm-Ostwald-Institut für Physikalische und Theoretische Chemie, Linnéstr. 2, 04103 Leipzig (Germany); Lenthe, Erik van [Scientific Computing & Modelling NV, De Boelelaan 1083, 1081 HV Amsterdam (Netherlands); Heine, Thomas [Wilhelm-Ostwald-Institut für Physikalische und Theoretische Chemie, Linnéstr. 2, 04103 Leipzig (Germany); Visscher, Lucas [Department of Theoretical Chemistry, Vrije Universiteit Amsterdam, De Boelelaan 1083, 1081 HV Amsterdam (Netherlands)

2016-05-14

We propose a new method of calculating electronically excited states that combines a density functional theory based ground state calculation with a linear response treatment that employs approximations used in the time-dependent density functional based tight binding (TD-DFTB) approach. The new method termed time-dependent density functional theory TD-DFT+TB does not rely on the DFTB parametrization and is therefore applicable to systems involving all combinations of elements. We show that the new method yields UV/Vis absorption spectra that are in excellent agreement with computationally much more expensive TD-DFT calculations. Errors in vertical excitation energies are reduced by a factor of two compared to TD-DFTB.

Effects of SiO2 nano-particles on tribological and mechanical properties of aluminum matrix composites by different dispersion methods

Science.gov (United States)

Azadi, Mahboobeh; Zolfaghari, Mehrdad; Rezanezhad, Saeid; Azadi, Mohammad

2018-05-01

This study has been presented with mechanical properties of aluminum matrix composites, reinforced by SiO2 nano-particles. The stir casting method was employed to produce various aluminum matrix composites. Different composites by varying the SiO2 nano-particle content (including 0.5 and 1 weight percents) and two dispersion methods (including ball-milling and pre-heating) were made. Then, the density, the hardness, the compression strength, the wear resistance and the microstructure of nano-composites have been studied in this research. Besides, the distribution of nano-particles in the aluminum matrix for all composites has been also evaluated by the field emission scanning electron microscopy (FESEM). Obtained results showed that the density, the elongation and the ultimate compressive strength of various nano-composites decreased by the presence of SiO2 nano-particles; however, the hardness, the wear resistance, the yield strength and the elastic modulus of composites increased by auditioning of nano-particles to the aluminum alloy. FESEM images indicated better wetting of the SiO2 reinforcement in the aluminum matrix, prepared by the pre-heating dispersion method, comparing to ball-milling. When SiO2 nano-particles were added to the aluminum alloy, the morphology of the Si phase and intermetallic phases changed, which enhanced mechanical properties. In addition, the wear mechanism plus the friction coefficient value were changed for various nano-composites with respect to the aluminum alloy.

Matrix-type multiple reciprocity boundary element method for solving three-dimensional two-group neutron diffusion equations

International Nuclear Information System (INIS)

Itagaki, Masafumi; Sahashi, Naoki.

1997-01-01

The multiple reciprocity boundary element method has been applied to three-dimensional two-group neutron diffusion problems. A matrix-type boundary integral equation has been derived to solve the first and the second group neutron diffusion equations simultaneously. The matrix-type fundamental solutions used here satisfy the equation which has a point source term and is adjoint to the neutron diffusion equations. A multiple reciprocity method has been employed to transform the matrix-type domain integral related to the fission source into an equivalent boundary one. The higher order fundamental solutions required for this formulation are composed of a series of two types of analytic functions. The eigenvalue itself is also calculated using only boundary integrals. Three-dimensional test calculations indicate that the present method provides stable and accurate solutions for criticality problems. (author)

Electron density interferometry measurement in laser-matter interaction

International Nuclear Information System (INIS)

Popovics-Chenais, C.

1981-05-01

This work is concerned with the laser-interferometry measurement of the electronic density in the corona and the conduction zone external part. Particularly, it is aimed at showing up density gradients and at their space-time localization. The first chapter recalls the density profile influence on the absorption principal mechanisms and the laser energy transport. In chapter two, the numerical and analytical hydrodynamic models describing the density profile are analysed. The influence on the density profile of the ponderomotive force associated to high oscillating electric fields is studied, together with the limited thermal conduction and suprathermal electron population. The mechanism action, in our measurement conditions, is numerically simulated. Calculations are made with experimental parameters. The measurement interaction conditions, together with the diagnostic method by high resolution laser interferometry are detailed. The results are analysed with the help of numerical simulation which is the experiment modeling. An overview of the mechanisms shown up by interferometric measurements and their correlation with other diagnostics is the conclusion of this work [fr

Investigation of the alpha cluster model and the density matrix expansion in ion-ion collision

International Nuclear Information System (INIS)

Rashdan, M.B.M.

1986-01-01

This thesis deals with the investigation of the alpha cluster model (ACM) of brink and studies of the accuracy of the density matrix expansion (DME) approximation in deriving the real part of the ion-ion optical potential. the ACM is applied to calculate the inelastic 0 1 + →2 1 + charge form factor for electron scattering by 12 C to investigate the validity of this model for 12 C nucleus. it is found that the experimental curve can be fitted over the entire range of the momentum transfer by a generator - coordinate state for the 2 1 + state that consist of a superposition of two triangular ACM states with two different cluster separations and the same oscillator parameter

Spin theory of the density functional: reduced matrices and density functions

International Nuclear Information System (INIS)

Pavlov, R.; Delchev, Y.; Pavlova, K.; Maruani, J.

1993-01-01

Expressions for the reduced matrices and density functions of N-fermion systems of arbitrary order s (1<=s<=N) are derived within the frame of rigorous spin approach to the density functional theory (DFT). Using the local-scale transformation method and taking into account the particle spin it is shown that the reduced matrices and density functions are functionals of the total one-fermion density. Similar dependence is found for the distribution density of s-particle aggregates. Generalization and applicability of DFT to the case of s-particle ensembles and aggregates is discussed. 14 refs

Multireference quantum chemistry through a joint density matrix renormalization group and canonical transformation theory.

Science.gov (United States)

Yanai, Takeshi; Kurashige, Yuki; Neuscamman, Eric; Chan, Garnet Kin-Lic

2010-01-14

We describe the joint application of the density matrix renormalization group and canonical transformation theory to multireference quantum chemistry. The density matrix renormalization group provides the ability to describe static correlation in large active spaces, while the canonical transformation theory provides a high-order description of the dynamic correlation effects. We demonstrate the joint theory in two benchmark systems designed to test the dynamic and static correlation capabilities of the methods, namely, (i) total correlation energies in long polyenes and (ii) the isomerization curve of the [Cu(2)O(2)](2+) core. The largest complete active spaces and atomic orbital basis sets treated by the joint DMRG-CT theory in these systems correspond to a (24e,24o) active space and 268 atomic orbitals in the polyenes and a (28e,32o) active space and 278 atomic orbitals in [Cu(2)O(2)](2+).

On the statistical interpretation of quantum mechanics: evolution of the density matrix

International Nuclear Information System (INIS)

Benzecri, J.P.

1986-01-01

Without attempting to identify ontological interpretation with a mathematical structure, we reduce philosophical speculation to five theses. In the discussion of these, a central role is devoted to the mathematical problem of the evolution of the density matrix. This article relates to the first 3 of these 5 theses [fr

Histogram specification as a method of density modification

International Nuclear Information System (INIS)

Harrison, R.W.

1988-01-01

A new method for improving the quality and extending the resolution of Fourier maps is described. The method is based on a histogram analysis of the electron density. The distribution of electron density values in the map is forced to be 'ideal'. The 'ideal' distribution is assumed to be Gaussian. The application of the method to improve the electron density map for the protein Acinetobacter asparaginase, which is a tetrameric enzyme of molecular weight 140000 daltons, is described. (orig.)

Experimental Electron Density Distribution in Two Cocrystals of Betaines with p-Hydroxybenzoic Acid

Directory of Open Access Journals (Sweden)

Agata Owczarzak

2018-03-01

Full Text Available Experimental determination of electron density distribution in crystals by means of high-resolution X-ray diffraction allows, among others, for studying the details of intra- and inter-molecular interactions. In case of co-crystals, this method may help in finding the conditions of creating such species. The results of such analysis for two co-crystals containing betaines, namely trigonelline (TRG: nicotinic acid N-methylbetaine, IUPAC name: 1-methylpyridinium-3-carboxylate and N-methylpiperidine betaine (MPB: 1-methylpiperidinium-1-yl-carboxylate with p-hydroxybenzoic acid (HBA are reported. TRG-HBA crystallizes as a hydrate. For both of the co-crystals, high-quality diffraction data were collected up to sinθ/λ = 1.13 Å−1. Hansen-Coppens multipolar model was then applied for modelling the electron density distribution and Atoms-In-Molecules approach was used for detailed analysis of interactions in crystals. A number of intermolecular interactions was identified, ranging from strong O-H···O hydrogen bonds through C-H···O to C-H···π and π···π interactions. Correlations between the geometrical characteristics of the contacts and the features of their critical points were analyzed in detail. Atomic charges show that in zwitterionic species there are regions of opposite charges, rather than charges that are localized on certain atoms. In case of MPB-HBA, a significant charge transfer between the components of co-crystal (0.5 e was found, as opposed to TRG-HBA, where all of the components are almost neutral.

Proposed non-interferometric FIR electron density measuring scheme for Tokamaks

Energy Technology Data Exchange (ETDEWEB)

Dodel, G; Kunz, W [Stuttgart Univ. (TH) (Germany, F.R.). Inst. fuer Plasmaforschung

1979-08-01

Extension of FIR polarimetry to electron density measurements in Tokamaks is suggested as a possible alternative for devices in which FIR interferometry is not applicable or difficult to handle due to reduced accessibility or strong mechanical vibrations. The method is numerically simulated. The relative experimental simplicity compared with interferometry has to be paid for with symmetry assumptions which enter into the evaluation process.

Interaction effects in liquids with low electron densities

International Nuclear Information System (INIS)

Warren, W.W. Jr.

1987-01-01

The author discusses two complementary classes of systems in which strong electron-electron or electron-ion interactions appear at low electron densities. The first are the expanded liquid alkali metals (cesium) in which electron correlation effects have a profound effect on the magnetic properties on the metallic side of the metal-nonmetal transition. The second group are molten alkali halides containing low densities of localized electrons introduced, say, by dissolution of small amounts of excess metal. (Auth.)

Response matrix Monte Carlo based on a general geometry local calculation for electron transport

International Nuclear Information System (INIS)

Ballinger, C.T.; Rathkopf, J.A.; Martin, W.R.

1991-01-01

A Response Matrix Monte Carlo (RMMC) method has been developed for solving electron transport problems. This method was born of the need to have a reliable, computationally efficient transport method for low energy electrons (below a few hundred keV) in all materials. Today, condensed history methods are used which reduce the computation time by modeling the combined effect of many collisions but fail at low energy because of the assumptions required to characterize the electron scattering. Analog Monte Carlo simulations are prohibitively expensive since electrons undergo coulombic scattering with little state change after a collision. The RMMC method attempts to combine the accuracy of an analog Monte Carlo simulation with the speed of the condensed history methods. Like condensed history, the RMMC method uses probability distributions functions (PDFs) to describe the energy and direction of the electron after several collisions. However, unlike the condensed history method the PDFs are based on an analog Monte Carlo simulation over a small region. Condensed history theories require assumptions about the electron scattering to derive the PDFs for direction and energy. Thus the RMMC method samples from PDFs which more accurately represent the electron random walk. Results show good agreement between the RMMC method and analog Monte Carlo. 13 refs., 8 figs

Matrix method for acoustic levitation simulation.

Science.gov (United States)

Andrade, Marco A B; Perez, Nicolas; Buiochi, Flavio; Adamowski, Julio C

2011-08-01

A matrix method is presented for simulating acoustic levitators. A typical acoustic levitator consists of an ultrasonic transducer and a reflector. The matrix method is used to determine the potential for acoustic radiation force that acts on a small sphere in the standing wave field produced by the levitator. The method is based on the Rayleigh integral and it takes into account the multiple reflections that occur between the transducer and the reflector. The potential for acoustic radiation force obtained by the matrix method is validated by comparing the matrix method results with those obtained by the finite element method when using an axisymmetric model of a single-axis acoustic levitator. After validation, the method is applied in the simulation of a noncontact manipulation system consisting of two 37.9-kHz Langevin-type transducers and a plane reflector. The manipulation system allows control of the horizontal position of a small levitated sphere from -6 mm to 6 mm, which is done by changing the phase difference between the two transducers. The horizontal position of the sphere predicted by the matrix method agrees with the horizontal positions measured experimentally with a charge-coupled device camera. The main advantage of the matrix method is that it allows simulation of non-symmetric acoustic levitators without requiring much computational effort.

TREATMENT OF NONADIABATIC TRANSITIONS BY DENSITY-MATRIX EVOLUTION AND MOLECULAR-DYNAMICS SIMULATIONS

NARCIS (Netherlands)

MAVRI, J; BERENDSEN, HJC

1994-01-01

A density matrix evolution (DME) method (H.J.C. Berendsen and J. Mavri, J. Phys. Chem., 97 (1993) 13469) to simulate the dynamics of quantum systems embedded in a classical environment is presented. The DME method allows treatment of nonadiabatic transitions. As numerical examples the collinear

Histogram specification as a method of density modification

Energy Technology Data Exchange (ETDEWEB)

Harrison, R.W.

1988-12-01

A new method for improving the quality and extending the resolution of Fourier maps is described. The method is based on a histogram analysis of the electron density. The distribution of electron density values in the map is forced to be 'ideal'. The 'ideal' distribution is assumed to be Gaussian. The application of the method to improve the electron density map for the protein Acinetobacter asparaginase, which is a tetrameric enzyme of molecular weight 140000 daltons, is described.

Factorization method for simulating QCD at finite density

International Nuclear Information System (INIS)

Nishimura, Jun

2003-01-01

We propose a new method for simulating QCD at finite density. The method is based on a general factorization property of distribution functions of observables, and it is therefore applicable to any system with a complex action. The so-called overlap problem is completely eliminated by the use of constrained simulations. We test this method in a Random Matrix Theory for finite density QCD, where we are able to reproduce the exact results for the quark number density. (author)

Homocomposites of chopped fluorinated polyethylene fiber with low-density polyethylene matrix

International Nuclear Information System (INIS)

Maity, J.; Jacob, C.; Das, C.K.; Alam, S.; Singh, R.P.

2008-01-01

Conventional composites are generally prepared by adding reinforcing agent to a matrix and the matrix wherein the reinforcing agents are different in chemical composition with the later having superior mechanical properties. This work presents the preparation and properties of homocomposites consisting of a low-density polyethylene (LDPE) matrix and an ultra high molecular weight polyethylene (UHMWPE) fiber reinforcing phase. Direct fluorination is an important surface modification process by which only a thin upper layer is modified, the bulk properties of the polymer remaining unchanged. In this work, surface fluorination of UHMWPE fiber was done and then fiber characterization was performed. It was observed that after fluorination the fiber surface became rough. Composites were then prepared using both fluorinated and non-fluorinated polyethylene fiber with a low-density polyethylene (LDPE) matrix to prepare single polymer composites. It was found that the thermal stability and mechanical properties were improved for fluorinated fiber composites. X-ray diffraction (XRD) analysis showed that the crystallinity of the composites increased and it is maximum for fluorinated fiber composites. Tensile strength (TS) and modulus also increased while elongation at break (EB) decreased for fiber composites and was a maximum for fluorinated fiber composites. Scanning electron microscopic analysis indicates that that the distribution of fiber into the matrix is homogeneous. It also indicates the better adhesion between the matrix and the reinforcing agent for modified fiber composites. We also did surface fluorination of the prepared composites and base polymer for knowing its application to different fields such as printability wettability, etc. To determine the various properties such as printability, wettability and adhesion properties, contact angle measurement was done. It was observed that the surface energies of surface modified composites and base polymer increases

Calculation of the electronic and magnetic structures of 3d impurities in the Hcp Fe matrix

International Nuclear Information System (INIS)

Franca, Fernando

1995-01-01

In this work we investigate the local magnetic properties and the electronic structure of HCP Fe, as well introducing transition metals atoms 3d (Cs, Ti, Cr, Mn, Co, Ni, Cu, Zn) in HCP iron matrix. We employed the discrete variational method (DVM), which is an orbital molecular method which incorporate the Hartree-Fock-Slater theory and the linear combination of atomic orbitals (LCAO), in the self-consistent charge approximation and the local density approximation of Von Barth and Hedin to the exchange-correlation potential. We used the embedded cluster model to investigate the electronic structure and the local magnetic properties for the central atom of a cluster of 27 atoms immersed in the microcrystal representing the HCP Fe. (author)

An approximate method for calculating electron-phonon matrix element of a disordered transition metal and relevant comments on superconductivity

International Nuclear Information System (INIS)

Zhang, L.

1981-08-01

A method based on the tight-binding approximation is developed to calculate the electron-phonon matrix element for the disordered transition metals. With the method as a basis the experimental Tsub(c) data of the amorphous transition metal superconductors are re-analysed. Some comments on the superconductivity of the disordered materials are given

Efficient k⋅p method for the calculation of total energy and electronic density of states

OpenAIRE

Iannuzzi, Marcella; Parrinello, Michele

2001-01-01

An efficient method for calculating the electronic structure in large systems with a fully converged BZ sampling is presented. The method is based on a k.p-like approximation developed in the framework of the density functional perturbation theory. The reliability and efficiency of the method are demostrated in test calculations on Ar and Si supercells

Correlated random-phase approximation from densities and in-medium matrix elements

Energy Technology Data Exchange (ETDEWEB)

Trippel, Richard; Roth, Robert [Institut fuer Kernphysik, Technische Universitaet Darmstadt (Germany)

2016-07-01

The random-phase approximation (RPA) as well as the second RPA (SRPA) are established tools for the study of collective excitations in nuclei. Addressing the well known lack of correlations, we derived a universal framework for a fully correlated RPA based on the use of one- and two-body densities. We apply densities from coupled cluster theory and investigate the impact of correlations. As an alternative approach to correlations we use matrix elements transformed via in-medium similarity renormalization group (IM-SRG) in combination with RPA and SRPA. We find that within SRPA the use of IM-SRG matrix elements leads to the disappearance of instabilities of low-lying states. For the calculations we use normal-ordered two- plus three-body interactions derived from chiral effective field theory. We apply different Hamiltonians to a number of doubly-magic nuclei and calculate electric transition strengths.

Electronic response and longitudinal phonons of a charge-density-wave distorted linear chain

International Nuclear Information System (INIS)

Giuliani, G.

1978-01-01

The longitudinal-phonon spectrum of an incommensurate charge-density-wave distorted linear chain at T = 0 K are calculated. This is done by direct numerical evaluation of the full static-electronic-response matrix. The electronic band structure assumed for this purpose is that of a mean-field theory 1-D Peierls insulator. The present results show how, within this simplified, but self-consistent picture, the phase and amplitude modes connect to, and interact with, the ordinary longitudinal-phonon branch. Effects due to our inclusion of (0,2ksub(F)) scattering along with the usual (-2ksub(F), 2ksub(F)) are also pointed out. An alternative approximate expression for the 1-D electronic-response matrix is also given. (author)

Orbital order and effective mass enhancement in t2 g two-dimensional electron gases

Science.gov (United States)

Tolsma, John; Principi, Alessandro; Polini, Marco; MacDonald, Allan

2015-03-01

It is now possible to prepare d-electron two-dimensional electron gas systems that are confined near oxide heterojunctions and contain t2 g electrons with a density much smaller than one electron per metal atom. I will discuss a generic model that captures all qualitative features of electron-electron interaction physics in t2 g two-dimensional electron gas systems, and the use of a GW approximation to explore t2 g quasiparticle properties in this new context. t2 g electron gases contain a high density isotropic light mass xy component and low-density xz and yz anisotropic components with light and heavy masses in orthogonal directions. The high density light mass band screens interactions within the heavy bands. As a result the wave vector dependence of the self-energy is reduced and the effective mass is increased. When the density in the heavy bands is low, the difference in anisotropy between the two heavy bands favors orbital order. When orbital order does not occur, interactions still reshape the heavy-band Fermi surfaces. I will discuss these results in the context of recently reported magnetotransport experiments.

Exact solution of a coupled spin–electron linear chain composed of localized Ising spins and mobile electrons

International Nuclear Information System (INIS)

Čisárová, Jana; Strečka, Jozef

2014-01-01

Exact solution of a coupled spin–electron linear chain composed of localized Ising spins and mobile electrons is found. The investigated spin–electron model is exactly solvable by the use of a transfer-matrix method after tracing out the degrees of freedom of mobile electrons delocalized over a couple of interstitial (decorating) sites. The exact ground-state phase diagram reveals an existence of five phases with different number of mobile electrons per unit cell, two of which are ferromagnetic, two are paramagnetic and one is antiferromagnetic. We have studied in particular the dependencies of compressibility and specific heat on temperature and electron density. - Highlights: • A coupled spin–electron chain composed of Ising spins and mobile electrons is exactly solved. • Quantum paramagnetic, ferromagnetic and antiferromagnetic ground states are found. • A compressibility shows a non-monotonous dependence on temperature and electron density. • Thermal dependences of specific heat display two distinct peaks

Annihilation probability density and other applications of the Schwinger multichannel method to the positron and electron scattering

International Nuclear Information System (INIS)

Varella, Marcio Teixeira do Nascimento

2001-12-01

We have calculated annihilation probability densities (APD) for positron collisions against He atom and H 2 molecule. It was found that direct annihilation prevails at low energies, while annihilation following virtual positronium (Ps) formation is the dominant mechanism at higher energies. In room-temperature collisions (10 -2 eV) the APD spread over a considerable extension, being quite similar to the electronic densities of the targets. The capture of the positron in an electronic Feshbach resonance strongly enhanced the annihilation rate in e + -H 2 collisions. We also discuss strategies to improve the calculation of the annihilation parameter (Z eff ), after debugging the computational codes of the Schwinger Multichannel Method (SMC). Finally, we consider the inclusion of the Ps formation channel in the SMC and show that effective configurations (pseudo eigenstates of the Hamiltonian of the collision ) are able to significantly reduce the computational effort in positron scattering calculations. Cross sections for electron scattering by polyatomic molecules were obtained in three different approximations: static-exchange (SE); tatic-exchange-plus-polarization (SEP); and multichannel coupling. The calculations for polar targets were improved through the rotational resolution of scattering amplitudes in which the SMC was combined with the first Born approximation (FBA). In general, elastic cross sections (SE and SEP approximations) showed good agreement with available experimental data for several targets. Multichannel calculations for e - -H 2 O scattering, on the other hand, presented spurious structures at the electronic excitation thresholds (author)

The finite element response matrix method

International Nuclear Information System (INIS)

Nakata, H.; Martin, W.R.

1983-02-01

A new technique is developed with an alternative formulation of the response matrix method implemented with the finite element scheme. Two types of response matrices are generated from the Galerkin solution to the weak form of the diffusion equation subject to an arbitrary current and source. The piecewise polynomials are defined in two levels, the first for the local (assembly) calculations and the second for the global (core) response matrix calculations. This finite element response matrix technique was tested in two 2-dimensional test problems, 2D-IAEA benchmark problem and Biblis benchmark problem, with satisfatory results. The computational time, whereas the current code is not extensively optimized, is of the same order of the well estabilished coarse mesh codes. Furthermore, the application of the finite element technique in an alternative formulation of response matrix method permits the method to easily incorporate additional capabilities such as treatment of spatially dependent cross-sections, arbitrary geometrical configurations, and high heterogeneous assemblies. (Author) [pt

Grand canonical electronic density-functional theory: Algorithms and applications to electrochemistry

International Nuclear Information System (INIS)

Sundararaman, Ravishankar; Goddard, William A. III; Arias, Tomas A.

2017-01-01

First-principles calculations combining density-functional theory and continuum solvation models enable realistic theoretical modeling and design of electrochemical systems. When a reaction proceeds in such systems, the number of electrons in the portion of the system treated quantum mechanically changes continuously, with a balancing charge appearing in the continuum electrolyte. A grand-canonical ensemble of electrons at a chemical potential set by the electrode potential is therefore the ideal description of such systems that directly mimics the experimental condition. We present two distinct algorithms: a self-consistent field method and a direct variational free energy minimization method using auxiliary Hamiltonians (GC-AuxH), to solve the Kohn-Sham equations of electronic density-functional theory directly in the grand canonical ensemble at fixed potential. Both methods substantially improve performance compared to a sequence of conventional fixed-number calculations targeting the desired potential, with the GC-AuxH method additionally exhibiting reliable and smooth exponential convergence of the grand free energy. Lastly, we apply grand-canonical density-functional theory to the under-potential deposition of copper on platinum from chloride-containing electrolytes and show that chloride desorption, not partial copper monolayer formation, is responsible for the second voltammetric peak.

Grand canonical electronic density-functional theory: Algorithms and applications to electrochemistry

Science.gov (United States)

Sundararaman, Ravishankar; Goddard, William A.; Arias, Tomas A.

2017-03-01

First-principles calculations combining density-functional theory and continuum solvation models enable realistic theoretical modeling and design of electrochemical systems. When a reaction proceeds in such systems, the number of electrons in the portion of the system treated quantum mechanically changes continuously, with a balancing charge appearing in the continuum electrolyte. A grand-canonical ensemble of electrons at a chemical potential set by the electrode potential is therefore the ideal description of such systems that directly mimics the experimental condition. We present two distinct algorithms: a self-consistent field method and a direct variational free energy minimization method using auxiliary Hamiltonians (GC-AuxH), to solve the Kohn-Sham equations of electronic density-functional theory directly in the grand canonical ensemble at fixed potential. Both methods substantially improve performance compared to a sequence of conventional fixed-number calculations targeting the desired potential, with the GC-AuxH method additionally exhibiting reliable and smooth exponential convergence of the grand free energy. Finally, we apply grand-canonical density-functional theory to the under-potential deposition of copper on platinum from chloride-containing electrolytes and show that chloride desorption, not partial copper monolayer formation, is responsible for the second voltammetric peak.

Analysis of the half-projected Hartree--Fock function: density matrix, natural orbitals, and configuration interaction equivalence

International Nuclear Information System (INIS)

Smeyers, Y.G.; Delgado-Barrio, G.

1976-01-01

The half-projected Hartree--Fock function for singlet states (HPHF) is analyzed in terms of natural electronic configurations. For this purpose the HPHF spinless density matrix and its natural orbitals are first deduced. It is found that the HPHF function does not contain any contribution from odd-times excited configurations. It is seen in addition, in the case of the singlet ground states, this function is approximately equivalent to two closed-shell configurations, although the nature of the excited one depends on the nuclear geometry. An example is given in the case of the LiH ground state. Finally, the application of this model for studying systems of more than two atoms is criticized

Numerical methods for strongly correlated many-body systems with bosonic degrees of freedom

International Nuclear Information System (INIS)

Dorfner, Florian Guenter

2017-01-01

Recent experimental advances allow the observation of electronic relaxation processes in solid-state systems in real time. After an initial excitation with an optical pulse, the relaxation depends on the microscopic interactions present in the system. The interaction of electrons with lattice degrees of freedom - the phonons - is ubiquitous in solids and, thus, it represents one of the most important relaxation channels. An analytic description of relaxation dynamics is hard to come by and very few exact solutions exist even for the equilibrium situation. Numerical methods are, in principle, able to solve the problem in both, equilibrium and out-of-equilibrium situations. However, wavefunction-based methods like exact diagonalization or the density matrix renormalization group method scale unfavorably in the number of local basis states. For electron-phonon coupled systems, the situation is especially severe because the local basis dimension can get very large depending on model parameters or in far-from-equilibrium situations. For groundstate problems, two independent strategies exist for density matrix renormalization group methods: the strictly single-site density matrix renormalization group method that scales linearly in the local dimension and the use of a local basis optimization scheme which truncates the local basis to a subset of the eigenstates of the local reduced density matrix with the largest eigenvalues - the optimal mode basis. In this thesis, we combine these two strategies in an improved algorithm which reduces the scaling from linear in the local dimension of the phonon occupation number basis to linear in the dimension of a smaller optimal mode basis. We demonstrate the improved scaling of this method on the example of the Holstein polaron and the half-filled Hubbard-Holstein model. We further describe an algorithm that combines the time-evolving block decimation method with a local basis optimization to lower the scaling with the local

Numerical methods for strongly correlated many-body systems with bosonic degrees of freedom

Energy Technology Data Exchange (ETDEWEB)

Dorfner, Florian Guenter

2017-02-23

Recent experimental advances allow the observation of electronic relaxation processes in solid-state systems in real time. After an initial excitation with an optical pulse, the relaxation depends on the microscopic interactions present in the system. The interaction of electrons with lattice degrees of freedom - the phonons - is ubiquitous in solids and, thus, it represents one of the most important relaxation channels. An analytic description of relaxation dynamics is hard to come by and very few exact solutions exist even for the equilibrium situation. Numerical methods are, in principle, able to solve the problem in both, equilibrium and out-of-equilibrium situations. However, wavefunction-based methods like exact diagonalization or the density matrix renormalization group method scale unfavorably in the number of local basis states. For electron-phonon coupled systems, the situation is especially severe because the local basis dimension can get very large depending on model parameters or in far-from-equilibrium situations. For groundstate problems, two independent strategies exist for density matrix renormalization group methods: the strictly single-site density matrix renormalization group method that scales linearly in the local dimension and the use of a local basis optimization scheme which truncates the local basis to a subset of the eigenstates of the local reduced density matrix with the largest eigenvalues - the optimal mode basis. In this thesis, we combine these two strategies in an improved algorithm which reduces the scaling from linear in the local dimension of the phonon occupation number basis to linear in the dimension of a smaller optimal mode basis. We demonstrate the improved scaling of this method on the example of the Holstein polaron and the half-filled Hubbard-Holstein model. We further describe an algorithm that combines the time-evolving block decimation method with a local basis optimization to lower the scaling with the local

Extending the precision and efficiency of the all-electron full-potential linearized augmented plane-wave density-functional theory method

International Nuclear Information System (INIS)

Michalicek, Gregor

2015-01-01

Density functional theory (DFT) is the most widely-used first-principles theory for analyzing, describing and predicting the properties of solids based on the fundamental laws of quantum mechanics. The success of the theory is a consequence of powerful approximations to the unknown exchange and correlation energy of the interacting electrons and of sophisticated electronic structure methods that enable the computation of the density functional equations on a computer. A widely used electronic structure method is the full-potential linearized augmented plane-wave (FLAPW) method, that is considered to be one of the most precise methods of its kind and often referred to as a standard. Challenged by the demand of treating chemically and structurally increasingly more complex solids, in this thesis this method is revisited and extended along two different directions: (i) precision and (ii) efficiency. In the full-potential linearized augmented plane-wave method the space of a solid is partitioned into nearly touching spheres, centered at each atom, and the remaining interstitial region between the spheres. The Kohn-Sham orbitals, which are used to construct the electron density, the essential quantity in DFT, are expanded into a linearized augmented plane-wave basis, which consists of plane waves in the interstitial region and angular momentum dependent radial functions in the spheres. In this thesis it is shown that for certain types of materials, e.g., materials with very broad electron bands or large band gaps, or materials that allow the usage of large space-filling spheres, the variational freedom of the basis in the spheres has to be extended in order to represent the Kohn-Sham orbitals with high precision over a large energy spread. Two kinds of additional radial functions confined to the spheres, so-called local orbitals, are evaluated and found to successfully eliminate this error. A new efficient basis set is developed, named linearized augmented lattice

Observation of electron plasma waves in plasma of two-temperature electrons

International Nuclear Information System (INIS)

Ikezawa, Shunjiro; Nakamura, Yoshiharu.

1981-01-01

Propagation of electron plasma waves in a large and unmagnetized plasma containing two Maxwellian distributions of electrons is studied experimentally. Two kinds of plasma sources which supply electrons of different temperature are used. The temperature ratio is about 3 and the density ratio of hot to cool electrons is varied from 0 to 0.5. A small contamination of hot electrons enhances the Landau damping of the principal mode known as the Bohm-Gross mode. When the density of hot electrons is larger than about 0.2, two modes are observed. The results agree with theoretical dispersion relations when excitation efficiencies of the modes are considered. (author)

Study of electron-molecule collisions via the finite-element method and R-matrix propagation technique: Model exchange

International Nuclear Information System (INIS)

Abdolsalami, F.; Abdolsalami, M.; Gomez, P.

1994-01-01

We have applied the finite-element method to electron-molecule collisions. All the calculations are done in the body frame within the fixed-nuclei approximation. A model potential, which is added to the static and polarization potential, has been used to represent the exchange effect. The method is applied to electron-H 2 scattering and the eigenphase sums and the cross sections obtained are in very good agreement with the corresponding results from the linear-algebraic approach. Finite-element calculations of the R matrix in the region where the static and exchange interactions are strong, however, has about one-half to one-fourth of the memory requirement of the linear-algebraic technique

Polypropylene/Polyaniline Nanofiber/Reduced Graphene Oxide Nanocomposite with Enhanced Electrical, Dielectric, and Ferroelectric Properties for a High Energy Density Capacitor.

Science.gov (United States)

Cho, Sunghun; Kim, Minkyu; Lee, Jun Seop; Jang, Jyongsik

2015-10-14

This work demonstrates a ternary nanocomposite system, composed of polypropylene (PP), redoped PANI (r-PANI) nanofibers, and reduced graphene oxides (RGOs), for use in a high energy density capacitor. r-PANI nanofibers were fabricated by the combination methods of chemical oxidation polymerization and secondary doping processes, resulting in higher conductivity (σ≈156 S cm(-1)) than that of the primarily doped PANI nanofibers (σ≈16 S cm(-1)). RGO sheets with high electron mobility and thermal stability can enhance the conductivity of r-PANI/RGO (σ≈220 S cm(-1)) and thermal stability of PP matrix. These findings could be extended to combine the advantages of r-PANI nanofibers and RGO sheets for developing an efficient means of preparing PP/r-PANI/RGO nanocomposite. When the r-PANI/RGO cofillers (10 vol %) were added to PP matrix, the resulting PP/r-PANI/RGO nanocomposite exhibited high dielectric constant (ε'≈51.8) with small dielectric loss (ε″≈9.3×10(-3)). Furthermore, the PP/r-PANI/RGO nanocomposite was used for an energy-harvesting device, which demonstrated high energy density (Ue≈12.6 J cm(-3)) and breakdown strength (E≈5.86×10(3) kV cm(-1)).

From Real Materials to Model Hamiltonians With Density Matrix Downfolding

Directory of Open Access Journals (Sweden)

Huihuo Zheng

2018-05-01

Full Text Available Due to advances in computer hardware and new algorithms, it is now possible to perform highly accurate many-body simulations of realistic materials with all their intrinsic complications. The success of these simulations leaves us with a conundrum: how do we extract useful physical models and insight from these simulations? In this article, we present a formal theory of downfolding–extracting an effective Hamiltonian from first-principles calculations. The theory maps the downfolding problem into fitting information derived from wave functions sampled from a low-energy subspace of the full Hilbert space. Since this fitting process most commonly uses reduced density matrices, we term it density matrix downfolding (DMD.

The Dirac operator on a finite domain and the R-matrix method

International Nuclear Information System (INIS)

Grant, I P

2008-01-01

Relativistic effects in electron-atom collisions and photo-excitation and -ionization processes increase in importance as the atomic number of the target atom grows and spin-dependent effects increase. A relativistic treatment in which electron motion is described using the Dirac Hamiltonian is then desirable. A version of the popular nonrelativistic R-matrix package incorporating terms from the Breit-Pauli Hamiltonian has been used for modelling such processes for some years. The fully relativistic Dirac R-matrix method has been less popular, but is becoming increasingly relevant for applications to heavy ion targets, where the need to use relativistic wavefunctions is more obvious. The Dirac R-matrix method has been controversial ever since it was first proposed by Goertzel (1948 Phys. Rev. 73 1463-6), and it is therefore important to confirm that recent elaborate and costly applications of the method, such as, Badnell et al (2004 J. Phys. B: At. Mol. Phys. 37 4589) and Ballance and Griffin (2007 J. Phys. B: At. Mol. Opt. Phys. 40 247-58), rest on secure foundations. The first part of this paper analyses the structure of the two-point boundary-value problem for the Dirac operator on a finite domain, from which we construct a unified derivation of the Schroedinger (nonrelativistic) and Dirac (relativistic) R-matrix methods. Suggestions that the usual relativistic theory is not well founded are shown to be without foundation

Microscopically based energy density functionals for nuclei using the density matrix expansion. II. Full optimization and validation

Science.gov (United States)

Navarro Pérez, R.; Schunck, N.; Dyhdalo, A.; Furnstahl, R. J.; Bogner, S. K.

2018-05-01

Background: Energy density functional methods provide a generic framework to compute properties of atomic nuclei starting from models of nuclear potentials and the rules of quantum mechanics. Until now, the overwhelming majority of functionals have been constructed either from empirical nuclear potentials such as the Skyrme or Gogny forces, or from systematic gradient-like expansions in the spirit of the density functional theory for atoms. Purpose: We seek to obtain a usable form of the nuclear energy density functional that is rooted in the modern theory of nuclear forces. We thus consider a functional obtained from the density matrix expansion of local nuclear potentials from chiral effective field theory. We propose a parametrization of this functional carefully calibrated and validated on selected ground-state properties that is suitable for large-scale calculations of nuclear properties. Methods: Our energy functional comprises two main components. The first component is a non-local functional of the density and corresponds to the direct part (Hartree term) of the expectation value of local chiral potentials on a Slater determinant. Contributions to the mean field and the energy of this term are computed by expanding the spatial, finite-range components of the chiral potential onto Gaussian functions. The second component is a local functional of the density and is obtained by applying the density matrix expansion to the exchange part (Fock term) of the expectation value of the local chiral potential. We apply the UNEDF2 optimization protocol to determine the coupling constants of this energy functional. Results: We obtain a set of microscopically constrained functionals for local chiral potentials from leading order up to next-to-next-to-leading order with and without three-body forces and contributions from Δ excitations. These functionals are validated on the calculation of nuclear and neutron matter, nuclear mass tables, single-particle shell structure

Symmetrized density matrix renormalization group algorithm for low-lying excited states of conjugated carbon systems: Application to 1,12-benzoperylene and polychrysene

Science.gov (United States)

Prodhan, Suryoday; Ramasesha, S.

2018-05-01

The symmetry adapted density matrix renormalization group (SDMRG) technique has been an efficient method for studying low-lying eigenstates in one- and quasi-one-dimensional electronic systems. However, the SDMRG method had bottlenecks involving the construction of linearly independent symmetry adapted basis states as the symmetry matrices in the DMRG basis were not sparse. We have developed a modified algorithm to overcome this bottleneck. The new method incorporates end-to-end interchange symmetry (C2) , electron-hole symmetry (J ) , and parity or spin-flip symmetry (P ) in these calculations. The one-to-one correspondence between direct-product basis states in the DMRG Hilbert space for these symmetry operations renders the symmetry matrices in the new basis with maximum sparseness, just one nonzero matrix element per row. Using methods similar to those employed in the exact diagonalization technique for Pariser-Parr-Pople (PPP) models, developed in the 1980s, it is possible to construct orthogonal SDMRG basis states while bypassing the slow step of the Gram-Schmidt orthonormalization procedure. The method together with the PPP model which incorporates long-range electronic correlations is employed to study the correlated excited-state spectra of 1,12-benzoperylene and a narrow mixed graphene nanoribbon with a chrysene molecule as the building unit, comprising both zigzag and cove-edge structures.

Simulation of the electron cloud density in BEPC II

International Nuclear Information System (INIS)

Liu Yudong; Guo Zhiyuan; Wang Jiuqing

2004-01-01

Electron Cloud Instability (ECI) may take place in positron storage ring when the machine is operated with multi-bunch positron beam. According to the actual shape of the vacuum chamber in the BEPC II, a program has been developed. With the code, authors can calculate the electron density in the chamber with different length of antechamber and the different secondary electron yield respectively. By the simulation, the possibility to put clearing electrodes in the chamber to reduce the electron density in the central region of the chamber is investigated. The simulation provides meaningful and important results for the BEPC II project and electron cloud instability research

Self-consistent embedding of density-matrix renormalization group wavefunctions in a density functional environment.

Science.gov (United States)

Dresselhaus, Thomas; Neugebauer, Johannes; Knecht, Stefan; Keller, Sebastian; Ma, Yingjin; Reiher, Markus

2015-01-28

We present the first implementation of a density matrix renormalization group algorithm embedded in an environment described by density functional theory. The frozen density embedding scheme is used with a freeze-and-thaw strategy for a self-consistent polarization of the orbital-optimized wavefunction and the environmental densities with respect to each other.

Response matrix method for large LMFBR analysis

International Nuclear Information System (INIS)

King, M.J.

1977-06-01

The feasibility of using response matrix techniques for computational models of large LMFBRs is examined. Since finite-difference methods based on diffusion theory have generally found a place in fast-reactor codes, a brief review of their general matrix foundation is given first in order to contrast it to the general strategy of response matrix methods. Then, in order to present the general method of response matrix technique, two illustrative examples are given. Matrix algorithms arising in the application to large LMFBRs are discussed, and the potential of the response matrix method is explored for a variety of computational problems. Principal properties of the matrices involved are derived with a view to application of numerical methods of solution. The Jacobi iterative method as applied to the current-balance eigenvalue problem is discussed

Relativistic density matrix in the diagonal momentum representation. Bose-gas

International Nuclear Information System (INIS)

Makhlin, A.N.; Sinyukov, Yu.M.

1984-01-01

The relativistic-invariance treatment of the ideal Bose-system arising from the diagonal momentum representation for the density matrix is developed. The average occupation members and their correlators for statistical systems in arbitrary inertial frames are found on the equal-time hypersurfaces. The relativistic partition function method for the calculation of thermodynamic properties of gases moving as a whole is constructed

A photoemission moments model using density functional and transfer matrix methods applied to coating layers on surfaces: Theory

Science.gov (United States)

Jensen, Kevin L.; Finkenstadt, Daniel; Shabaev, Andrew; Lambrakos, Samuel G.; Moody, Nathan A.; Petillo, John J.; Yamaguchi, Hisato; Liu, Fangze

2018-01-01

Recent experimental measurements of a bulk material covered with a small number of graphene layers reported by Yamaguchi et al. [NPJ 2D Mater. Appl. 1, 12 (2017)] (on bialkali) and Liu et al. [Appl. Phys. Lett. 110, 041607 (2017)] (on copper) and the needs of emission models in beam optics codes have lead to substantial changes in a Moments model of photoemission. The changes account for (i) a barrier profile and density of states factor based on density functional theory (DFT) evaluations, (ii) a Drude-Lorentz model of the optical constants and laser penetration depth, and (iii) a transmission probability evaluated by an Airy Transfer Matrix Approach. Importantly, the DFT results lead to a surface barrier profile of a shape similar to both resonant barriers and reflectionless wells: the associated quantum mechanical transmission probabilities are shown to be comparable to those recently required to enable the Moments (and Three Step) model to match experimental data but for reasons very different than the assumption by conventional wisdom that a barrier is responsible. The substantial modifications of the Moments model components, motivated by computational materials methods, are developed. The results prepare the Moments model for use in treating heterostructures and discrete energy level systems (e.g., quantum dots) proposed for decoupling the opposing metrics of performance that undermine the performance of advanced light sources like the x-ray Free Electron Laser. The consequences of the modified components on quantum yield, emittance, and emission models needed by beam optics codes are discussed.

Effective one-body potential of DFT plus correlated kinetic energy density for two-electron spherical model atoms

International Nuclear Information System (INIS)

March, N.H.; Ludena, Eduardo V.

2004-01-01

For three model problems concerning two-electron spin-compensated ground states with spherical density, the third-order linear homogeneous differential equation constructed for the determination of ρ(r) is used here in conjunction with the von Weizsacker functional to characterize the one-body potential of density functional theory (DFT). Correlated von Weizsacker-type terms are compared to the exact DFT functional

Electron density profile measurements from hydrogen line intensity ratio method in Versatile Experimental Spherical Torus

Energy Technology Data Exchange (ETDEWEB)

Kim, YooSung; Shi, Yue-Jiang, E-mail: yjshi@snu.ac.kr; Yang, Jeong-hun; Kim, SeongCheol; Kim, Young-Gi; Dang, Jeong-Jeung; Yang, Seongmoo; Jo, Jungmin; Chung, Kyoung-Jae [Department of Nuclear Engineering, Seoul National University, Seoul 151-744 (Korea, Republic of); Oh, Soo-Ghee [Division of Energy Systems Research, Ajou University, Suwon 442-749 (Korea, Republic of); Hwang, Y. S. [Department of Nuclear Engineering, Seoul National University, Seoul 151-744 (Korea, Republic of); Center for Advanced Research in Fusion Reactor Engineering, Seoul National University, Seoul 151-744 (Korea, Republic of)

2016-11-15

Electron density profiles of versatile experiment spherical torus plasmas are measured by using a hydrogen line intensity ratio method. A fast-frame visible camera with appropriate bandpass filters is used to detect images of Balmer line intensities. The unique optical system makes it possible to take images of H{sub α} and H{sub β} radiation simultaneously, with only one camera. The frame rate is 1000 fps and the spatial resolution of the system is about 0.5 cm. One-dimensional local emissivity profiles have been obtained from the toroidal line of sight with viewing dumps. An initial result for the electron density profile is presented and is in reasonable agreement with values measured by a triple Langmuir probe.

Scanning laser densitometry and color perimetry demonstrate reduced photopigment density and sensitivity in two patients with retinal degeneration.

Science.gov (United States)

Tornow, R P; Stilling, R; Zrenner, E

1999-10-01

To test the feasibility of scanning laser densitometry with a modified Rodenstock scanning laser ophthalmoscope (SLO) to measure the rod and cone photopigment distribution in patients with retinal diseases. Scanning laser densitometry was performed using a modified Rodenstock scanning laser ophthalmoscope. The distribution of the photopigments was calculated from dark adapted and bleached images taken with the 514 nm laser of the SLO. This wavelength is absorbed by rod and cone photopigments. Discrimination is possible due to their different spatial distribution. Additionally, to measure retinal sensitivity profiles, dark adapted two color static perimetry with a Tübinger manual perimeter was performed along the horizontal meridian with 1 degree spacing. A patient with retinitis pigmentosa had slightly reduced photopigment density within the central +/- 5 degrees but no detectable photopigment for eccentricities beyond 5 degrees. A patient with cone dystrophy had nearly normal pigment density beyond +/- 5 degrees, but considerably reduced photopigment density within the central +/- 5 degrees. Within the central +/- 5 degrees, the patient with retinitis pigmentosa had normal sensitivity for the red stimulus and reduced sensitivity for the green stimulus. There was no measurable function beyond 7 degrees. The patient with cone dystrophy had normal sensitivity for the green stimulus outside the foveal center and reduced sensitivity for the red stimulus at the foveal center. The results of color perimetry for this patient with a central scotoma were probably influenced by eccentric fixation. Scanning laser densitometry with a modified Rodenstock SLO is a useful method to assess the human photopigment distribution. Densitometry results were confirmed by dark adapted two color static perimetry. Photopigment distribution and retinal sensitivity profiles can be measured with high spatial resolution. This may help to measure exactly the temporal development of retinal

Analysis of Nonlinear Dynamics by Square Matrix Method

Energy Technology Data Exchange (ETDEWEB)

Yu, Li Hua [Brookhaven National Lab. (BNL), Upton, NY (United States). Energy and Photon Sciences Directorate. National Synchrotron Light Source II

2016-07-25

The nonlinear dynamics of a system with periodic structure can be analyzed using a square matrix. In this paper, we show that because the special property of the square matrix constructed for nonlinear dynamics, we can reduce the dimension of the matrix from the original large number for high order calculation to low dimension in the first step of the analysis. Then a stable Jordan decomposition is obtained with much lower dimension. The transformation to Jordan form provides an excellent action-angle approximation to the solution of the nonlinear dynamics, in good agreement with trajectories and tune obtained from tracking. And more importantly, the deviation from constancy of the new action-angle variable provides a measure of the stability of the phase space trajectories and their tunes. Thus the square matrix provides a novel method to optimize the nonlinear dynamic system. The method is illustrated by many examples of comparison between theory and numerical simulation. Finally, in particular, we show that the square matrix method can be used for optimization to reduce the nonlinearity of a system.

Precision Electron Density Measurements in the SSX MHD Wind Tunnel

Science.gov (United States)

Suen-Lewis, Emma M.; Barbano, Luke J.; Shrock, Jaron E.; Kaur, Manjit; Schaffner, David A.; Brown, Michael R.

2017-10-01

We characterize fluctuations of the line averaged electron density of Taylor states produced by the magnetized coaxial plasma gun of the SSX device using a 632.8 nm HeNe laser interferometer. The analysis method uses the electron density dependence of the refractive index of the plasma to determine the electron density of the Taylor states. Typical magnetic field and density values in the SSX device approach about B ≅ 0.3 T and n = 0 . 4 ×1016 cm-3 . Analysis is improved from previous density measurement methods by developing a post-processing method to remove relative phase error between interferometer outputs and to account for approximately linear phase drift due to low-frequency mechanical vibrations of the interferometer. Precision density measurements coupled with local measurements of the magnetic field will allow us to characterize the wave composition of SSX plasma via density vs. magnetic field correlation analysis, and compare the wave composition of SSX plasma with that of the solar wind. Preliminary results indicate that density and magnetic field appear negatively correlated. Work supported by DOE ARPA-E ALPHA program.

Five- and six-electron harmonium atoms: Highly accurate electronic properties and their application to benchmarking of approximate 1-matrix functionals

Science.gov (United States)

Cioslowski, Jerzy; Strasburger, Krzysztof

2018-04-01

Electronic properties of several states of the five- and six-electron harmonium atoms are obtained from large-scale calculations employing explicitly correlated basis functions. The high accuracy of the computed energies (including their components), natural spinorbitals, and their occupation numbers makes them suitable for testing, calibration, and benchmarking of approximate formalisms of quantum chemistry and solid state physics. In the case of the five-electron species, the availability of the new data for a wide range of the confinement strengths ω allows for confirmation and generalization of the previously reached conclusions concerning the performance of the presently known approximations for the electron-electron repulsion energy in terms of the 1-matrix that are at heart of the density matrix functional theory (DMFT). On the other hand, the properties of the three low-lying states of the six-electron harmonium atom, computed at ω = 500 and ω = 1000, uncover deficiencies of the 1-matrix functionals not revealed by previous studies. In general, the previously published assessment of the present implementations of DMFT being of poor accuracy is found to hold. Extending the present work to harmonically confined systems with even more electrons is most likely counterproductive as the steep increase in computational cost required to maintain sufficient accuracy of the calculated properties is not expected to be matched by the benefits of additional information gathered from the resulting benchmarks.

Time dependent density matrix theory and effective interaction

Energy Technology Data Exchange (ETDEWEB)

Tohyama, Mitsuru [Kyorin Univ., Mitaka, Tokyo (Japan). School of Medicine

1998-07-01

A correlated ground state of {sup 16}O and an E2 giant resonance built on it are calculated using an extended version of the time-dependent Hartree-Fock theory called the time-dependent density-matrix theory (TDDM). The Skyrme force is used in the calculation of both a mean field and two-body correlations. It is found that TDDM gives reasonable ground-state correlations and a large spreading width of the E2 giant resonance when single-particle states in the continuum are treated appropriately. (author)

Fast-electron self-collimation in a plasma density gradient

International Nuclear Information System (INIS)

Yang, X. H.; Borghesi, M.; Robinson, A. P. L.

2012-01-01

A theoretical and numerical study of fast electron transport in solid and compressed fast ignition relevant targets is presented. The principal aim of the study is to assess how localized increases in the target density (e.g., by engineering of the density profile) can enhance magnetic field generation and thus pinching of the fast electron beam through reducing the rate of temperature rise. The extent to which this might benefit fast ignition is discussed.

Density functional application to strongly correlated electron systems

International Nuclear Information System (INIS)

Eschrig, H.; Koepernik, K.; Chaplygin, I.

2003-01-01

The local spin density approximation plus onsite Coulomb repulsion approach (LSDA+U) to density functional theory is carefully reanalyzed. Its possible link to single-particle Green's function theory is occasionally discussed. A simple and elegant derivation of the important sum rules for the on-site interaction matrix elements linking them to the values of U and J is presented. All necessary expressions for an implementation of LSDA+U into a non-orthogonal basis solver for the Kohn-Sham equations are given, and implementation into the full-potential local-orbital solver (Phys. Rev. B 59 (1999) 1743) is made. Results of application to several planar cuprate structures are reported in detail and conclusions on the interpretation of the physics of the electronic structure of the cuprates are drawn

Sparse Matrix for ECG Identification with Two-Lead Features

Directory of Open Access Journals (Sweden)

Kuo-Kun Tseng

2015-01-01

Full Text Available Electrocardiograph (ECG human identification has the potential to improve biometric security. However, improvements in ECG identification and feature extraction are required. Previous work has focused on single lead ECG signals. Our work proposes a new algorithm for human identification by mapping two-lead ECG signals onto a two-dimensional matrix then employing a sparse matrix method to process the matrix. And that is the first application of sparse matrix techniques for ECG identification. Moreover, the results of our experiments demonstrate the benefits of our approach over existing methods.

The finite element response Matrix method

International Nuclear Information System (INIS)

Nakata, H.; Martin, W.R.

1983-01-01

A new method for global reactor core calculations is described. This method is based on a unique formulation of the response matrix method, implemented with a higher order finite element method. The unique aspects of this approach are twofold. First, there are two levels to the overall calculational scheme: the local or assembly level and the global or core level. Second, the response matrix scheme, which is formulated at both levels, consists of two separate response matrices rather than one response matrix as is generally the case. These separate response matrices are seen to be quite beneficial for the criticality eigenvalue calculation, because they are independent of k /SUB eff/. The response matrices are generated from a Galerkin finite element solution to the weak form of the diffusion equation, subject to an arbitrary incoming current and an arbitrary distributed source. Calculational results are reported for two test problems, the two-dimensional International Atomic Energy Agency benchmark problem and a two-dimensional pressurized water reactor test problem (Biblis reactor), and they compare well with standard coarse mesh methods with respect to accuracy and efficiency. Moreover, the accuracy (and capability) is comparable to fine mesh for a fraction of the computational cost. Extension of the method to treat heterogeneous assemblies and spatial depletion effects is discussed

Numerical Methods for a Multicomponent Two-Phase Interface Model with Geometric Mean Influence Parameters

KAUST Repository

Kou, Jisheng

2015-07-16

In this paper, we consider an interface model for multicomponent two-phase fluids with geometric mean influence parameters, which is popularly used to model and predict surface tension in practical applications. For this model, there are two major challenges in theoretical analysis and numerical simulation: the first one is that the influence parameter matrix is not positive definite; the second one is the complicated structure of the energy function, which requires us to find out a physically consistent treatment. To overcome these two challenging problems, we reduce the formulation of the energy function by employing a linear transformation and a weighted molar density, and furthermore, we propose a local minimum grand potential energy condition to establish the relation between the weighted molar density and mixture compositions. From this, we prove the existence of the solution under proper conditions and prove the maximum principle of the weighted molar density. For numerical simulation, we propose a modified Newton\\'s method for solving this nonlinear model and analyze its properties; we also analyze a finite element method with a physical-based adaptive mesh-refinement technique. Numerical examples are tested to verify the theoretical results and the efficiency of the proposed methods.

Numerical Methods for a Multicomponent Two-Phase Interface Model with Geometric Mean Influence Parameters

KAUST Repository

Kou, Jisheng; Sun, Shuyu

2015-01-01

In this paper, we consider an interface model for multicomponent two-phase fluids with geometric mean influence parameters, which is popularly used to model and predict surface tension in practical applications. For this model, there are two major challenges in theoretical analysis and numerical simulation: the first one is that the influence parameter matrix is not positive definite; the second one is the complicated structure of the energy function, which requires us to find out a physically consistent treatment. To overcome these two challenging problems, we reduce the formulation of the energy function by employing a linear transformation and a weighted molar density, and furthermore, we propose a local minimum grand potential energy condition to establish the relation between the weighted molar density and mixture compositions. From this, we prove the existence of the solution under proper conditions and prove the maximum principle of the weighted molar density. For numerical simulation, we propose a modified Newton's method for solving this nonlinear model and analyze its properties; we also analyze a finite element method with a physical-based adaptive mesh-refinement technique. Numerical examples are tested to verify the theoretical results and the efficiency of the proposed methods.

On the statistical interpretation of quantum mechanics: evolution of the density matrix

International Nuclear Information System (INIS)

Benzecri, J.-P.

1986-01-01

Using two classical examples (the Young slit experiment and coherent and incoherent crystal diffraction of neutrons) we show in a general framework, that for a system viewed as consisting of two components, depolarisation of the density matrix by one of these can result from the application of the Schroedinger equation to the global system [fr

Density functional theory of electron transfer beyond the Born-Oppenheimer approximation: Case study of LiF

Science.gov (United States)

Li, Chen; Requist, Ryan; Gross, E. K. U.

2018-02-01

We perform model calculations for a stretched LiF molecule, demonstrating that nonadiabatic charge transfer effects can be accurately and seamlessly described within a density functional framework. In alkali halides like LiF, there is an abrupt change in the ground state electronic distribution due to an electron transfer at a critical bond length R = Rc, where an avoided crossing of the lowest adiabatic potential energy surfaces calls the validity of the Born-Oppenheimer approximation into doubt. Modeling the R-dependent electronic structure of LiF within a two-site Hubbard model, we find that nonadiabatic electron-nuclear coupling produces a sizable elongation of the critical Rc by 0.5 bohr. This effect is very accurately captured by a simple and rigorously derived correction, with an M-1 prefactor, to the exchange-correlation potential in density functional theory, M = reduced nuclear mass. Since this nonadiabatic term depends on gradients of the nuclear wave function and conditional electronic density, ∇Rχ(R) and ∇Rn(r, R), it couples the Kohn-Sham equations at neighboring R points. Motivated by an observed localization of nonadiabatic effects in nuclear configuration space, we propose a local conditional density approximation—an approximation that reduces the search for nonadiabatic density functionals to the search for a single function y(n).

Density-matrix renormalization group method for the conductance of one-dimensional correlated systems using the Kubo formula

Science.gov (United States)

Bischoff, Jan-Moritz; Jeckelmann, Eric

2017-11-01

We improve the density-matrix renormalization group (DMRG) evaluation of the Kubo formula for the zero-temperature linear conductance of one-dimensional correlated systems. The dynamical DMRG is used to compute the linear response of a finite system to an applied ac source-drain voltage; then the low-frequency finite-system response is extrapolated to the thermodynamic limit to obtain the dc conductance of an infinite system. The method is demonstrated on the one-dimensional spinless fermion model at half filling. Our method is able to replicate several predictions of the Luttinger liquid theory such as the renormalization of the conductance in a homogeneous conductor, the universal effects of a single barrier, and the resonant tunneling through a double barrier.

Matrix elements of intraband transitions in quantum dot intermediate band solar cells: the influence of quantum dot presence on the extended-state electron wave-functions

International Nuclear Information System (INIS)

Nozawa, Tomohiro; Arakawa, Yasuhiko

2014-01-01

The intraband transitions which are essential for quantum dot intermediate band solar cells (QD IBSCs) are theoretically investigated by estimating the matrix elements from a ground bound state, which is often regarded as an intermediate band (IB), to conduction band (CB) states for a structure with a quantum dot (QD) embedded in a matrix (a QD/matrix structure). We have found that the QD pushes away the electron envelope functions (probability densities) from the QD region in almost all quantum states above the matrix CB minimum. As a result, the matrix elements of the intraband transitions in the QD/matrix structure are largely reduced, compared to those calculated assuming the envelope functions of free electrons (i.e., plane-wave envelope functions) in a matrix structure as the final states of the intraband transitions. The result indicates the strong influence of the QD itself on the intraband transitions from the IB to the CB states in QD IBSC devices. This work will help in better understanding the problem of the intraband transitions and give new insight, that is, engineering of quantum states is indispensable for the realization of QD IBSCs with high solar energy conversion efficiencies. (paper)

Probing the Milky Way electron density using multi-messenger astronomy

Science.gov (United States)

Breivik, Katelyn; Larson, Shane

2015-04-01

Multi-messenger observations of ultra-compact binaries in both gravitational waves and electromagnetic radiation supply highly complementary information, providing new ways of characterizing the internal dynamics of these systems, as well as new probes of the galaxy itself. Electron density models, used in pulsar distance measurements via the electron dispersion measure, are currently not well constrained. Simultaneous radio and gravitational wave observations of pulsars in binaries provide a method of measuring the average electron density along the line of sight to the pulsar, thus giving a new method for constraining current electron density models. We present this method and assess its viability with simulations of the compact binary component of the Milky Way using the public domain binary evolution code, BSE. This work is supported by NASA Award NNX13AM10G.

Links between matrix bulk density, macropore characteristics and hydraulic behavior of soils

DEFF Research Database (Denmark)

Katuwal, Sheela; Møldrup, Per; Lamandé, Mathieu

2013-01-01

characteristics on soil hydraulic functions has rarely been studied. With the objective of studying the links between these parameters we quantified macropore characteristics of intact soil columns (19 cm diameter x 20 cm high) from two agricultural field sites (Silstrup and Faardrup) in Denmark using coarse...... resolution X-ray CT and linked them with laboratory measurements of air permeability and leaching experiment. In addition to macropore characteristics, we also quantified the CT-number of the matrix as a measure of the bulk density of the matrix, i.e., excluding macropores in the soil. Soils from the two...

Stability of Sarma phases in density imbalanced electron-hole bilayer systems

International Nuclear Information System (INIS)

Subasi, A. L.; Tanatar, B.; Pieri, P.; Senatore, G.

2010-01-01

We study excitonic condensation in an electron-hole bilayer system with unequal layer densities at zero temperature. Using mean-field theory we solve the Bardeen-Cooper-Schrieffer (BCS) gap equations numerically and investigate the effects of intralayer interactions. The electron-hole system evolves from BCS in the weak coupling limit to Bose-Einstein condensation (BEC) in the strong coupling limit. We analyze the stability of the Sarma phase with k,-k pairing by calculating the superfluid mass density and also by checking the compressibility matrix. We find that with bare Coulomb interactions the superfluid density is always positive in the Sarma phase, due to a peculiar momentum structure of the gap function originating from the singular behavior of the Coulomb potential at zero momentum and the presence of a sharp Fermi surface. Introducing a simple model for screening, we find that the superfluid density becomes negative in some regions of the phase diagram, corresponding to an instability toward a Fulde-Ferrel-Larkin-Ovchinnikov-type superfluid phase. Thus, intralayer interaction and screening together can lead to a rich phase diagram in the BCS-BEC crossover regime in electron-hole bilayer systems.

Uncertainty relations and reduced density matrices: Mapping many-body quantum mechanics onto four particles

Science.gov (United States)

Mazziotti, David A.; Erdahl, Robert M.

2001-04-01

For the description of ground-state correlation phenomena an accurate mapping of many-body quantum mechanics onto four particles is developed. The energy for a quantum system with no more than two-particle interactions may be expressed in terms of a two-particle reduced density matrix (2-RDM), but variational optimization of the 2-RDM requires that it corresponds to an N-particle wave function. We derive N-representability conditions on the 2-RDM that guarantee the validity of the uncertainty relations for all operators with two-particle interactions. One of these conditions is shown to be necessary and sufficient to make the RDM solutions of the dispersion condition equivalent to those from the contracted Schrödinger equation (CSE) [Mazziotti, Phys. Rev. A 57, 4219 (1998)]. In general, the CSE is a stronger N-representability condition than the dispersion condition because the CSE implies the dispersion condition as well as additional N-representability constraints from the Hellmann-Feynman theorem. Energy minimization subject to the representability constraints is performed for a boson model with 10, 30, and 75 particles. Even when traditional wave-function methods fail at large perturbations, the present method yields correlation energies within 2%.

Deriving the coronal hole electron temperature: electron density dependent ionization / recombination considerations

International Nuclear Information System (INIS)

Doyle, John Gerard; Perez-Suarez, David; Singh, Avninda; Chapman, Steven; Bryans, Paul; Summers, Hugh; Savin, Daniel Wolf

2010-01-01

Comparison of appropriate theoretically derived line ratios with observational data can yield estimates of a plasma's physical parameters, such as electron density or temperature. The usual practice in the calculation of the line ratio is the assumption of excitation by electrons/protons followed by radiative decay. Furthermore, it is normal to use the so-called coronal approximation, i.e. one only considers ionization and recombination to and from the ground-state. A more accurate treatment is to include ionization/recombination to and from metastable levels. Here, we apply this to two lines from adjacent ionization stages, Mg IX 368 A and Mg X 625 A, which has been shown to be a very useful temperature diagnostic. At densities typical of coronal hole conditions, the difference between the electron temperature derived assuming the zero density limit compared with the electron density dependent ionization/recombination is small. This, however, is not the case for flares where the electron density is orders of magnitude larger. The derived temperature for the coronal hole at solar maximum is around 1.04 MK compared to just below 0.82 MK at solar minimum.

Assembling phosphorene flexagons for 2D electron-density-guided nanopatterning and nanofabrication.

Science.gov (United States)

Kang, Kisung; Jang, Woosun; Soon, Aloysius

2017-07-27

To build upon the rich structural diversity in the ever-increasing polymorphic phases of two-dimensional phosphorene, we propose different assembly methods (namely, the "bottom-up" and "top-down" approaches) that involve four commonly reported parent phases (i.e. the α-, β-, γ-, and δ-phosphorene) in combination with the lately reported remarkably low-energy one-dimensional defects in α-phosphorene. In doing so, we generate various periodically repeated phosphorene patterns in these so-called phosphorene flexagons and present their local electron density (via simulated scanning tunneling microscopy (STM) images). These interesting electron density patterns seen in the flexagons (mimicking symmetry patterns that one may typically see in a kaleidoscope) may assist as potential 2D templates where electron-density-guided nanopatterning and nanofabrication in complex organized nanoarchitectures are important.

Density effects on electronic configurations in dense plasmas

Science.gov (United States)

Faussurier, Gérald; Blancard, Christophe

2018-02-01

We present a quantum mechanical model to describe the density effects on electronic configurations inside a plasma environment. Two different approaches are given by starting from a quantum average-atom model. Illustrations are shown for an aluminum plasma in local thermodynamic equilibrium at solid density and at a temperature of 100 eV and in the thermodynamic conditions of a recent experiment designed to characterize the effects of the ionization potential depression treatment. Our approach compares well with experiment and is consistent in that case with the approach of Stewart and Pyatt to describe the ionization potential depression rather than with the method of Ecker and Kröll.

NATO Advanced Study Institute on Electron and Magnetization Densities in Molecules and Crystals

CERN Document Server

1980-01-01

The interest of describing the ground state properties of a system in terms of one electron density (or its two spin components) is obvious, in particular due to the simple physical significance of this function. Recent experimental progress in diffraction made the measurement of charge and magnetization densities in crystalline solids possible, with an accuracy at least as good as theoretical accuracy. Theoretical developments of the many-body problem have proved the extreme importance of the one electron density function and presently, accurate methods of band structure determination become available. Parallel to the diffraction techniques, other domains of research (inelastic scattering, resonance, molecular spectroscopy) deal with quantities directly related to the one particle density. But the two types of studies do not interfere enough and one should obviously gain more information by interpreting all experiments that are related to the density together. It became necessary to have an International Sch...

Assessing the effect of electron density in photon dose calculations

International Nuclear Information System (INIS)

Seco, J.; Evans, P. M.

2006-01-01

Photon dose calculation algorithms (such as the pencil beam and collapsed cone, CC) model the attenuation of a primary photon beam in media other than water, by using pathlength scaling based on the relative mass density of the media to water. In this study, we assess if differences in the electron density between the water and media, with different atomic composition, can influence the accuracy of conventional photon dose calculations algorithms. A comparison is performed between an electron-density scaling method and the standard mass-density scaling method for (i) tissues present in the human body (such as bone, muscle, etc.), and for (ii) water-equivalent plastics, used in radiotherapy dosimetry and quality assurance. We demonstrate that the important material property that should be taken into account by photon dose algorithms is the electron density, and not the mass density. The mass-density scaling method is shown to overestimate, relative to electron-density predictions, the primary photon fluence for tissues in the human body and water-equivalent plastics, where 6%-7% and 10% differences were observed respectively for bone and air. However, in the case of patients, differences are expected to be smaller due to the large complexity of a treatment plan and of the patient anatomy and atomic composition and of the smaller thickness of bone/air that incident photon beams of a treatment plan may have to traverse. Differences have also been observed for conventional dose algorithms, such as CC, where an overestimate of the lung dose occurs, when irradiating lung tumors. The incorrect lung dose can be attributed to the incorrect modeling of the photon beam attenuation through the rib cage (thickness of 2-3 cm in bone upstream of the lung tumor) and through the lung and the oversimplified modeling of electron transport in convolution algorithms. In the present study, the overestimation of the primary photon fluence, using the mass-density scaling method, was shown

A Golub-Kahan-type reduction method for matrix pairs

NARCIS (Netherlands)

Hochstenbach, M.E.; Reichel, L.; Yu, X.

2015-01-01

We describe a novel method for reducing a pair of large matrices {A;B} to a pair of small matrices {H;K}. The method is an extension of Golub-Kahan bidiagonalization to matrix pairs, and simplifies to the latter method when B is the identity matrix. Applications to Tikhonov regularization of large

A Golub-Kahan-type reduction method for matrix pairs

NARCIS (Netherlands)

Hochstenbach, M.E.; Reichel, L.; Yu, X.

2015-01-01

We describe a novel method for reducing a pair of large matrices {A,B} to a pair of small matrices {H,K}. The method is an extension of Golub–Kahan bidiagonalization to matrix pairs, and simplifies to the latter method when B is the identity matrix. Applications to Tikhonov regularization of large

Quantitative evaluation of the matrix effect in bioanalytical methods based on LC-MS: A comparison of two approaches.

Science.gov (United States)

Rudzki, Piotr J; Gniazdowska, Elżbieta; Buś-Kwaśnik, Katarzyna

2018-06-05

Liquid chromatography coupled to mass spectrometry (LC-MS) is a powerful tool for studying pharmacokinetics and toxicokinetics. Reliable bioanalysis requires the characterization of the matrix effect, i.e. influence of the endogenous or exogenous compounds on the analyte signal intensity. We have compared two methods for the quantitation of matrix effect. The CVs(%) of internal standard normalized matrix factors recommended by the European Medicines Agency were evaluated against internal standard normalized relative matrix effects derived from Matuszewski et al. (2003). Both methods use post-extraction spiked samples, but matrix factors require also neat solutions. We have tested both approaches using analytes of diverse chemical structures. The study did not reveal relevant differences in the results obtained with both calculation methods. After normalization with the internal standard, the CV(%) of the matrix factor was on average 0.5% higher than the corresponding relative matrix effect. The method adopted by the European Medicines Agency seems to be slightly more conservative in the analyzed datasets. Nine analytes of different structures enabled a general overview of the problem, still, further studies are encouraged to confirm our observations. Copyright © 2018 Elsevier B.V. All rights reserved.

A method to measure the suprathermal density distribution by electron cyclotron emission

International Nuclear Information System (INIS)

Tutter, M.

1986-05-01

Electron cyclotron emission spectra of suprathermal electrons in a thermal main plasma are calculated. It is shown that for direction of observation oblique to the magnetic field, which decays in direction to the receiver, one may obtain information on the spatial density distribution of the suprathermal electrons from those spectra. (orig.)

Large-N limit of the two-Hermitian-matrix model by the hidden BRST method

International Nuclear Information System (INIS)

Alfaro, J.

1993-01-01

This paper discusses the large-N limit of the two-Hermitian-matrix model in zero dimensions, using the hidden Becchi-Rouet-Stora-Tyutin method. A system of integral equations previously found is solved, showing that it contained the exact solution of the model in leading order of large N

Sparse maps—A systematic infrastructure for reduced-scaling electronic structure methods. I. An efficient and simple linear scaling local MP2 method that uses an intermediate basis of pair natural orbitals

Energy Technology Data Exchange (ETDEWEB)

Pinski, Peter; Riplinger, Christoph; Neese, Frank, E-mail: evaleev@vt.edu, E-mail: frank.neese@cec.mpg.de [Max Planck Institute for Chemical Energy Conversion, Stiftstr. 34-36, D-45470 Mülheim an der Ruhr (Germany); Valeev, Edward F., E-mail: evaleev@vt.edu, E-mail: frank.neese@cec.mpg.de [Department of Chemistry, Virginia Tech, Blacksburg, Virginia 24061 (United States)

2015-07-21

In this work, a systematic infrastructure is described that formalizes concepts implicit in previous work and greatly simplifies computer implementation of reduced-scaling electronic structure methods. The key concept is sparse representation of tensors using chains of sparse maps between two index sets. Sparse map representation can be viewed as a generalization of compressed sparse row, a common representation of a sparse matrix, to tensor data. By combining few elementary operations on sparse maps (inversion, chaining, intersection, etc.), complex algorithms can be developed, illustrated here by a linear-scaling transformation of three-center Coulomb integrals based on our compact code library that implements sparse maps and operations on them. The sparsity of the three-center integrals arises from spatial locality of the basis functions and domain density fitting approximation. A novel feature of our approach is the use of differential overlap integrals computed in linear-scaling fashion for screening products of basis functions. Finally, a robust linear scaling domain based local pair natural orbital second-order Möller-Plesset (DLPNO-MP2) method is described based on the sparse map infrastructure that only depends on a minimal number of cutoff parameters that can be systematically tightened to approach 100% of the canonical MP2 correlation energy. With default truncation thresholds, DLPNO-MP2 recovers more than 99.9% of the canonical resolution of the identity MP2 (RI-MP2) energy while still showing a very early crossover with respect to the computational effort. Based on extensive benchmark calculations, relative energies are reproduced with an error of typically <0.2 kcal/mol. The efficiency of the local MP2 (LMP2) method can be drastically improved by carrying out the LMP2 iterations in a basis of pair natural orbitals. While the present work focuses on local electron correlation, it is of much broader applicability to computation with sparse tensors in

Preionization electron density measurement by collecting electric charge

International Nuclear Information System (INIS)

Giordano, G.; Letardi, T.

1988-01-01

A method using electron collection for preionization-electron number density measurements is presented. A cathode-potential drop model is used to describe the measurement principle. There is good agreement between the model and the experimental result

Multiple resonance compensation for betatron coupling and its equivalence with matrix method

CERN Document Server

De Ninno, G

1999-01-01

Analyses of betatron coupling can be broadly divided into two categories: the matrix approach that decouples the single-turn matrix to reveal the normal modes and the hamiltonian approach that evaluates the coupling in terms of the action of resonances in perturbation theory. The latter is often regarded as being less exact but good for physical insight. The common opinion is that the correction of the two closest sum and difference resonances to the working point is sufficient to reduce the off-axis terms in the 4X4 single-turn matrix, but this is only partially true. The reason for this is explained, and a method is developed that sums to infinity all coupling resonances and, in this way, obtains results equivalent to the matrix approach. The two approaches is discussed with reference to the dynamic aperture. Finally, the extension of the summation method to resonances of all orders is outlined and the relative importance of a single resonance compared to all resonances of a given order is analytically desc...

The Stark-crossing method for the simultaneous determination of the electron temperature and density in plasmas

International Nuclear Information System (INIS)

Torres, J; Carabano, O; Fernandez, M; Rubio, S; Alvarez, R; Rodero, A; Lao, C; Quintero, M C; Gamero, A; Sola, A

2006-01-01

The use of the Stark broadening of Balmer lines spontaneously emitted by atmospheric-pressure plasmas as a method to determine both the electron density and temperature in high-pressure plasmas is discussed in this paper. This method is applied to argon and helium plasmas produced in microwave discharges. Especially for Ar plasmas, valuable and reliable results are obtained

Existence of time-dependent density-functional theory for open electronic systems: time-dependent holographic electron density theorem.

Science.gov (United States)

Zheng, Xiao; Yam, ChiYung; Wang, Fan; Chen, GuanHua

2011-08-28

We present the time-dependent holographic electron density theorem (TD-HEDT), which lays the foundation of time-dependent density-functional theory (TDDFT) for open electronic systems. For any finite electronic system, the TD-HEDT formally establishes a one-to-one correspondence between the electron density inside any finite subsystem and the time-dependent external potential. As a result, any electronic property of an open system in principle can be determined uniquely by the electron density function inside the open region. Implications of the TD-HEDT on the practicality of TDDFT are also discussed.

Two-component scattering model and the electron density spectrum

Science.gov (United States)

Zhou, A. Z.; Tan, J. Y.; Esamdin, A.; Wu, X. J.

2010-02-01

In this paper, we discuss a rigorous treatment of the refractive scintillation caused by a two-component interstellar scattering medium and a Kolmogorov form of density spectrum. It is assumed that the interstellar scattering medium is composed of a thin-screen interstellar medium (ISM) and an extended interstellar medium. We consider the case that the scattering of the thin screen concentrates in a thin layer represented by a δ function distribution and that the scattering density of the extended irregular medium satisfies the Gaussian distribution. We investigate and develop equations for the flux density structure function corresponding to this two-component ISM geometry in the scattering density distribution and compare our result with the observations. We conclude that the refractive scintillation caused by this two-component ISM scattering gives a more satisfactory explanation for the observed flux density variation than does the single extended medium model. The level of refractive scintillation is strongly sensitive to the distribution of scattering material along the line of sight (LOS). The theoretical modulation indices are comparatively less sensitive to the scattering strength of the thin-screen medium, but they critically depend on the distance from the observer to the thin screen. The logarithmic slope of the structure function is sensitive to the scattering strength of the thin-screen medium, but is relatively insensitive to the thin-screen location. Therefore, the proposed model can be applied to interpret the structure functions of flux density observed in pulsar PSR B2111 + 46 and PSR B0136 + 57. The result suggests that the medium consists of a discontinuous distribution of plasma turbulence embedded in the interstellar medium. Thus our work provides some insight into the distribution of the scattering along the LOS to the pulsar PSR B2111 + 46 and PSR B0136 + 57.

Analytic energy gradients for the coupled-cluster singles and doubles method with the density-fitting approximation

International Nuclear Information System (INIS)

Bozkaya, Uğur; Sherrill, C. David

2016-01-01

An efficient implementation is presented for analytic gradients of the coupled-cluster singles and doubles (CCSD) method with the density-fitting approximation, denoted DF-CCSD. Frozen core terms are also included. When applied to a set of alkanes, the DF-CCSD analytic gradients are significantly accelerated compared to conventional CCSD for larger molecules. The efficiency of our DF-CCSD algorithm arises from the acceleration of several different terms, which are designated as the “gradient terms”: computation of particle density matrices (PDMs), generalized Fock-matrix (GFM), solution of the Z-vector equation, formation of the relaxed PDMs and GFM, back-transformation of PDMs and GFM to the atomic orbital (AO) basis, and evaluation of gradients in the AO basis. For the largest member of the alkane set (C 10 H 22 ), the computational times for the gradient terms (with the cc-pVTZ basis set) are 2582.6 (CCSD) and 310.7 (DF-CCSD) min, respectively, a speed up of more than 8-folds. For gradient related terms, the DF approach avoids the usage of four-index electron repulsion integrals. Based on our previous study [U. Bozkaya, J. Chem. Phys. 141, 124108 (2014)], our formalism completely avoids construction or storage of the 4-index two-particle density matrix (TPDM), using instead 2- and 3-index TPDMs. The DF approach introduces negligible errors for equilibrium bond lengths and harmonic vibrational frequencies.

MO-F-CAMPUS-J-04: Tissue Segmentation-Based MR Electron Density Mapping Method for MR-Only Radiation Treatment Planning of Brain

Energy Technology Data Exchange (ETDEWEB)

Yu, H [Sunnybrook Health Sciences Centre, Toronto, Ontario (Canada); Lee, Y [Sunnybrook Odette Cancer Centre, Toronto, Ontario (Canada); Ruschin, M [Odette Cancer Centre, Toronto, ON (Canada); Karam, I [Sunnybrook Odette Cancer Center, Toronto, Ontario (Canada); Sahgal, A [University of Toronto, Toronto, ON (Canada)

2015-06-15

Purpose: Automatically derive electron density of tissues using MR images and generate a pseudo-CT for MR-only treatment planning of brain tumours. Methods: 20 stereotactic radiosurgery (SRS) patients’ T1-weighted MR images and CT images were retrospectively acquired. First, a semi-automated tissue segmentation algorithm was developed to differentiate tissues with similar MR intensities and large differences in electron densities. The method started with approximately 12 slices of manually contoured spatial regions containing sinuses and airways, then air, bone, brain, cerebrospinal fluid (CSF) and eyes were automatically segmented using edge detection and anatomical information including location, shape, tissue uniformity and relative intensity distribution. Next, soft tissues - muscle and fat were segmented based on their relative intensity histogram. Finally, intensities of voxels in each segmented tissue were mapped into their electron density range to generate pseudo-CT by linearly fitting their relative intensity histograms. Co-registered CT was used as a ground truth. The bone segmentations of pseudo-CT were compared with those of co-registered CT obtained by using a 300HU threshold. The average distances between voxels on external edges of the skull of pseudo-CT and CT in three axial, coronal and sagittal slices with the largest width of skull were calculated. The mean absolute electron density (in Hounsfield unit) difference of voxels in each segmented tissues was calculated. Results: The average of distances between voxels on external skull from pseudo-CT and CT were 0.6±1.1mm (mean±1SD). The mean absolute electron density differences for bone, brain, CSF, muscle and fat are 78±114 HU, and 21±8 HU, 14±29 HU, 57±37 HU, and 31±63 HU, respectively. Conclusion: The semi-automated MR electron density mapping technique was developed using T1-weighted MR images. The generated pseudo-CT is comparable to that of CT in terms of anatomical position of

Density-matrix approach for the electroluminescence of molecules in a scanning tunneling microscope.

Science.gov (United States)

Tian, Guangjun; Liu, Ji-Cai; Luo, Yi

2011-04-29

The electroluminescence (EL) of molecules confined inside a nanocavity in the scanning tunneling microscope possesses many intriguing but unexplained features. We present here a general theoretical approach based on the density-matrix formalism to describe the EL from molecules near a metal surface induced by both electron tunneling and localized surface plasmon excitations simultaneously. It reveals the underlying physical mechanism for the external bias dependent EL. The important role played by the localized surface plasmon on the EL is highlighted. Calculations for porphyrin derivatives have reproduced corresponding experimental spectra and nicely explained the observed unusual large variation of emission spectral profiles. This general theoretical approach can find many applications in the design of molecular electronic and photonic devices.

Direct measurement of electron density in microdischarge at atmospheric pressure by Stark broadening

International Nuclear Information System (INIS)

Dong Lifang; Ran Junxia; Mao Zhiguo

2005-01-01

We present a method and results for measurement of electron density in atmospheric-pressure dielectric barrier discharge. The electron density of microdischarge in atmospheric pressure argon is measured by using the spectral line profile method. The asymmetrical deconvolution is used to obtain Stark broadening. The results show that the electron density in single filamentary microdischarge at atmospheric pressure argon is 3.05x10 15 cm -3 if the electron temperature is 10,000 K. The result is in good agreement with the simulation. The electron density in dielectric barrier discharge increases with the increase of applied voltage

Electron-impact study of PO2 using the R-matrix method

Science.gov (United States)

Bharadvaja, Anand; Kaur, Savinder; Baluja, K. L.

2013-06-01

The R-matrix approach is used to study the electron scattering from PO2 radical at low electron impact energies. The elastic scattering phenomenon is studied in static-exchange, one-state and many-states close-coupling approximation. The elastic differential cross sections, corresponding momentum-transfer cross sections, and collision frequency are calculated in the one-state configuration interaction approximation only. Calculations reveal a stable bound state of PO2- having symmetry 1A1, a configuration of ⋯8a12,⋯2b12,⋯5b22,⋯1a22, and vertical electron affinity of 2.94 eV. The excited state of anion PO2- having symmetry 3B1 is also just bound relative to the ground state of PO2 at its equilibrium geometry. The shape, core-excited, and Feshbach resonances are analyzed in different symmetries up to 7 eV. The partial waves up to l=4 are used to represent continuum electron. The converged cross sections are obtained for the partial waves having l greater than 4 by applying Born correction. Certain interesting spectroscopic properties of radical are also reported.

Equation satisfied by electron-electron mutual Coulomb repulsion energy density functional

OpenAIRE

Joubert, Daniel P.

2011-01-01

The electron-electron mutual Coulomb repulsion energy density functional satisfies an equation that links functionals and functional derivatives at N-electron and (N-1)-electron densities for densities determined from the same adiabatic scaled external potential for the N-electron system.

Electron density measurement for steady state plasmas

International Nuclear Information System (INIS)

Kawano, Yasunori; Chiba, Shinichi; Inoue, Akira

2000-01-01

Electron density of a large tokamak has been measured successfully by the tangential CO 2 laser polarimeter developed in JT-60U. The tangential Faraday rotation angles of two different wavelength of 9.27 and 10.6 μm provided the electron density independently. Two-color polarimeter concept for elimination of Faraday rotation at vacuum windows is verified for the first time. A system stability for long time operation up to ∼10 hours is confirmed. A fluctuation of a signal baseline is observed with a period of ∼3 hours and an amplitude of 0.4 - 0.7deg. In order to improve the polarimeter, an application of diamond window for reduction of the Faraday rotation at vacuum windows and another two-color polarimeter concept for elimination of mechanical rotation component are proposed. (author)

Robust Estimation of Electron Density From Anatomic Magnetic Resonance Imaging of the Brain Using a Unifying Multi-Atlas Approach

Energy Technology Data Exchange (ETDEWEB)

Ren, Shangjie [Tianjin Key Laboratory of Process Measurement and Control, School of Electrical Engineering and Automation, Tianjin University, Tianjin (China); Department of Radiation Oncology, Stanford University School of Medicine, Palo Alto, California (United States); Hara, Wendy; Wang, Lei; Buyyounouski, Mark K.; Le, Quynh-Thu; Xing, Lei [Department of Radiation Oncology, Stanford University School of Medicine, Palo Alto, California (United States); Li, Ruijiang, E-mail: rli2@stanford.edu [Department of Radiation Oncology, Stanford University School of Medicine, Palo Alto, California (United States)

2017-03-15

Purpose: To develop a reliable method to estimate electron density based on anatomic magnetic resonance imaging (MRI) of the brain. Methods and Materials: We proposed a unifying multi-atlas approach for electron density estimation based on standard T1- and T2-weighted MRI. First, a composite atlas was constructed through a voxelwise matching process using multiple atlases, with the goal of mitigating effects of inherent anatomic variations between patients. Next we computed for each voxel 2 kinds of conditional probabilities: (1) electron density given its image intensity on T1- and T2-weighted MR images; and (2) electron density given its spatial location in a reference anatomy, obtained by deformable image registration. These were combined into a unifying posterior probability density function using the Bayesian formalism, which provided the optimal estimates for electron density. We evaluated the method on 10 patients using leave-one-patient-out cross-validation. Receiver operating characteristic analyses for detecting different tissue types were performed. Results: The proposed method significantly reduced the errors in electron density estimation, with a mean absolute Hounsfield unit error of 119, compared with 140 and 144 (P<.0001) using conventional T1-weighted intensity and geometry-based approaches, respectively. For detection of bony anatomy, the proposed method achieved an 89% area under the curve, 86% sensitivity, 88% specificity, and 90% accuracy, which improved upon intensity and geometry-based approaches (area under the curve: 79% and 80%, respectively). Conclusion: The proposed multi-atlas approach provides robust electron density estimation and bone detection based on anatomic MRI. If validated on a larger population, our work could enable the use of MRI as a primary modality for radiation treatment planning.

Tuning the conductivity threshold and carrier density of two-dimensional electron gas at oxide interfaces through interface engineering

Directory of Open Access Journals (Sweden)

H. J. Harsan Ma

2015-08-01

Full Text Available The two-dimensional electron gas (2DEG formed at the perovskite oxides heterostructures is of great interest because of its potential applications in oxides electronics and nanoscale multifunctional devices. A canonical example is the 2DEG at the interface between a polar oxide LaAlO3 (LAO and non-polar SrTiO3 (STO. Here, the LAO polar oxide can be regarded as the modulating or doping layer and is expected to define the electronic properties of 2DEG at the LAO/STO interface. However, to practically implement the 2DEG in electronics and device design, desired properties such as tunable 2D carrier density are necessary. Here, we report the tuning of conductivity threshold, carrier density and electronic properties of 2DEG in LAO/STO heterostructures by insertion of a La0.5Sr0.5TiO3 (LSTO layer of varying thicknesses, and thus modulating the amount of polarization of the oxide over layers. Our experimental result shows an enhancement of carrier density up to a value of about five times higher than that observed at the LAO/STO interface. A complete thickness dependent metal-insulator phase diagram is obtained by varying the thickness of LAO and LSTO providing an estimate for the critical thickness needed for the metallic phase. The observations are discussed in terms of electronic reconstruction induced by polar oxides.

Electron density analysis of 1-butyl-3-methylimidazolium chloride ionic liquid.

Science.gov (United States)

del Olmo, Lourdes; Morera-Boado, Cercis; López, Rafael; García de la Vega, José M

2014-06-01

An analysis of the electron density of different conformers of the 1-butyl-3-methylimidazolium chloride (bmimCl) ionic liquid by using DFT through the BVP86 density functional has been obtained within the framework of Bader's atom in molecules (AIM), localized orbital locator (LOL), natural bond orbital (NBO), and deformed atoms in molecules (DAM). We also present an analysis of the reduced density gradients that deliver the non-covalent interaction regions and allow to understand the nature of intermolecular interactions. The most polar conformer can be characterized as ionic by AIM, LOL, and DAM methods while the most stable and the least polar shows shared-type interactions. The NBO method allows to comprehend what causes the stabilization of the most stable conformer based on analysis of the second-order perturbative energy and the charge transferred among the natural orbitals involved in the interaction.

Shack-Hartmann Electron Densitometer (SHED): An Optical System for Diagnosing Free Electron Density in Laser-Produced Plasmas

Science.gov (United States)

2016-11-01

Free Electron Density in Laser-Produced Plasmas by Anthony R Valenzuela Approved for public release; distribution is...AND SUBTITLE Shack-Hartmann Electron Densitometer (SHED): An Optical System for Diagnosing Free Electron Density in Laser-Produced Plasmas 5a...SUPPLEMENTARY NOTES 14. ABSTRACT The Shack-Hartmann Electron Densitometer is a novel method to diagnose ultrashort pulse laser–produced plasmas

Effective atomic numbers and electron density of dosimetric material

Directory of Open Access Journals (Sweden)

Kaginelli S

2009-01-01

Full Text Available A novel method for determination of mass attenuation coefficient of x-rays employing NaI (Tl detector system and radioactive sources is described.in this paper. A rigid geometry arrangement and gating of the spectrometer at FWHM position and selection of absorber foils are all done following detailed investigation, to minimize the effect of small angle scattering and multiple scattering on the mass attenuation coefficient, m/r, value. Firstly, for standardization purposes the mass attenuation coefficients of elemental foils such as Aluminum, Copper, Molybdenum, Tantalum and Lead are measured and then, this method is utilized for dosimetric interested material (sulfates. The experimental mass attenuation coefficient values are compared with the theoretical values to find good agreement between the theory and experiment within one to two per cent. The effective atomic numbers of the biological substitute material are calculated by sum rule and from the graph. The electron density of dosimetric material is calculated using the effective atomic number. The study has discussed in detail the attenuation coefficient, effective atomic number and electron density of dosimetric material/biological substitutes.

Anomalous evolution of Ar metastable density with electron density in high density Ar discharge

International Nuclear Information System (INIS)

Park, Min; Chang, Hong-Young; You, Shin-Jae; Kim, Jung-Hyung; Shin, Yong-Hyeon

2011-01-01

Recently, an anomalous evolution of argon metastable density with plasma discharge power (electron density) was reported [A. M. Daltrini, S. A. Moshkalev, T. J. Morgan, R. B. Piejak, and W. G. Graham, Appl. Phys. Lett. 92, 061504 (2008)]. Although the importance of the metastable atom and its density has been reported in a lot of literature, however, a basic physics behind the anomalous evolution of metastable density has not been clearly understood yet. In this study, we investigated a simple global model to elucidate the underlying physics of the anomalous evolution of argon metastable density with the electron density. On the basis of the proposed simple model, we reproduced the anomalous evolution of the metastable density and disclosed the detailed physics for the anomalous result. Drastic changes of dominant mechanisms for the population and depopulation processes of Ar metastable atoms with electron density, which take place even in relatively low electron density regime, is the clue to understand the result.

Two-mode Gaussian density matrices and squeezing of photons

International Nuclear Information System (INIS)

Tucci, R.R.

1992-01-01

In this paper, the authors generalize to 2-mode states the 1-mode state results obtained in a previous paper. The authors study 2-mode Gaussian density matrices. The authors find a linear transformation which maps the two annihilation operators, one for each mode, into two new annihilation operators that are uncorrelated and unsqueezed. This allows the authors to express the density matrix as a product of two 1-mode density matrices. The authors find general conditions under which 2-mode Gaussian density matrices become pure states. Possible pure states include the 2-mode squeezed pure states commonly mentioned in the literature, plus other pure states never mentioned before. The authors discuss the entropy and thermodynamic laws (Second Law, Fundamental Equation, and Gibbs-Duhem Equation) for the 2-mode states being considered

Calculation of dynamic and electronic properties of perfect and defect crystals by semiempirical quantum mechanical methods

International Nuclear Information System (INIS)

Zunger, A.

1975-07-01

Semiempirical all-valence-electron LCAO methods, that were previously used to study the electronic structure of molecules are applied to three problems in solid state physics: the electronic band structure of covalent crystals, point defect problems in solids and lattice dynamical study of molecular crystals. Calculation methods for the electronic band structure of regular solids are introduced and problems regarding the computation of the density matrix in solids are discussed. Three models for treating the electronic eigenvalue problem in the solid, within the proposed calculation schemes, are discussed and the proposed models and calculation schemes are applied to the calculation of the electronic structure of several solids belonging to different crystal types. The calculation models also describe electronic properties of deep defects in covalent insulating crystals. The possible usefulness of the semieipirical LCAO methods in determining the first order intermolecular interaction potential in solids and an improved model for treating the lattice dynamics and related thermodynamical properties of molecular solids are presented. The improved lattice dynamical is used to compute phonon dispersion curves, phonon density of states, stable unit cell structure, lattice heat capacity and thermal crystal parameters, in α and γ-N 2 crystals, using the N 2 -N 2 intermolecular interaction potential that has been computed from the semiempirical LCAO methods. (B.G.)

Efficient mixing scheme for self-consistent all-electron charge density

Science.gov (United States)

Shishidou, Tatsuya; Weinert, Michael

2015-03-01

In standard ab initio density-functional theory calculations, the charge density ρ is gradually updated using the ``input'' and ``output'' densities of the current and previous iteration steps. To accelerate the convergence, Pulay mixing has been widely used with great success. It expresses an ``optimal'' input density ρopt and its ``residual'' Ropt by a linear combination of the densities of the iteration sequences. In large-scale metallic systems, however, the long range nature of Coulomb interaction often causes the ``charge sloshing'' phenomenon and significantly impacts the convergence. Two treatments, represented in reciprocal space, are known to suppress the sloshing: (i) the inverse Kerker metric for Pulay optimization and (ii) Kerker-type preconditioning in mixing Ropt. In all-electron methods, where the charge density does not have a converging Fourier representation, treatments equivalent or similar to (i) and (ii) have not been described so far. In this work, we show that, by going through the calculation of Hartree potential, one can accomplish the procedures (i) and (ii) without entering the reciprocal space. Test calculations are done with a FLAPW method.

Spin Density Matrix Elements in exclusive production of ω mesons at Hermes

Directory of Open Access Journals (Sweden)

Marianski B.

2014-03-01

Full Text Available Spin density matrix elements have been determined for exclusive ω meson production on hydrogen and deuterium targets, in the kinematic region of 1.0 < Q2 < 10.0 GeV2, 3.0 < W < 6.3 GeV and –t' < 0.2 GeV2. The data, from which SDMEs are determined, were accumulated with the HERMES forward spectrometer during the running period of 1996 to 2007 using the 27.6 GeV electron or positron beam of HERA. A sizable contribution of unnatural parity exchange amplitudes is found for exclusive ω meson production.

The method of the transfer matrix applied to the study of the electronic properties of the Kronig- Penney model with structural disorder

International Nuclear Information System (INIS)

Miranda S, Anabel; Landauro S, Carlos

2008-01-01

In the present work the transfer matrix method is employed to study the electronic properties of the Kronig-Penney model including disorder in the periodic system. The results show that although the electronic properties are very similar to the corresponding periodic case, disorder in the system produces a decrease of the transmission in the whole range of energies which indicates clearly a reduction of the electronic transport (conductivity) due to the disorder in the system. (author)

Two methods for isolating the lung area of a CT scan for density information

International Nuclear Information System (INIS)

Hedlund, L.W.; Anderson, R.F.; Goulding, P.L.; Beck, J.W.; Effmann, E.L.; Putman, C.E.

1982-01-01

Extracting density information from irregularly shaped tissue areas of CT scans requires automated methods when many scans are involved. We describe two computer methods that automatically isolate the lung area of a CT scan. Each starts from a single, operator specified point in the lung. The first method follows the steep density gradient boundary between lung and adjacent tissues; this tracking method is useful for estimating the overall density and total area of lung in a scan because all pixels within the lung area are available for statistical sampling. The second method finds all contiguous pixels of lung that are within the CT number range of air to water and are not a part of strong density gradient edges; this method is useful for estimating density and area of the lung parenchyma. Structures within the lung area that are surrounded by strong density gradient edges, such as large blood vessels, airways and nodules, are excluded from the lung sample while lung areas with diffuse borders, such as an area of mild or moderate edema, are retained. Both methods were tested on scans from an animal model of pulmonary edema and were found to be effective in isolating normal and diseased lungs. These methods are also suitable for isolating other organ areas of CT scans that are bounded by density gradient edges

The J-Matrix Method Developments and Applications

CERN Document Server

Alhaidari, Abdulaziz D; Heller, Eric J; Abdelmonem, Mohamed S

2008-01-01

This volume aims to provide the fundamental knowledge to appreciate the advantages of the J-matrix method and to encourage its use and further development. The J-matrix method is an algebraic method of quantum scattering with substantial success in atomic and nuclear physics. The accuracy and convergence property of the method compares favourably with other successful scattering calculation methods. Despite its thirty-year long history new applications are being found for the J-matrix method. This book gives a brief account of the recent developments and some selected applications of the method in atomic and nuclear physics. New findings are reported in which experimental results are compared to theoretical calculations. Modifications, improvements and extensions of the method are discussed using the language of the J-matrix. The volume starts with a Foreword by the two co-founders of the method, E.J. Heller and H.A. Yamani and it contains contributions from 24 prominent international researchers.

Electronic structure and electron momentum density in TiSi

Energy Technology Data Exchange (ETDEWEB)

Ghaleb, A.M. [Department of Physics, College of Science, University of Kirkuk, Kirkuk (Iraq); Mohammad, F.M. [Department of Physics, College of Science, University of Tikreet, Tikreet (Iraq); Sahariya, Jagrati [Department of Physics, University College of Science, M.L. Sukhadia University, Udaipur 313001, Rajasthan (India); Sharma, Mukesh [Physics Division, Forensic Science Laboratory, Jaipur, Rajasthan (India); Ahuja, B.L., E-mail: blahuja@yahoo.com [Department of Physics, University College of Science, M.L. Sukhadia University, Udaipur 313001, Rajasthan (India)

2013-03-01

We report the electron momentum density in titanium monosilicide using {sup 241}Am Compton spectrometer. Experimental Compton profile has been compared with the theoretical profiles computed using linear combination of atomic orbitals (LCAO). The energy bands, density of states and Fermi surface structures of TiSi are reported using the LCAO and the full potential linearized augmented plane wave methods. Theoretical anisotropies in directional Compton profiles are interpreted in terms of energy bands. To confirm the conducting behavior, we also report the real space analysis of experimental Compton profile of TiSi.

Calculation of the two-electron Darwin term using explicitly correlated wave functions

International Nuclear Information System (INIS)

Middendorf, Nils; Höfener, Sebastian; Klopper, Wim; Helgaker, Trygve

2012-01-01

Graphical abstract: The two-electron Darwin term is computed analytically at the MP2-F12 level of theory using density fitted integrals. Highlights: ► Two-electron Darwin term computed analytically at the MP2-F12 level. ► Darwin two-electron integrals computed using density fitting techniques. ► Two-electron Darwin term dominated by singlet pair contributions. ► Much improved basis set convergence is achieved with F12 methods. ► Interference correction works well for the two-electron Darwin term. - Abstract: This article is concerned with the calculation of the two-electron Darwin term (D2). At the level of explicitly correlated second-order perturbation theory (MP2-F12), the D2 term is obtained as an analytic energy derivative; at the level of explicitly correlated coupled-cluster theory, it is obtained from finite differences. To avoid the calculation of four-center integrals, a density-fitting approximation is applied to the D2 two-electron integrals without loss of accuracy, even though the absolute value of the D2 term is typically about 0.1 mE h . Explicitly correlated methods provide a qualitatively correct description of the short-range region around the Coulomb hole, even for small orbital basis sets. Therefore, explicitly correlated wave functions remedy the otherwise extremely slow convergence of the D2 contribution with respect to the basis-set size, yielding more accurate results than those obtained by two-point basis-set extrapolation. Moreover, we show that the interference correction of Petersson’s complete-basis-set model chemistry can be used to compute a D2 basis-set correction at the MP2-F12 level to improve standard coupled-cluster singles-and-doubles results.

Evaluation of macromolecular electron-density map quality using the correlation of local r.m.s. density

International Nuclear Information System (INIS)

Terwilliger, Thomas C.; Berendzen, Joel

1999-01-01

The correlation of local r.m.s. density is shown to be a good measure of the presence of distinct solvent and macromolecule regions in macromolecular electron-density maps. It has recently been shown that the standard deviation of local r.m.s. electron density is a good indicator of the presence of distinct regions of solvent and protein in macromolecular electron-density maps [Terwilliger & Berendzen (1999 ▶). Acta Cryst. D55, 501–505]. Here, it is demonstrated that a complementary measure, the correlation of local r.m.s. density in adjacent regions on the unit cell, is also a good measure of the presence of distinct solvent and protein regions. The correlation of local r.m.s. density is essentially a measure of how contiguous the solvent (and protein) regions are in the electron-density map. This statistic can be calculated in real space or in reciprocal space and has potential uses in evaluation of heavy-atom solutions in the MIR and MAD methods as well as for evaluation of trial phase sets in ab initio phasing procedures

Projected quasiparticle theory for molecular electronic structure

Science.gov (United States)

Scuseria, Gustavo E.; Jiménez-Hoyos, Carlos A.; Henderson, Thomas M.; Samanta, Kousik; Ellis, Jason K.

2011-09-01

We derive and implement symmetry-projected Hartree-Fock-Bogoliubov (HFB) equations and apply them to the molecular electronic structure problem. All symmetries (particle number, spin, spatial, and complex conjugation) are deliberately broken and restored in a self-consistent variation-after-projection approach. We show that the resulting method yields a comprehensive black-box treatment of static correlations with effective one-electron (mean-field) computational cost. The ensuing wave function is of multireference character and permeates the entire Hilbert space of the problem. The energy expression is different from regular HFB theory but remains a functional of an independent quasiparticle density matrix. All reduced density matrices are expressible as an integration of transition density matrices over a gauge grid. We present several proof-of-principle examples demonstrating the compelling power of projected quasiparticle theory for quantum chemistry.

Porting of the DBCSR library for Sparse Matrix-Matrix Multiplications to Intel Xeon Phi systems

OpenAIRE

Bethune, Iain; Gloess, Andeas; Hutter, Juerg; Lazzaro, Alfio; Pabst, Hans; Reid, Fiona

2017-01-01

Multiplication of two sparse matrices is a key operation in the simulation of the electronic structure of systems containing thousands of atoms and electrons. The highly optimized sparse linear algebra library DBCSR (Distributed Block Compressed Sparse Row) has been specifically designed to efficiently perform such sparse matrix-matrix multiplications. This library is the basic building block for linear scaling electronic structure theory and low scaling correlated methods in CP2K. It is para...

Density-dependent electron transport and precise modeling of GaN high electron mobility transistors

Energy Technology Data Exchange (ETDEWEB)

Bajaj, Sanyam, E-mail: bajaj.10@osu.edu; Shoron, Omor F.; Park, Pil Sung; Krishnamoorthy, Sriram; Akyol, Fatih; Hung, Ting-Hsiang [Department of Electrical and Computer Engineering, The Ohio State University, Columbus, Ohio 43210 (United States); Reza, Shahed; Chumbes, Eduardo M. [Raytheon Integrated Defense Systems, Andover, Massachusetts 01810 (United States); Khurgin, Jacob [Department of Electrical and Computer Engineering, Johns Hopkins University, Baltimore, Maryland 21218 (United States); Rajan, Siddharth [Department of Electrical and Computer Engineering, The Ohio State University, Columbus, Ohio 43210 (United States); Department of Material Science and Engineering, The Ohio State University, Columbus, Ohio 43210 (United States)

2015-10-12

We report on the direct measurement of two-dimensional sheet charge density dependence of electron transport in AlGaN/GaN high electron mobility transistors (HEMTs). Pulsed IV measurements established increasing electron velocities with decreasing sheet charge densities, resulting in saturation velocity of 1.9 × 10{sup 7 }cm/s at a low sheet charge density of 7.8 × 10{sup 11 }cm{sup −2}. An optical phonon emission-based electron velocity model for GaN is also presented. It accommodates stimulated longitudinal optical (LO) phonon emission which clamps the electron velocity with strong electron-phonon interaction and long LO phonon lifetime in GaN. A comparison with the measured density-dependent saturation velocity shows that it captures the dependence rather well. Finally, the experimental result is applied in TCAD-based device simulator to predict DC and small signal characteristics of a reported GaN HEMT. Good agreement between the simulated and reported experimental results validated the measurement presented in this report and established accurate modeling of GaN HEMTs.

Density-dependent electron transport and precise modeling of GaN high electron mobility transistors

International Nuclear Information System (INIS)

Bajaj, Sanyam; Shoron, Omor F.; Park, Pil Sung; Krishnamoorthy, Sriram; Akyol, Fatih; Hung, Ting-Hsiang; Reza, Shahed; Chumbes, Eduardo M.; Khurgin, Jacob; Rajan, Siddharth

2015-01-01

We report on the direct measurement of two-dimensional sheet charge density dependence of electron transport in AlGaN/GaN high electron mobility transistors (HEMTs). Pulsed IV measurements established increasing electron velocities with decreasing sheet charge densities, resulting in saturation velocity of 1.9 × 10 7  cm/s at a low sheet charge density of 7.8 × 10 11  cm −2 . An optical phonon emission-based electron velocity model for GaN is also presented. It accommodates stimulated longitudinal optical (LO) phonon emission which clamps the electron velocity with strong electron-phonon interaction and long LO phonon lifetime in GaN. A comparison with the measured density-dependent saturation velocity shows that it captures the dependence rather well. Finally, the experimental result is applied in TCAD-based device simulator to predict DC and small signal characteristics of a reported GaN HEMT. Good agreement between the simulated and reported experimental results validated the measurement presented in this report and established accurate modeling of GaN HEMTs

Properties of short-range and long-range correlation energy density functionals from electron-electron coalescence

International Nuclear Information System (INIS)

Gori-Giorgi, Paola; Savin, Andreas

2006-01-01

The combination of density-functional theory with other approaches to the many-electron problem through the separation of the electron-electron interaction into a short-range and a long-range contribution is a promising method, which is raising more and more interest in recent years. In this work some properties of the corresponding correlation energy functionals are derived by studying the electron-electron coalescence condition for a modified (long-range-only) interaction. A general relation for the on-top (zero electron-electron distance) pair density is derived, and its usefulness is discussed with some examples. For the special case of the uniform electron gas, a simple parametrization of the on-top pair density for a long-range only interaction is presented and supported by calculations within the ''extended Overhauser model.'' The results of this work can be used to build self-interaction corrected short-range correlation energy functionals

Efficient Computation of Sparse Matrix Functions for Large-Scale Electronic Structure Calculations: The CheSS Library.

Science.gov (United States)

Mohr, Stephan; Dawson, William; Wagner, Michael; Caliste, Damien; Nakajima, Takahito; Genovese, Luigi

2017-10-10

We present CheSS, the "Chebyshev Sparse Solvers" library, which has been designed to solve typical problems arising in large-scale electronic structure calculations using localized basis sets. The library is based on a flexible and efficient expansion in terms of Chebyshev polynomials and presently features the calculation of the density matrix, the calculation of matrix powers for arbitrary powers, and the extraction of eigenvalues in a selected interval. CheSS is able to exploit the sparsity of the matrices and scales linearly with respect to the number of nonzero entries, making it well-suited for large-scale calculations. The approach is particularly adapted for setups leading to small spectral widths of the involved matrices and outperforms alternative methods in this regime. By coupling CheSS to the DFT code BigDFT, we show that such a favorable setup is indeed possible in practice. In addition, the approach based on Chebyshev polynomials can be massively parallelized, and CheSS exhibits excellent scaling up to thousands of cores even for relatively small matrix sizes.

Postsynaptic density protein 95 in the striosome and matrix compartments of the human neostriatum.

Directory of Open Access Journals (Sweden)

Ryoma eMorigaki

2015-11-01

Full Text Available The human neostriatum consists of two functional subdivisions referred to as the striosome (patch and matrix compartments. The striosome-matrix dopamine systems play a central role in cortico-thalamo-basal ganglia circuits, and their involvement is thought to underlie the genesis of multiple movement and behavioral disorders, and of drug addiction. Human neuropathology also has shown that striosomes and matrix have differential vulnerability patterns in several striatal neurodegenerative diseases. Postsynaptic density protein 95 (PSD-95, also known as DLG4, is a major scaffolding protein in the postsynaptic densities of dendritic spines. PSD-95 is now known to negatively regulate not only N-methyl-D-aspartate glutamate signaling, but also dopamine D1 signals at sites of postsynaptic transmission. Accordingly, a neuroprotective role for PSD-95 against dopamine D1 receptor (D1R-mediated neurotoxicity in striatal neurodegeneration also has been suggested. Here, we used a highly sensitive immunohistochemistry technique to show that in the human neostriatum, PSD-95 is differentially concentrated in the striosome and matrix compartments, with a higher density of PSD-95 labeling in the matrix compartment than in the striosomes. This compartment-specific distribution of PSD-95 was strikingly complementary to that of D1R. In addition to the possible involvement of PSD-95-mediated synaptic function in compartment-specific dopamine signals, we suggest that the striosomes might be more susceptible to D1R-mediated neurotoxicity than the matrix compartment. This notion may provide new insight into the compartment-specific vulnerability of MSNs in striatal neurodegeneration.

A novel loss reduced modulation strategy for matrix converters

DEFF Research Database (Denmark)

Helle, Lars; Munk-Nielsen, Stig

2001-01-01

This paper presents a new modulation strategy for three-phase to three-phase matrix converters. The new modulation strategy is applicable whenever the output voltage reference is below half the input voltage. By applying this new modulation method, the switching losses are reduced by 15-35% compa...

Transition densities with electron scattering

International Nuclear Information System (INIS)

Heisenberg, J.

1985-01-01

This paper reviews the ground state and transition charge densities in nuclei via electron scattering. Using electrons as a spectroscopic tool in nuclear physics, these transition densities can be determined with high precision, also in the nuclear interior. These densities generally ask for a microscopic interpretation in terms of contributions from individual nucleons. The results for single particle transitions confirm the picture of particle-phonon coupling. (Auth.)

Departures from local thermodynamic equilibrium in cutting arc plasmas derived from electron and gas density measurements using a two-wavelength quantitative Schlieren technique

International Nuclear Information System (INIS)

Prevosto, L.; Mancinelli, B.; Artana, G.; Kelly, H.

2011-01-01

A two-wavelength quantitative Schlieren technique that allows inferring the electron and gas densities of axisymmetric arc plasmas without imposing any assumption regarding statistical equilibrium models is reported. This technique was applied to the study of local thermodynamic equilibrium (LTE) departures within the core of a 30 A high-energy density cutting arc. In order to derive the electron and heavy particle temperatures from the inferred density profiles, a generalized two-temperature Saha equation together with the plasma equation of state and the quasineutrality condition were employed. Factors such as arc fluctuations that influence the accuracy of the measurements and the validity of the assumptions used to derive the plasma species temperature were considered. Significant deviations from chemical equilibrium as well as kinetic equilibrium were found at elevated electron temperatures and gas densities toward the arc core edge. An electron temperature profile nearly constant through the arc core with a value of about 14000-15000 K, well decoupled from the heavy particle temperature of about 1500 K at the arc core edge, was inferred.

Optimal pinnate leaf-like network/matrix structure for enhanced conductive cooling

International Nuclear Information System (INIS)

Hu, Liguo; Zhou, Han; Zhu, Hanxing; Fan, Tongxiang; Zhang, Di

2015-01-01

Highlights: • We present a pinnate leaf-like network/matrix structure for conductive cooling. • We study the effect of matrix thickness on network conductive cooling performance. • Matrix thickness determines optimal distance between collection channels in network. • We determine the optimal network architecture from a global perspective. • Optimal network greatly reduces the maximum temperature difference in the network. - Abstract: Heat generated in electronic devices has to be effectively removed because excessive temperature strongly impairs their performance and reliability. Embedding a high thermal conductivity network into an electronic device is an effective method to conduct the generated heat to the outside. In this study, inspired by the pinnate leaf, we present a pinnate leaf-like network embedded in the matrix (i.e., electronic device) to cool the matrix by conduction and develop a method to construct the optimal network. In this method, we first investigate the effect of the matrix thickness on the conductive cooling performance of the network, and then optimize the network architecture from a global perspective so that to minimize the maximum temperature difference between the heat sink and the matrix. The results indicate that the matrix thickness determines the optimal distance of the neighboring collection channels in the network, which minimizes the maximum temperature difference between the matrix and the network, and that the optimal network greatly reduces the maximum temperature difference in the network. The results can serve as a design guide for efficient conductive cooling of electronic devices

Interaction of two solitary waves in quantum electron-positron-ion plasma

International Nuclear Information System (INIS)

Xu Yanxia; Lin Maimai; Shi Yuren; Duan Wenshan; Liu Zongming; Chen Jianmin

2011-01-01

The collision between two ion-acoustic solitary waves with arbitrary colliding angle θ in an unmagnetized, ultracold quantum three-component e-p-i plasma has been investigated. By using the extended Poincare-Lighthill-Kuo (PLK) perturbation method, we obtain the KdV equations and the analytical phase shifts after the collision of two solitary waves in this three-component plasma. The effects of the quantum parameter H, the ratio of Fermi positron temperature to Fermi electron temperature σ, the ratio of Fermi positron number density to Fermi electron number density μ, and the ratio of Fermi ion temperature to Fermi electron temperature ρ on the phase shifts are studied. It is found that these parameters can significantly influence the phase shifts of the solitons.

Electron beam cross-linking of natural rubber/linear-low density polyethylene blends

International Nuclear Information System (INIS)

Ahmad, A.; Mohd, D. H.; Abdullah, I.

2005-01-01

Effects of electron beam irradiation on the mechanical properties and morphological structure of natural rubber/linear-low density polyethylene blend was investigated The natural rubber/linear-low density polyethylene blend was prepared by melt blending in a Haake internal mixer at 140 d ig C , rotor speed of 50 rpm, and in 15 min Liquid natural rubber was incorporated into the blend as a compatibilizer Samples in the form of 1 mm sheets were exposed to 50-300 kGy of electron beam irradiation and analyzed for swelling index and gel content, tensile strength, and surface morphology. The result Indicated that gel content and mechanical properties of the samples increased with radiation dosage. The honey-comb structure of the surface morphology in low dosage irradiated samples slowly transformed into a continuous matrix on increasing radiation dose The variation of mechanical and physical properties was due to Increase in cross-linking density in the rubber and plastic phases and rubber-plastic Interaction on irradiation

Future directions in electron--ion collision physics

International Nuclear Information System (INIS)

Reed, K.J.; Griffin, D.C.

1992-01-01

This report discusses the following topics: Summary of session on synergistic co-ordination of theory and experiment; synergism between experiment and theory in atomic physics; comparison of theory and experiment for electron-ion excitation and ionization; summary of session on new theoretical and computational methods; new theoretical and computational methods-r-matrix calculations; the coulomb three-body problem: a progress report; summary of session on needs and applications for electron-ion collisional data; electron-ion collisions in the plasma edge; needs and applications of theoretical data for electron impact excitation; summary of session on relativistic effects, indirect effects, resonance, etc; direct and resonant processes in electron-ion collisions; relativistic calculations of electron impact ionization and dielectronic recombination cross section for highly charged ions; electron-ion recombination in the close-coupling approximation; modified resonance amplitudes with strongly correlated channels; a density-matrix approach to the broadening of spectral lines by autoionization, radiative transitions and electron-ion collisions; towards a time-dependent description of electron-atom/ion collisions two electron systems; and comments on inclusion of the generalized bright interaction in electron impact excitation of highly charged ions

On extending Kohn-Sham density functionals to systems with fractional number of electrons.

Science.gov (United States)

Li, Chen; Lu, Jianfeng; Yang, Weitao

2017-06-07

We analyze four ways of formulating the Kohn-Sham (KS) density functionals with a fractional number of electrons, through extending the constrained search space from the Kohn-Sham and the generalized Kohn-Sham (GKS) non-interacting v-representable density domain for integer systems to four different sets of densities for fractional systems. In particular, these density sets are (I) ensemble interacting N-representable densities, (II) ensemble non-interacting N-representable densities, (III) non-interacting densities by the Janak construction, and (IV) non-interacting densities whose composing orbitals satisfy the Aufbau occupation principle. By proving the equivalence of the underlying first order reduced density matrices associated with these densities, we show that sets (I), (II), and (III) are equivalent, and all reduce to the Janak construction. Moreover, for functionals with the ensemble v-representable assumption at the minimizer, (III) reduces to (IV) and thus justifies the previous use of the Aufbau protocol within the (G)KS framework in the study of the ground state of fractional electron systems, as defined in the grand canonical ensemble at zero temperature. By further analyzing the Aufbau solution for different density functional approximations (DFAs) in the (G)KS scheme, we rigorously prove that there can be one and only one fractional occupation for the Hartree Fock functional, while there can be multiple fractional occupations for general DFAs in the presence of degeneracy. This has been confirmed by numerical calculations using the local density approximation as a representative of general DFAs. This work thus clarifies important issues on density functional theory calculations for fractional electron systems.

Energy-momentum density of graphite by electron-momentum spectroscopy

International Nuclear Information System (INIS)

Vos, M.; Fang, Z.; Canney, S.; Kheifets, A.; McCarthy, I.E.; Weigold, E.

1996-11-01

The energy-resolved electron momentum density of graphite has been measured along a series of well-defined directions using electron momentum spectroscopy (EMS). This is the first measurement of this kind performed on a single-crystal target with a thoroughly controlled orientation which clearly demonstrates the different nature of the σ and π bands in graphite. Good agreement between the calculated density and the measured one is found, further establishing that fact that EMS yields more direct and complete information on the valence electronic structure that any other method. 12 refs., 2 figs

Theory of sum-frequency generation spectroscopy of adsorbed molecules using the density matrix method-broadband vibrational sum-frequency generation and applications

International Nuclear Information System (INIS)

Bonn, M; Ueba, H; Wolf, M

2005-01-01

A generalized theory of frequency- and time-resolved vibrational sum-frequency generation (SFG) spectroscopy of adsorbates at surfaces is presented using the density matrix formalism. Our theoretical treatment is specifically aimed at addressing issues that accompany the relatively novel SFG approach using broadband infrared pulses. The ultrashort duration of these pulses makes them ideally suited for time-resolved investigations, for which we present a complete theoretical treatment. A second key characteristic of these pulses is their large bandwidth and high intensity, which allow for highly non-linear effects, including vibrational ladder climbing of surface vibrations. We derive general expressions relating the density matrix to SFG spectra, and apply these expressions to specific experimental results by solving the coupled optical Bloch equations of the density matrix elements. Thus, we can theoretically reproduce recent experimentally demonstrated hot band SFG spectra using femtosecond broadband infrared excitation of carbon monoxide (CO) on a Ru(001) surface

M3: Matrix Multiplication on MapReduce

DEFF Research Database (Denmark)

Silvestri, Francesco; Ceccarello, Matteo

2015-01-01

M3 is an Hadoop library for performing dense and sparse matrix multiplication in MapReduce. The library is based on multi-round algorithms exploiting the 3D decomposition of the problem.......M3 is an Hadoop library for performing dense and sparse matrix multiplication in MapReduce. The library is based on multi-round algorithms exploiting the 3D decomposition of the problem....

Exact many-body dynamics with stochastic one-body density matrix evolution

International Nuclear Information System (INIS)

Lacroix, D.

2004-05-01

In this article, we discuss some properties of the exact treatment of the many-body problem with stochastic Schroedinger equation (SSE). Starting from the SSE theory, an equivalent reformulation is proposed in terms of quantum jumps in the density matrix space. The technical details of the derivation a stochastic version of the Liouville von Neumann equation are given. It is shown that the exact Many-Body problem could be replaced by an ensemble of one-body density evolution, where each density matrix evolves according to its own mean-field augmented by a one-body noise. (author)

Oriented Polar Molecules in a Solid Inert-Gas Matrix: A Proposed Method for Measuring the Electric Dipole Moment of the Electron

Directory of Open Access Journals (Sweden)

A. C. Vutha

2018-01-01

Full Text Available We propose a very sensitive method for measuring the electric dipole moment of the electron using polar molecules embedded in a cryogenic solid matrix of inert-gas atoms. The polar molecules can be oriented in the z ^ -direction by an applied electric field, as has recently been demonstrated by Park et al. The trapped molecules are prepared into a state that has its electron spin perpendicular to z ^ , and a magnetic field along z ^ causes precession of this spin. An electron electric dipole moment d e would affect this precession due to the up to 100 GV/cm effective electric field produced by the polar molecule. The large number of polar molecules that can be embedded in a matrix, along with the expected long coherence times for the precession, allows for the possibility of measuring d e to an accuracy that surpasses current measurements by many orders of magnitude. Because the matrix can inhibit molecular rotations and lock the orientation of the polar molecules, it may not be necessary to have an electric field present during the precession. The proposed technique can be applied using a variety of polar molecules and inert gases, which, along with other experimental variables, should allow for careful study of systematic uncertainties in the measurement.

Nuclear-electronic orbital reduced explicitly correlated Hartree-Fock approach: Restricted basis sets and open-shell systems

International Nuclear Information System (INIS)

Brorsen, Kurt R.; Sirjoosingh, Andrew; Pak, Michael V.; Hammes-Schiffer, Sharon

2015-01-01

The nuclear electronic orbital (NEO) reduced explicitly correlated Hartree-Fock (RXCHF) approach couples select electronic orbitals to the nuclear orbital via Gaussian-type geminal functions. This approach is extended to enable the use of a restricted basis set for the explicitly correlated electronic orbitals and an open-shell treatment for the other electronic orbitals. The working equations are derived and the implementation is discussed for both extensions. The RXCHF method with a restricted basis set is applied to HCN and FHF − and is shown to agree quantitatively with results from RXCHF calculations with a full basis set. The number of many-particle integrals that must be calculated for these two molecules is reduced by over an order of magnitude with essentially no loss in accuracy, and the reduction factor will increase substantially for larger systems. Typically, the computational cost of RXCHF calculations with restricted basis sets will scale in terms of the number of basis functions centered on the quantum nucleus and the covalently bonded neighbor(s). In addition, the RXCHF method with an odd number of electrons that are not explicitly correlated to the nuclear orbital is implemented using a restricted open-shell formalism for these electrons. This method is applied to HCN + , and the nuclear densities are in qualitative agreement with grid-based calculations. Future work will focus on the significance of nonadiabatic effects in molecular systems and the further enhancement of the NEO-RXCHF approach to accurately describe such effects

Nuclear-electronic orbital reduced explicitly correlated Hartree-Fock approach: Restricted basis sets and open-shell systems

Energy Technology Data Exchange (ETDEWEB)

Brorsen, Kurt R.; Sirjoosingh, Andrew; Pak, Michael V.; Hammes-Schiffer, Sharon, E-mail: shs3@illinois.edu [Department of Chemistry, University of Illinois at Urbana-Champaign, 600 South Mathews Ave., Urbana, Illinois 61801 (United States)

2015-06-07

The nuclear electronic orbital (NEO) reduced explicitly correlated Hartree-Fock (RXCHF) approach couples select electronic orbitals to the nuclear orbital via Gaussian-type geminal functions. This approach is extended to enable the use of a restricted basis set for the explicitly correlated electronic orbitals and an open-shell treatment for the other electronic orbitals. The working equations are derived and the implementation is discussed for both extensions. The RXCHF method with a restricted basis set is applied to HCN and FHF{sup −} and is shown to agree quantitatively with results from RXCHF calculations with a full basis set. The number of many-particle integrals that must be calculated for these two molecules is reduced by over an order of magnitude with essentially no loss in accuracy, and the reduction factor will increase substantially for larger systems. Typically, the computational cost of RXCHF calculations with restricted basis sets will scale in terms of the number of basis functions centered on the quantum nucleus and the covalently bonded neighbor(s). In addition, the RXCHF method with an odd number of electrons that are not explicitly correlated to the nuclear orbital is implemented using a restricted open-shell formalism for these electrons. This method is applied to HCN{sup +}, and the nuclear densities are in qualitative agreement with grid-based calculations. Future work will focus on the significance of nonadiabatic effects in molecular systems and the further enhancement of the NEO-RXCHF approach to accurately describe such effects.

Density-density functionals and effective potentials in many-body electronic structure calculations

International Nuclear Information System (INIS)

Reboredo, Fernando A.; Kent, Paul R.

2008-01-01

We demonstrate the existence of different density-density functionals designed to retain selected properties of the many-body ground state in a non-interacting solution starting from the standard density functional theory ground state. We focus on diffusion quantum Monte Carlo applications that require trial wave functions with optimal Fermion nodes. The theory is extensible and can be used to understand current practices in several electronic structure methods within a generalized density functional framework. The theory justifies and stimulates the search of optimal empirical density functionals and effective potentials for accurate calculations of the properties of real materials, but also cautions on the limits of their applicability. The concepts are tested and validated with a near-analytic model.

A two-loop test of M(atrix) theory

International Nuclear Information System (INIS)

Becker, K.

1997-01-01

We consider the scattering of two Dirichlet zero-branes in M(atrix) theory. Using the formulation of M(atrix) theory in terms of ten-dimensional super Yang-Mills theory dimensionally reduced to (0+1) dimensions, we obtain the effective (velocity-dependent) potential describing these particles. At one loop we obtain the well-known result for the leading order of the effective potential V eff ∝v 4 /r 7 , where v and r are the relative velocity and distance between the two zero-branes, respectively. A calculation of the effective potential at two loops shows that no renormalizations of the v 4 term of the effective potential occur at this order. (orig.)

Wave functions and low-order density matrices for a class of two-electron 'artificial atoms' embracing Hookean and Moshinsky models

International Nuclear Information System (INIS)

Holas, A.; Howard, I.A.; March, N.H.

2003-01-01

A class of model two-electron 'artificial atoms' is proposed which embraces both Hookean and Moshinsky models. Particle densities and spinless first-order density matrices are obtained for this class of models. These quantities and the interacting system kinetic energy can be calculated using the ground-state solution of an explicit single-particle radial Schroedinger equation

Methods for converging correlation energies within the dielectric matrix formalism

Science.gov (United States)

Dixit, Anant; Claudot, Julien; Gould, Tim; Lebègue, Sébastien; Rocca, Dario

2018-03-01

Within the dielectric matrix formalism, the random-phase approximation (RPA) and analogous methods that include exchange effects are promising approaches to overcome some of the limitations of traditional density functional theory approximations. The RPA-type methods however have a significantly higher computational cost, and, similarly to correlated quantum-chemical methods, are characterized by a slow basis set convergence. In this work we analyzed two different schemes to converge the correlation energy, one based on a more traditional complete basis set extrapolation and one that converges energy differences by accounting for the size-consistency property. These two approaches have been systematically tested on the A24 test set, for six points on the potential-energy surface of the methane-formaldehyde complex, and for reaction energies involving the breaking and formation of covalent bonds. While both methods converge to similar results at similar rates, the computation of size-consistent energy differences has the advantage of not relying on the choice of a specific extrapolation model.

The multifacet graphically contracted function method. I. Formulation and implementation

International Nuclear Information System (INIS)

Shepard, Ron; Brozell, Scott R.; Gidofalvi, Gergely

2014-01-01

The basic formulation for the multifacet generalization of the graphically contracted function (MFGCF) electronic structure method is presented. The analysis includes the discussion of linear dependency and redundancy of the arc factor parameters, the computation of reduced density matrices, Hamiltonian matrix construction, spin-density matrix construction, the computation of optimization gradients for single-state and state-averaged calculations, graphical wave function analysis, and the efficient computation of configuration state function and Slater determinant expansion coefficients. Timings are given for Hamiltonian matrix element and analytic optimization gradient computations for a range of model problems for full-CI Shavitt graphs, and it is observed that both the energy and the gradient computation scale as O(N 2 n 4 ) for N electrons and n orbitals. The important arithmetic operations are within dense matrix-matrix product computational kernels, resulting in a computationally efficient procedure. An initial implementation of the method is used to present applications to several challenging chemical systems, including N 2 dissociation, cubic H 8 dissociation, the symmetric dissociation of H 2 O, and the insertion of Be into H 2 . The results are compared to the exact full-CI values and also to those of the previous single-facet GCF expansion form

The multifacet graphically contracted function method. I. Formulation and implementation

Science.gov (United States)

Shepard, Ron; Gidofalvi, Gergely; Brozell, Scott R.

2014-08-01

The basic formulation for the multifacet generalization of the graphically contracted function (MFGCF) electronic structure method is presented. The analysis includes the discussion of linear dependency and redundancy of the arc factor parameters, the computation of reduced density matrices, Hamiltonian matrix construction, spin-density matrix construction, the computation of optimization gradients for single-state and state-averaged calculations, graphical wave function analysis, and the efficient computation of configuration state function and Slater determinant expansion coefficients. Timings are given for Hamiltonian matrix element and analytic optimization gradient computations for a range of model problems for full-CI Shavitt graphs, and it is observed that both the energy and the gradient computation scale as O(N2n4) for N electrons and n orbitals. The important arithmetic operations are within dense matrix-matrix product computational kernels, resulting in a computationally efficient procedure. An initial implementation of the method is used to present applications to several challenging chemical systems, including N2 dissociation, cubic H8 dissociation, the symmetric dissociation of H2O, and the insertion of Be into H2. The results are compared to the exact full-CI values and also to those of the previous single-facet GCF expansion form.

Mechanisms of the electron density depletion in the SAR arc region

Directory of Open Access Journals (Sweden)

A. V. Pavlov

1996-02-01

Full Text Available This study compares the measurements of electron density and temperature and the integral airglow intensity at 630 nm in the SAR arc region and slightly south of this (obtained by the Isis 2 spacecraft during the 18 December 1971 magnetic storm, with the model results obtained using the time dependent one-dimensional mathematical model of the Earth\\'s ionosphere and plasmasphere. The explicit expression in the third Enskog approximation for the electron thermal conductivity coefficient in the multicomponent mixture of ionized gases and a simplified calculation method for this coefficient presents an opportunity to calculate more exactly the electron temperature and density and 630 nm emission within SAR arc region are used in the model. Collisions between N2 and hot thermal electrons in the SAR arc region produce vibrationally excited nitrogen molecules. It appears that the loss rate of O+(4S due to reactions with the vibrationally excited nitrogen is enough to explain electron density depression by a factor of two at F-region heights and the topside ionosphere density variations within the SAR arc if the erosion of plasma within geomagnetic field tubes, during the main phase of the geomagnetic storm and subsequent filling of geomagnetic tubes during the recovery phase, are considered. To explain the disagreement by a factor 1.5 between the observed and modeled SAR arc electron densities an additional plasma drift velocity ~–30 m s–1 in the ion continuity equations is needed during the recovery phase. This additional plasma drift velocity is likely caused by the transition from convecting to corotating flux tubes on the equatorward wall of the trough. The electron densities and temperatures and 630 nm integral intensity at the SAR arc and slightly south of this region as measured for the 18 December 1971 magnetic storm were correctly described by the model without perpendicular electric fields. Within this model framework the effect of the

Two-Electron Transfer Pathways.

Science.gov (United States)

Lin, Jiaxing; Balamurugan, D; Zhang, Peng; Skourtis, Spiros S; Beratan, David N

2015-06-18

The frontiers of electron-transfer chemistry demand that we develop theoretical frameworks to describe the delivery of multiple electrons, atoms, and ions in molecular systems. When electrons move over long distances through high barriers, where the probability for thermal population of oxidized or reduced bridge-localized states is very small, the electrons will tunnel from the donor (D) to acceptor (A), facilitated by bridge-mediated superexchange interactions. If the stable donor and acceptor redox states on D and A differ by two electrons, it is possible that the electrons will propagate coherently from D to A. While structure-function relations for single-electron superexchange in molecules are well established, strategies to manipulate the coherent flow of multiple electrons are largely unknown. In contrast to one-electron superexchange, two-electron superexchange involves both one- and two-electron virtual intermediate states, the number of virtual intermediates increases very rapidly with system size, and multiple classes of pathways interfere with one another. In the study described here, we developed simple superexchange models for two-electron transfer. We explored how the bridge structure and energetics influence multielectron superexchange, and we compared two-electron superexchange interactions to single-electron superexchange. Multielectron superexchange introduces interference between singly and doubly oxidized (or reduced) bridge virtual states, so that even simple linear donor-bridge-acceptor systems have pathway topologies that resemble those seen for one-electron superexchange through bridges with multiple parallel pathways. The simple model systems studied here exhibit a richness that is amenable to experimental exploration by manipulating the multiple pathways, pathway crosstalk, and changes in the number of donor and acceptor species. The features that emerge from these studies may assist in developing new strategies to deliver multiple

Convergent J-matrix calculation of the Poet-Temkin model of electron-hydrogen scattering

International Nuclear Information System (INIS)

Konovalov, D.A.; McCarthy, I.E.

1994-01-01

It is shown that the Poet-Temkin model of electron-hydrogen scattering could be solved to any required accuracy using the J-matrix method. The convergence in the basis size is achieved to an accuracy of better than 2% with the inclusion of 37 basis L 2 functions. Previously observed pseudoresonances in the J-matrix calculation naturally disappear with an increase in basis size. No averaging technique is necessary to smooth the convergent J-matrix results. (Author)

Matrix density effects on the mechanical properties of SiC fiber-reinforced silicon nitride matrix properties

Science.gov (United States)

Bhatt, Ramakrishna T.; Kiser, Lames D.

1990-01-01

The room temperature mechanical properties were measured for SiC fiber reinforced reaction-bonded silicon nitride composites (SiC/RBSN) of different densities. The composites consisted of approx. 30 vol percent uniaxially aligned 142 micron diameter SiC fibers (Textron SCS-6) in a reaction-bonded Si3N4 matrix. The composite density was varied by changing the consolidation pressure during RBSN processing and by hot isostatically pressing the SiC/RBSN composites. Results indicate that as the consolidation pressure was increased from 27 to 138 MPa, the average pore size of the nitrided composites decreased from 0.04 to 0.02 microns and the composite density increased from 2.07 to 2.45 gm/cc. Nonetheless, these improvements resulted in only small increases in the first matrix cracking stress, primary elastic modulus, and ultimate tensile strength values of the composites. In contrast, HIP consolidation of SiC/RBSN resulted in a fully dense material whose first matrix cracking stress and elastic modulus were approx. 15 and 50 percent higher, respectively, and ultimate tensile strength values were approx. 40 percent lower than those for unHIPed SiC/RBSN composites. The modulus behavior for all specimens can be explained by simple rule-of-mixture theory. Also, the loss in ultimate strength for the HIPed composites appears to be related to a degradation in fiber strength at the HIP temperature. However, the density effect on matrix fracture strength was much less than would be expected based on typical monolithic Si3N4 behavior, suggesting that composite theory is indeed operating. Possible practical implications of these observations are discussed.

Electron density in surface barrier discharge emerging at argon/water interface: quantification for streamers and leaders

Science.gov (United States)

Cvetanović, Nikola; Galmiz, Oleksandr; Synek, Petr; Zemánek, Miroslav; Brablec, Antonín; Hoder, Tomáš

2018-02-01

Optical emission spectroscopy, fast intensified CCD imaging and electrical measurements were applied to investigate the basic plasma parameters of surface barrier discharge emerging from a conductive water electrode. The discharge was generated at the triple-line interface of atmospheric pressure argon gas and conductive water solution at the fused silica dielectrics using a sinusoidal high-voltage waveform. The spectroscopic methods of atomic line broadening and molecular spectroscopy were used to determine the electron densities and the gas temperature in the active plasma. These parameters were obtained for both applied voltage polarities and resolved spatially. Two different spectral signatures were identified in the spatially resolved spectra resulting in electron densities differing by two orders of magnitude. It is shown that two discharge mechanisms take a place: the streamer and the leader one, with electron densities of 1014 and 1016 cm-3, respectively. This spectroscopic evidence is supported by the combined diagnostics of electrical current measurements and phase-resolved intensified CCD camera imaging.

Communication: Wigner functions in action-angle variables, Bohr-Sommerfeld quantization, the Heisenberg correspondence principle, and a symmetrical quasi-classical approach to the full electronic density matrix

International Nuclear Information System (INIS)

Miller, William H.; Cotton, Stephen J.

2016-01-01

It is pointed out that the classical phase space distribution in action-angle (a-a) variables obtained from a Wigner function depends on how the calculation is carried out: if one computes the standard Wigner function in Cartesian variables (p, x), and then replaces p and x by their expressions in terms of a-a variables, one obtains a different result than if the Wigner function is computed directly in terms of the a-a variables. Furthermore, the latter procedure gives a result more consistent with classical and semiclassical theory—e.g., by incorporating the Bohr-Sommerfeld quantization condition (quantum states defined by integer values of the action variable) as well as the Heisenberg correspondence principle for matrix elements of an operator between such states—and has also been shown to be more accurate when applied to electronically non-adiabatic applications as implemented within the recently developed symmetrical quasi-classical (SQC) Meyer-Miller (MM) approach. Moreover, use of the Wigner function (obtained directly) in a-a variables shows how our standard SQC/MM approach can be used to obtain off-diagonal elements of the electronic density matrix by processing in a different way the same set of trajectories already used (in the SQC/MM methodology) to obtain the diagonal elements.

Communication: Wigner functions in action-angle variables, Bohr-Sommerfeld quantization, the Heisenberg correspondence principle, and a symmetrical quasi-classical approach to the full electronic density matrix.

Science.gov (United States)

Miller, William H; Cotton, Stephen J

2016-08-28

It is pointed out that the classical phase space distribution in action-angle (a-a) variables obtained from a Wigner function depends on how the calculation is carried out: if one computes the standard Wigner function in Cartesian variables (p, x), and then replaces p and x by their expressions in terms of a-a variables, one obtains a different result than if the Wigner function is computed directly in terms of the a-a variables. Furthermore, the latter procedure gives a result more consistent with classical and semiclassical theory-e.g., by incorporating the Bohr-Sommerfeld quantization condition (quantum states defined by integer values of the action variable) as well as the Heisenberg correspondence principle for matrix elements of an operator between such states-and has also been shown to be more accurate when applied to electronically non-adiabatic applications as implemented within the recently developed symmetrical quasi-classical (SQC) Meyer-Miller (MM) approach. Moreover, use of the Wigner function (obtained directly) in a-a variables shows how our standard SQC/MM approach can be used to obtain off-diagonal elements of the electronic density matrix by processing in a different way the same set of trajectories already used (in the SQC/MM methodology) to obtain the diagonal elements.

Communication: Wigner functions in action-angle variables, Bohr-Sommerfeld quantization, the Heisenberg correspondence principle, and a symmetrical quasi-classical approach to the full electronic density matrix

Energy Technology Data Exchange (ETDEWEB)

Miller, William H., E-mail: millerwh@berkeley.edu; Cotton, Stephen J., E-mail: StephenJCotton47@gmail.com [Department of Chemistry and Kenneth S. Pitzer Center for Theoretical Chemistry, University of California, and Chemical Sciences Division, Lawrence Berkeley National Laboratory, Berkeley, California 94720 (United States)

2016-08-28

It is pointed out that the classical phase space distribution in action-angle (a-a) variables obtained from a Wigner function depends on how the calculation is carried out: if one computes the standard Wigner function in Cartesian variables (p, x), and then replaces p and x by their expressions in terms of a-a variables, one obtains a different result than if the Wigner function is computed directly in terms of the a-a variables. Furthermore, the latter procedure gives a result more consistent with classical and semiclassical theory—e.g., by incorporating the Bohr-Sommerfeld quantization condition (quantum states defined by integer values of the action variable) as well as the Heisenberg correspondence principle for matrix elements of an operator between such states—and has also been shown to be more accurate when applied to electronically non-adiabatic applications as implemented within the recently developed symmetrical quasi-classical (SQC) Meyer-Miller (MM) approach. Moreover, use of the Wigner function (obtained directly) in a-a variables shows how our standard SQC/MM approach can be used to obtain off-diagonal elements of the electronic density matrix by processing in a different way the same set of trajectories already used (in the SQC/MM methodology) to obtain the diagonal elements.

Density matrix-based time-dependent configuration interaction approach to ultrafast spin-flip dynamics

Science.gov (United States)

Wang, Huihui; Bokarev, Sergey I.; Aziz, Saadullah G.; Kühn, Oliver

2017-08-01

Recent developments in attosecond spectroscopy yield access to the correlated motion of electrons on their intrinsic timescales. Spin-flip dynamics is usually considered in the context of valence electronic states, where spin-orbit coupling is weak and processes related to the electron spin are usually driven by nuclear motion. However, for core-excited states, where the core-hole has a nonzero angular momentum, spin-orbit coupling is strong enough to drive spin-flips on a much shorter timescale. Using density matrix-based time-dependent restricted active space configuration interaction including spin-orbit coupling, we address an unprecedentedly short spin-crossover for the example of L-edge (2p→3d) excited states of a prototypical Fe(II) complex. This process occurs on a timescale, which is faster than that of Auger decay (∼4 fs) treated here explicitly. Modest variations of carrier frequency and pulse duration can lead to substantial changes in the spin-state yield, suggesting its control by soft X-ray light.

Non-zero temperature two-mode squeezing for time-dependent two-level systems

International Nuclear Information System (INIS)

Aliaga, J.; Gruver, J.L.; Proto, A.N.; Cerdeira, H.A.

1994-01-01