Science.gov

Sample records for quasiparticle self-consistent gw

  1. Quasi-Particle Self-Consistent GW for Molecules.

    PubMed

    Kaplan, F; Harding, M E; Seiler, C; Weigend, F; Evers, F; van Setten, M J

    2016-06-14

    We present the formalism and implementation of quasi-particle self-consistent GW (qsGW) and eigenvalue only quasi-particle self-consistent GW (evGW) adapted to standard quantum chemistry packages. Our implementation is benchmarked against high-level quantum chemistry computations (coupled-cluster theory) and experimental results using a representative set of molecules. Furthermore, we compare the qsGW approach for five molecules relevant for organic photovoltaics to self-consistent GW results (scGW) and analyze the effects of the self-consistency on the ground state density by comparing calculated dipole moments to their experimental values. We show that qsGW makes a significant improvement over conventional G0W0 and that partially self-consistent flavors (in particular evGW) can be excellent alternatives. PMID:27168352

  2. Finite temperature quasiparticle self-consistent GW approximation

    SciTech Connect

    Vanschilfgaarde, Mark; Leonard, Fran­cois; Desjarlais, Michael Paul; Kotani, Takao; Faleev, Sergey V

    2005-10-01

    We present a new ab initio method for electronic structure calculations of materials at finite temperature (FT) based on the all-electron quasiparticle self-consistent GW (QPscGW) approximation and Keldysh time-loop Green's function approach. We apply the method to Si, Ge, GaAs, InSb, and diamond and show that the band gaps of these materials universally decrease with temperature in contrast with the local density approximation (LDA) of density functional theory (DFT) where the band gaps universally increase. At temperatures of a few eV the difference between quasiparticle energies obtained in FT-QPscGW and FT-LDA approaches significantly reduces. This result suggests that existing simulations of very high temperature materials based on the FT-LDA are more justified then it might appear from well-known LDA band gap errors at zero-temperature.

  3. Electronic structure from relativistic quasiparticle self-consistent GW calculations

    NASA Astrophysics Data System (ADS)

    Blügel, Stefan

    Most theoretical studies of topological insulators (TIs) are based on tight-binding descriptions and density functional theory (DFT). But recently, many-body calculations within the GW approximation attract much attention in the study of these materials. We present an implementation of the quasiparticle self-consistent (QS) GW method where the spin-orbit coupling (SOC) is fully taken into account in each iteration rather than added a posteriori. Within the all-electron FLAPW formalism, we show DFT, one-shot GW , and QS GW calculations for several, well-known TIs. We present a comparison of the calculations to photoemission spectroscopy and show that the GW corrected bands agree much better with experiment. For example, we show that Bi2Se3 is a direct gap semiconductor, in contrast to what was believed for many years by interpreting experimental results on the basis of DFT and that small strains in Bi can lead to a semimetal-to-semiconductor or trivial-to-topological transitions. Quasiparticle calculations for low-dimensional systems are still very demanding. In order to study the topological surface states with an approach based on GW , we use Wannier functions to construct a Hamiltonian that reproduces the many-body band structure of the bulk, and that is used to construct a slab Hamiltonian. With this approach, we discuss the effect of quasiparticle corrections on the surface states of TIs and on the interaction between bulk and surface states Work was funded by the Virtual Institute for Topological Insulators of the Helmholtz Association and carried out in collaboration with Irene Aguilera, Gustav Bihlmayer, and Christoph Friedrich.

  4. Quasiparticle self-consistent GW study of cuprates: electronic structure, model parameters, and the two-band theory for Tc

    PubMed Central

    Jang, Seung Woo; Kotani, Takao; Kino, Hiori; Kuroki, Kazuhiko; Han, Myung Joon

    2015-01-01

    Despite decades of progress, an understanding of unconventional superconductivity still remains elusive. An important open question is about the material dependence of the superconducting properties. Using the quasiparticle self-consistent GW method, we re-examine the electronic structure of copper oxide high-Tc materials. We show that QSGW captures several important features, distinctive from the conventional LDA results. The energy level splitting between and is significantly enlarged and the van Hove singularity point is lowered. The calculated results compare better than LDA with recent experimental results from resonant inelastic xray scattering and angle resolved photoemission experiments. This agreement with the experiments supports the previously suggested two-band theory for the material dependence of the superconducting transition temperature, Tc. PMID:26206417

  5. Accurate energy bands calculated by the hybrid quasiparticle self-consistent GW method implemented in the ecalj package

    NASA Astrophysics Data System (ADS)

    Deguchi, Daiki; Sato, Kazunori; Kino, Hiori; Kotani, Takao

    2016-05-01

    We have recently implemented a new version of the quasiparticle self-consistent GW (QSGW) method in the ecalj package released at http://github.com/tkotani/ecalj. Since the new version of the ecalj package is numerically stable and more accurate than the previous versions, we can perform calculations easily without being bothered with tuning input parameters. Here we examine its ability to describe energy band properties, e.g., band-gap energy, eigenvalues at special points, and effective mass, for a variety of semiconductors and insulators. We treat C, Si, Ge, Sn, SiC (in 2H, 3C, and 4H structures), (Al, Ga, In) × (N, P, As, Sb), (Zn, Cd, Mg) × (O, S, Se, Te), SiO2, HfO2, ZrO2, SrTiO3, PbS, PbTe, MnO, NiO, and HgO. We propose that a hybrid QSGW method, where we mix 80% of QSGW and 20% of LDA, gives universally good agreement with experiments for these materials.

  6. Electronic structure calculations of delafossite Cu-based transparent conducting oxides CuMO2 (M =B,Al,Ga,In) by quasiparticle self-consistent GW approximation and Tran-Blaha's modified Becke-Johnson exchange potential

    NASA Astrophysics Data System (ADS)

    Thatribud, Abdulmutta; Pengpan, Teparksorn

    2014-09-01

    In this work, band gaps of the delafossite Cu-based transparent conducting oxides CuMO2 (M =B,Al,Ga,In) are calculated by density functional theory (DFT) implemented with many-body perturbation theory (MBPT) based on quasiparticle self-consistent GW approximation (QPscGW) and with Tran-Blaha's modified Becke-Johnson functional (DFT-TB09). Their band gaps are explicitly improved from DFT within local density approximation (LDA). Their optical absorption spectra are also calculated by solving Bethe-Salpeter equation (BSE) that includes the electron-hole correlation effect; they show strong excitonic peaks.

  7. Self-consistent GW calculations for semiconductors and insulators

    NASA Astrophysics Data System (ADS)

    Shishkin, M.; Kresse, G.

    2007-06-01

    We present GW calculations for small and large gap systems comprising typical semiconductors (Si, SiC, GaAs, GaN, ZnO, ZnS, CdS, and AlP), small gap semiconductors (PbS, PbSe, and PbTe), insulators (C, BN, MgO, and LiF), and noble gas solids (Ar and Ne). It is shown that the G0W0 approximation always yields too small band gaps. To improve agreement with experiment, the eigenvalues in the Green’s function G (GW0) and in the Green’s function and the dielectric matrix (GW) are updated until self-consistency is reached. The first approximation leads to excellent agreement with experiment, whereas an update of the eigenvalues in G and W gives too large band gaps for virtually all materials. From a pragmatic point of view, the GW0 approximation thus seems to be an accurate and still reasonably fast method for predicting quasiparticle energies in simple sp -bonded systems. We furthermore observe that the band gaps in materials with shallow d states (GaAs, GaN, and ZnO) are systematically underestimated. We propose that an inaccurate description of the static dielectric properties of these materials is responsible for the underestimation of the band gaps in GW0 , which is itself a result of the incomplete cancellation of the Hartree self-energy within the d shell by local or gradient corrected density functionals.

  8. Full self-consistency versus quasiparticle self-consistency in diagrammatic approaches: Exactly solvable two-site Hubbard model

    SciTech Connect

    Kutepov, A. L.

    2015-07-22

    Self-consistent solutions of Hedin's equations (HE) for the two-site Hubbard model (HM) have been studied. They have been found for three-point vertices of increasing complexity (Γ = 1 (GW approximation), Γ₁ from the first-order perturbation theory, and the exact vertex ΓE). Comparison is made between the cases when an additional quasiparticle (QP) approximation for Green's functions is applied during the self-consistent iterative solving of HE and when QP approximation is not applied. Results obtained with the exact vertex are directly related to the present open question—which approximation is more advantageous for future implementations, GW + DMFT or QPGW + DMFT. It is shown that in a regime of strong correlations only the originally proposed GW + DMFT scheme is able to provide reliable results. Vertex corrections based on Perturbation Theory systematically improve the GW results when full self-consistency is applied. The application of QP self-consistency combined with PT vertex corrections shows similar problems to the case when the exact vertex is applied combined with QP sc. An analysis of Ward Identity violation is performed for all studied in this work's approximations and its relation to the general accuracy of the schemes used is provided.

  9. Full self-consistency versus quasiparticle self-consistency in diagrammatic approaches: Exactly solvable two-site Hubbard model

    DOE PAGESBeta

    Kutepov, A. L.

    2015-07-22

    Self-consistent solutions of Hedin's equations (HE) for the two-site Hubbard model (HM) have been studied. They have been found for three-point vertices of increasing complexity (Γ = 1 (GW approximation), Γ₁ from the first-order perturbation theory, and the exact vertex ΓE). Comparison is made between the cases when an additional quasiparticle (QP) approximation for Green's functions is applied during the self-consistent iterative solving of HE and when QP approximation is not applied. Results obtained with the exact vertex are directly related to the present open question—which approximation is more advantageous for future implementations, GW + DMFT or QPGW + DMFT.more » It is shown that in a regime of strong correlations only the originally proposed GW + DMFT scheme is able to provide reliable results. Vertex corrections based on Perturbation Theory systematically improve the GW results when full self-consistency is applied. The application of QP self-consistency combined with PT vertex corrections shows similar problems to the case when the exact vertex is applied combined with QP sc. An analysis of Ward Identity violation is performed for all studied in this work's approximations and its relation to the general accuracy of the schemes used is provided.« less

  10. Bilayer graphene spectral function in the random phase approximation and self-consistent GW approximation

    NASA Astrophysics Data System (ADS)

    Sabashvili, Andro; Östlund, Stellan; Granath, Mats

    2013-08-01

    We calculate the single-particle spectral function for doped bilayer graphene in the low energy limit, described by two parabolic bands with zero band gap and long range Coulomb interaction. Calculations are done using thermal Green's functions in both the random phase approximation (RPA) and the fully self-consistent GW approximation. Consistent with previous studies RPA yields a spectral function which, apart from the Landau quasiparticle peaks, shows additional coherent features interpreted as plasmarons, i.e., composite electron-plasmon excitations. In the GW approximation the plasmaron becomes incoherent and peaks are replaced by much broader features. The deviation of the quasiparticle weight and mass renormalization from their noninteracting values is small which indicates that bilayer graphene is a weakly interacting system. The electron energy loss function, Im[-ɛq-1(ω)] shows a sharp plasmon mode in RPA which in the GW approximation becomes less coherent and thus consistent with the weaker plasmaron features in the corresponding single-particle spectral function.

  11. Effects of self-consistency and plasmon-pole models on GW calculations for closed-shell molecules

    NASA Astrophysics Data System (ADS)

    Lischner, Johannes; Sharifzadeh, Sahar; Deslippe, Jack; Neaton, Jeffrey B.; Louie, Steven G.

    2014-09-01

    We present theoretical calculations of quasiparticle energies in closed-shell molecules using the GW method. We compare three different approaches: a full-frequency G0W0 (FF-G0W0) method with density functional theory (DFT-PBE) used as a starting mean field; a full-frequency GW0 (FF-GW0) method where the interacting Green's function is approximated by replacing the DFT energies with self-consistent quasiparticle energies or Hartree-Fock energies; and a G0W0 method with a Hybertsen-Louie generalized plasmon-pole model (HL GPP-G0W0). While the latter two methods lead to good agreement with experimental ionization potentials and electron affinities for methane, ozone, and beryllium oxide molecules, FF-G0W0 results can differ by more than one electron volt from experiment. We trace this failure of the FF-G0W0 method to the occurrence of incorrect self-energy poles describing shake-up processes in the vicinity of the quasiparticle energies.

  12. Thermodynamic properties of hot nuclei within the self-consistent quasiparticle random-phase approximation

    SciTech Connect

    Hung, N. Quang; Dang, N. Dinh

    2010-10-15

    The thermodynamic properties of hot nuclei are described within the canonical and microcanonical ensemble approaches. These approaches are derived based on the solutions of the BCS and self-consistent quasiparticle random-phase approximation at zero temperature embedded into the canonical and microcanonical ensembles. The results obtained agree well with the recent data extracted from experimental level densities by the Oslo group for {sup 94}Mo, {sup 98}Mo, {sup 162}Dy, and {sup 172}Yb nuclei.

  13. Pairing within the self-consistent quasiparticle random-phase approximation at finite temperature

    SciTech Connect

    Dang, N. Dinh; Hung, N. Quang

    2008-06-15

    An approach to pairing in finite nuclei at nonzero temperature is proposed, which incorporates the effects due to the quasiparticle-number fluctuation (QNF) around Bardeen-Cooper-Schrieffer (BCS) mean field and dynamic coupling to quasiparticle-pair vibrations within the self-consistent quasiparticle random-phase approximation (SCQRPA). The numerical calculations of pairing gap, total energy, and heat capacity were carried out within a doubly folded multilevel model as well as realistic nuclei {sup 56}Fe and {sup 120}Sn. The results obtained show that, under the effect of QNF, in the region of moderate and strong couplings, the sharp transition between the superconducting and normal phases is smoothed out, resulting in a thermal pairing gap, which does not collapse at the BCS critical temperature, but has a tail, which extends to high temperature. The dynamic coupling of quasiparticles to SCQRPA vibrations significantly improves the agreement with the results of exact calculations and those obtained within the finite-temperature quantal Monte Carlo method for the total energy and heat capacity. It also causes a deviation of the quasiparticle occupation numbers from the Fermi-Dirac distributions for free fermions.

  14. Linear response of light deformed nuclei investigated by self-consistent quasiparticle random-phase approximation

    SciTech Connect

    Losa, C.; Doessing, T.; Pastore, A.; Vigezzi, E.; Broglia, R. A.

    2010-06-15

    We present a calculation of the properties of vibrational states in deformed, axially-symmetric even-even nuclei, within the framework of a fully self-consistent quasiparticle random phase approximation (QRPA). The same Skyrme energy density and density-dependent pairing functionals are used to calculate the mean field and the residual interaction in the particle-hole and particle-particle channels. We have tested our software in the case of spherical nuclei against fully self-consistent calculations published in the literature, finding excellent agreement. We investigate the consequences of neglecting the spin-orbit and Coulomb residual interactions in QRPA. Furthermore we discuss the improvement obtained in the QRPA result associated with the removal of spurious modes. Isoscalar and isovector responses in the deformed {sup 24-26}Mg, {sup 34}Mg isotopes are presented and compared to experimental findings.

  15. Maximally-localized Wannier functions for GW quasiparticles

    NASA Astrophysics Data System (ADS)

    Hamann, D. R.; Vanderbilt, David

    2009-03-01

    Recent efforts carrying the GW many-body approximation to self-consistency have given improved electronic structure results.^1 Maximally-localized Wannier functions^2 formed from the quasiparticle wave functions^3 provide an efficient and highly accurate basis for interpolating the SCGW bands from a coarse Brillouin-zone mesh to symmetry lines. Since the MLWF's correspond to chemists' bond orbitals, they potentially also provide insight into the qualitative effects of the improved treatment of correlations in SCGW compared to LDA. We report results on SrTiO3, solid Ar, and molecular CO. Band interpolation is accurate and effective for both solids. Small shifts in the degree of hybridization can be visualized for some of the SrTiO3 and CO MLWF's. In Ar, individual conduction-band Bloch functions were found to have large differences between LDA and SCGW.^1 However, a manifold of 9 d and spd-hybrid MLWF's which proved to be the minimum necessary for the lower conduction bands showed minimal differences in the two cases. A fully-functional interface to the WANNIER90 library within the SCGW-capable ABINIT code has been implemented and will be publicly available in the near future. 1. F. Bruneval et al., Phys. Rev. B 74, 045102 (2006). 2. N. Marzari and D. Vanderbilt, Phys. Rev. B 56, 12 847 (1997). 3. M. van Schilfgaarde et al., Phys. Rev. Lett. 96, 226402 (2006).

  16. A self-consistent GW approach to the van der Waals potential for a helium dimer.

    PubMed

    Shoji, Toru; Kuwahara, Riichi; Ono, Shota; Ohno, Kaoru

    2016-09-21

    van der Waals interaction between two helium (He) atoms is studied by calculating the total energy as a function of the He-He distance within the self-consistent GW approximation, which is expected to behave correctly in the long wavelength limit. In the Born-Oppenheimer (BO) approximation, the pair potential curve has its minimum value at 2.87 Å, which is somewhat larger than the local density approximation result, 2.40 Å, and is closer to previous quantum chemistry results. The expectation value for the interatomic distance, calculated by solving the Schrödinger equation for the two nuclei problem using the BO potential energy curve, is 30 Å, which is smaller but of the same order as previous experimental and theoretical results. PMID:27538378

  17. Quasiparticle electronic structure of Bi2Se3 via the sc-COHSEX+GW approach

    NASA Astrophysics Data System (ADS)

    Barker, Bradford A.; Deslippe, Jack; Yazyev, Oleg; Louie, Steven G.

    We present ab initio calculations of the quasiparticle electronic band structure of three-dimensional topological insulator material Bi2Se3 using the full spinor GW approach. The mean-field is initially computed at the DFT level in the local density approximation (LDA) using fully-relativistic pseudopotentials. We then improve the mean-field electronic structure by solving Dyson's equation in the static COHSEX approximation, self-consistently updating the eigenvalues, eigenvectors, and dielectric screening. After a few iterations, we then perform a GW calculation to determine the quasiparticle energies. We compare our calculated results to experimental values of the band gaps and effective masses. This work was supported by NSF Grant No. DMR15-1508412 and U.S. DOE under Contract No. DE-AC02-05CH11231. Computational resources have been provided by DOE at LBNL's NERSC facility and the NSF through XSEDE resources at NICS.

  18. Self-consistent quasiparticle random-phase approximation for a multilevel pairing model

    SciTech Connect

    Hung, N. Quang; Dang, N. Dinh

    2007-11-15

    Particle-number projection within the Lipkin-Nogami (LN) method is applied to the self-consistent quasiparticle random-phase approximation (SCQRPA), which is tested in an exactly solvable multilevel pairing model. The SCQRPA equations are numerically solved to find the energies of the ground and excited states at various numbers {omega} of doubly degenerate equidistant levels. The use of the LN method allows one to avoid the collapse of the BCS (QRPA) to obtain the energies of the ground and excited states as smooth functions of the interaction parameter G. The comparison between results given by different approximations such as the SCRPA, QRPA, LNQRPA, SCQRPA, and LNSCQRPA is carried out. Although the use of the LN method significantly improves the agreement with the exact results in the intermediate coupling region, we found that in the strong coupling region the SCQRPA results are closest to the exact ones.

  19. Off-Diagonal Self-Energy Terms and Partially Self-Consistency in GW Calculations for Single Molecules: Efficient Implementation and Quantitative Effects on Ionization Potentials.

    PubMed

    Kaplan, F; Weigend, F; Evers, F; van Setten, M J

    2015-11-10

    The GW method in its most widespread variant takes, as an input, Kohn-Sham (KS) single particle energies and single particle states and yields results for the single-particle excitation energies that are significantly improved over the bare KS estimates. Fundamental shortcomings of density functional theory (DFT) when applied to excitation energies as well as artifacts introduced by approximate exchange-correlation (XC) functionals are thus reduced. At its heart lies the quasi-particle (qp) equation, whose solution yields the corrected excitation energies and qp-wave functions. We propose an efficient approximation scheme to treat this equation based on second-order perturbation theory and self-consistent iteration schemes. We thus avoid solving (large) eigenvalue problems at the expense of a residual error that is comparable to the intrinsic uncertainty of the GW truncation scheme and is, in this sense, insignificant. PMID:26574312

  20. Screened hybrid and self-consistent GW calculations of cadmium/magnesium indium sulfide materials

    NASA Astrophysics Data System (ADS)

    Lucero, Melissa J.; Aguilera, Irene; Diaconu, Cristian V.; Palacios, Pablo; Wahnón, Perla; Scuseria, Gustavo E.

    2011-05-01

    The cadmium and magnesium indium sulfides are medium-gap semiconductors demonstrating a propensity to form intermediate band materials when doped with transition metals. The inherent structural diversity exhibited by M+2In2S4 thiospinels and related AB2X4 compounds often precludes definitive experimental determination of the band-gap width and type of transition. Employing a series of traditional semilocal functionals (e.g., the local spin density approximation; the Perdew, Burke, and Enzerhof functional; and the Tao, Perdew, Staroverov, and Scuseria functional) the screened hybrid of Heyd, Scuseria, and Ernzerhof (HSE), band gaps, projected densities of states, and band structures are calculated for the normal, full inverse, and intermediate configurations of [Cd/Mg]8In16S32. Band structures and band gaps are also obtained via self-consistent many-body methods, using the static Coulomb-hole and screened exchange approximation to GW as a starting point for perturbative G0W0 calculations. Comparison to experiment indicates that HSE provides an accurate, computationally efficient, and relatively rapid means for predicting band-gap properties in spinel-type photovoltaic materials.

  1. Quasiparticle electronic structure of bismuth telluride in the GW approximation

    NASA Astrophysics Data System (ADS)

    Kioupakis, Emmanouil; Tiago, Murilo L.; Louie, Steven G.

    2010-12-01

    The quasiparticle band structure of bismuth telluride (Bi2Te3) , an important thermoelectric material that exhibits topologically insulating surface states, is calculated from first principles in the GW approximation. The quasiparticle energies are evaluated in fine detail in the first Brillouin zone using a Wannier-function interpolation method, allowing the accurate determination of the location of the band extrema (which is in the mirror plane) as well as the values of the quasiparticle band gap (0.17 eV) and effective-mass tensors. Spin-orbit interaction effects were included. The valence band exhibits two distinct maxima in the mirror plane that differ by just 1 meV, giving rise to one direct and one indirect band gap of very similar magnitude. The effective-mass tensors are in reasonable agreement with experiment. The Wannier interpolation coefficients can be used for the tight-binding parametrization of the band structure. Our work elucidates the electronic structure of Bi2Te3 and sheds light on its exceptional thermoelectric and topologically insulating properties.

  2. Quasiparticle electronic structure of bismuth telluride in the GW approximation

    SciTech Connect

    Kioupakis, Emmanouil; Tiago, Murilo L; Louie, Steven G.

    2010-01-01

    The quasiparticle band structure of bismuth telluride Bi2Te3 , an important thermoelectric material that exhibits topologically insulating surface states, is calculated from first principles in the GW approximation. The quasiparticle energies are evaluated in fine detail in the first Brillouin zone using a Wannier-function interpo- lation method, allowing the accurate determination of the location of the band extrema which is in the mirror plane as well as the values of the quasiparticle band gap 0.17 eV and effective-mass tensors. Spin-orbit interaction effects were included. The valence band exhibits two distinct maxima in the mirror plane that differ by just 1 meV, giving rise to one direct and one indirect band gap of very similar magnitude. The effective- mass tensors are in reasonable agreement with experiment. The Wannier interpolation coefficients can be used for the tight-binding parametrization of the band structure. Our work elucidates the electronic structure of Bi2Te3 and sheds light on its exceptional thermoelectric and topologically insulating properties.

  3. Quasiparticle and Optical Excitations in Solid Ne and Ar: GW and BSE Approximations

    SciTech Connect

    Patterson, Charles H.; Galamic-Mulaomerovic, S.

    2007-12-26

    The GW approximation and the Bethe-Salpeter equation (BSE) have been used to calculate quasiparticle and optical excitations in solid Ne and Ar. Absolute positions of quasiparticle and quasihole energies are found to be in very good agreement with experimental values. Binding energies of Frenkel excitons for these systems calculated using the BSE are also in good agreement with experiment. Splitting of excitons into longitudinal and transverse modes is calculated and found to be approximately twice the experimentally measured value.

  4. Screening and nonlocal correlations in the extended Hubbard model from self-consistent combined GW and dynamical mean field theory

    NASA Astrophysics Data System (ADS)

    Ayral, Thomas; Biermann, Silke; Werner, Philipp

    2013-03-01

    We describe a recent implementation of the combined GW and dynamical mean field method (GW+DMFT) for the two-dimensional Hubbard model with onsite and nearest-neighbor repulsion. We clarify the relation of the GW+DMFT scheme to alternative approaches in the literature, and discuss the corresponding approximations to the free-energy functional of the model. Furthermore, we describe a numerically exact technique for the solution of the GW+DMFT equations, namely, the hybridization expansion continuous-time algorithm for impurity models with retarded interactions. We compute the low-temperature phase diagram of the half-filled extended Hubbard model, addressing the metal-insulator transition at small intersite interactions and the transition to a charge-ordered state for stronger intersite repulsions. GW+DMFT introduces a nontrivial momentum dependence into the many-body self-energy and polarization. We find that the charge fluctuations included in the present approach have a larger impact on the latter than on the former. Finally, within the GW+DMFT framework, as in extended DMFT, the intersite repulsion translates into a frequency dependence of the local effective interaction. We analyze this dependence and show how it affects the local spectral function.

  5. Comparing quasiparticle GW+DMFT and LDA+DMFT for the test bed material SrVO3

    NASA Astrophysics Data System (ADS)

    Taranto, C.; Kaltak, M.; Parragh, N.; Sangiovanni, G.; Kresse, G.; Toschi, A.; Held, K.

    2013-10-01

    We have implemented the quasiparticle GW+dynamical mean field theory (DMFT) approach in the Vienna ab initio simulation package. To this end, a quasiparticle Hermitization of the G0W0 self-energy a lá Kotani-Schilfgaarde is employed, and the interaction values are obtained from the locally unscreened random phase approximation (RPA) using a projection onto Wannier orbitals. We compare quasiparticle GW+DMFT and local density approximation (LDA)+DMFT against each other and against experiment for SrVO3. We observe a partial compensation of stronger electronic correlations due to the reduced GW bandwidth and weaker correlations due to a larger screening of the RPA interaction, so that the obtained spectra are quite similar and agree well with experiment. Noteworthy, the quasiparticle GW+DMFT better reproduces the position of the lower Hubbard side band.

  6. Self-consistent Calculation of the Quasi-particle Energy Spectrum of Sodium using the Correlated Hartree Fock Method

    NASA Astrophysics Data System (ADS)

    Ishihara, Takamitsu; Yamagami, Hiroshi; Yasuhara, Hiroshi

    2001-12-01

    Self-consistent band calculation of sodium is performed in the correlated Hartree Fock scheme proposed by Yasuhara and Takada [Phys. Rev. B 43 (1991) 7200], which contains information on the effective mass of the electron liquid in the form of a nonlocal spin-parallel potential, and the remaining information of the self-energy operator in the form of a local potential. The bandwidth of occupied states is somewhat increased under the influence of the non-local spin-parallel potential, compared with the free electron value. No significant difference can be found in the distortion of the Fermi surface between the present theory and the LDA.

  7. Self-consistent quasiparticle formulation of a multiphonon method and its application to the neutron-rich O20 nucleus

    NASA Astrophysics Data System (ADS)

    De Gregorio, G.; Knapp, F.; Lo Iudice, N.; Vesely, P.

    2016-04-01

    A Bogoliubov quasiparticle formulation of an equation-of-motion phonon method, suited for open-shell nuclei, is derived. Like its particle-hole version, it consists of deriving a set of equations of motions whose iterative solution generates an orthonormal basis of n -phonon states (n =0 ,1 ,2 ,... ), built of quasiparticle Tamm-Dancoff phonons, which simplifies the solution of the eigenvalue problem. The method is applied to the open-shell neutron-rich O20 for illustrative purposes. A Hartree-Fock-Bogoliubov canonical basis, derived from an intrinsic two-body optimized chiral Hamiltonian, is used to derive and solve the eigenvalue equations in a space encompassing a truncated two-phonon basis. The spurious admixtures induced by the violation of the particle number and the center-of-mass motion are eliminated to a large extent by a Gram-Schmidt orthogonalization procedure. The calculation takes into account the Pauli principle, is self-consistent, and is parameter free except for the energy cutoff used to truncate the two-phonon basis, which induces an increasing depression of the ground state through its strong coupling to the quasiparticle vacuum. Such a cutoff is fixed so as to reproduce the first 1- level. The two-phonon states are shown to enhance the level density of the low-energy spectrum, consistently with the data, and to induce a fragmentation of the E 1 strength which, while accounting for the very low E 1 transitions, is not sufficient to reproduce the experimental cross section in the intermediate energy region. This and other discrepancies suggest the need of including the three-phonon states. These are also expected to offset the action of the two phonons on the quasiparticle vacuum and, therefore, free the calculation from any parameter.

  8. Quasiparticle and Optical Properties of Mono- and Bi-layer SnS2: A First-Principles GW and GW +BSE Study

    NASA Astrophysics Data System (ADS)

    Wu, Meng; Qiu, Diana; Louie, Steven G.

    2015-03-01

    Unlike most semiconducting transition metal dichalcogenides, SnS2, another layered metal dichalcogenide, is calculated within density functional theory to be an indirect bandgap semiconductor in both its bulk and monolayer forms. Experimental characterization of mono- and bi-layer SnS2 has been performed, but the details of its quasiparticle and excitonic properties remain unclear. Thus, we employ ab initio GW and GW +BSE calculations to study the quasiparticle band structure and optical absorption spectrum, respectively, of mono- and bi-layer SnS2 with spin-orbit coupling included throughout the calculations. We further investigate the character of excitonic states contributing to the optical spectrum. This work was supported by NSF Grant No. DMR10-1006184 and the U.S. Department of Energy under Contract No. DE-AC02-05CH11231. Computational resources have been provided by DOE at Lawrence Berkeley National Laboratory's NERSC facility.

  9. Substrate-induced renormalization of the quasiparticle and optical gaps in monolayer transition metal dichalcogenides from GW and GW-BSE calculations

    NASA Astrophysics Data System (ADS)

    da Jornada, Felipe H.; Ong, Chin Shen; Qiu, Diana Y.; Louie, Steven G.

    There has been a considerable effort to experimentally characterize the electronic and optical properties of novel atomically thin 2D semiconductors, such as mono- and few-layer transition metal dichalcogenides (TMDs). However, the role that different substrates play in these experiments still remains unclear. From a theoretical perspective, it is hard to include the substrate in an ab initio framework, while in experiments, it is often difficult to suspend these samples. Here, we present a new method to compute the substrate effect on the quasiparticle and optical properties of quasi-2D materials based on state-of-the-art ab initio GW and GW plus Bethe-Salpeter equation (GW-BSE) methods. We compute the effects of different metallic and semiconducting substrates, and show that the quasiparticle gap and exciton binding energy can be dramatically reduced even with semiconducting substrates. This work was supported by the National Science Foundation under Grant No. DMR15-1508412 and the DOE under Contract No. DE-AC02-05CH11231.

  10. Quasi-particle band structure of potassium-doped few-layer black phosphorus with GW approximation

    NASA Astrophysics Data System (ADS)

    Kim, Han-Gyu; Baik, Seung Su; Choi, Hyoung Joon

    We calculate the quasi-particle band structure of pristine and potassium-doped black phosphorus (BP) by using the GW approximation. We obtain band gaps of pristine bulk and few-layer BP and compare them with the result of the density functional calculations and experimental measurements. For potassium-doped cases, we calculate the electronic band structure of potassium-doped few-layer BPs with various doping densities. We obtain the critical doping density for the band-gap closing, and the energy-band dispersions when the band gap is inverted. We discuss Dirac semimetal properties of doped few-layer BPs obtained by the GW approximation. This work was supported by NRF of Korea (Grant No. 2011-0018306) and KISTI supercomputing center (Project No. KSC-2015-C3-039).

  11. On the combined use of GW approximation and cumulant expansion in the calculations of quasiparticle spectra: The paradigm of Si valence bands

    NASA Astrophysics Data System (ADS)

    Gumhalter, Branko; Kovač, Vjekoslav; Caruso, Fabio; Lambert, Henry; Giustino, Feliciano

    2016-07-01

    Since the earliest implementations of the various GW approximations and cumulant expansion in the calculations of quasiparticle propagators and spectra, several attempts have been made to combine the advantageous properties and results of these two theoretical approaches. While the GW-plus-cumulant approach has proven successful in interpreting photoemission spectroscopy data in solids, the formal connection between the two methods has not been investigated in detail. By introducing a general bijective integral representation of the cumulants, we can rigorously identify at which point these two approximations can be connected for the paradigmatic model of quasiparticle interaction with the dielectric response of the system that has been extensively exploited in recent interpretations of the satellite structures in photoelectron spectra. We establish a protocol for consistent practical implementation of the thus established GW +cumulant scheme and illustrate it by comprehensive state-of-the-art first-principles calculations of intrinsic angle-resolved photoemission spectra from Si valence bands.

  12. All-electron GW quasiparticle band structures of group 14 nitride compounds

    NASA Astrophysics Data System (ADS)

    Chu, Iek-Heng; Kozhevnikov, Anton; Schulthess, Thomas C.; Cheng, Hai-Ping

    2014-07-01

    We have investigated the group 14 nitrides (M3N4) in the spinel phase (γ-M3N4 with M = C, Si, Ge, and Sn) and β phase (β-M3N4 with M = Si, Ge, and Sn) using density functional theory with the local density approximation and the GW approximation. The Kohn-Sham energies of these systems have been first calculated within the framework of full-potential linearized augmented plane waves (LAPW) and then corrected using single-shot G0W0 calculations, which we have implemented in the modified version of the Elk full-potential LAPW code. Direct band gaps at the Γ point have been found for spinel-type nitrides γ-M3N4 with M = Si, Ge, and Sn. The corresponding GW-corrected band gaps agree with experiment. We have also found that the GW calculations with and without the plasmon-pole approximation give very similar results, even when the system contains semi-core d electrons. These spinel-type nitrides are novel materials for potential optoelectronic applications because of their direct and tunable band gaps.

  13. All-electron GW quasiparticle band structures of group 14 nitride compounds

    SciTech Connect

    Chu, Iek-Heng; Cheng, Hai-Ping; Kozhevnikov, Anton; Schulthess, Thomas C.

    2014-07-28

    We have investigated the group 14 nitrides (M{sub 3}N{sub 4}) in the spinel phase (γ-M{sub 3}N{sub 4} with M = C, Si, Ge, and Sn) and β phase (β-M{sub 3}N{sub 4} with M = Si, Ge, and Sn) using density functional theory with the local density approximation and the GW approximation. The Kohn-Sham energies of these systems have been first calculated within the framework of full-potential linearized augmented plane waves (LAPW) and then corrected using single-shot G{sub 0}W{sub 0} calculations, which we have implemented in the modified version of the Elk full-potential LAPW code. Direct band gaps at the Γ point have been found for spinel-type nitrides γ-M{sub 3}N{sub 4} with M = Si, Ge, and Sn. The corresponding GW-corrected band gaps agree with experiment. We have also found that the GW calculations with and without the plasmon-pole approximation give very similar results, even when the system contains semi-core d electrons. These spinel-type nitrides are novel materials for potential optoelectronic applications because of their direct and tunable band gaps.

  14. All-electron GW quasiparticle band structures of group 14 nitride compounds

    NASA Astrophysics Data System (ADS)

    Chu, Iek-Heng; Kozhenikov, Anton; Schulthess, Thomas; Cheng, Hai-Ping

    2014-03-01

    We have investigated the group 14 nitrides (M3N4) in both the spinel phase (with M =C, Si, Ge and Sn) and the beta phase (with M =Si, Ge and Sn) using density functional theory (DFT) with the local density approximation (LDA). The Kohn-Sham energies of these systems are first calculated within the framework of full-potential LAPW and then corrected using single-shot G0W0 calculations, which we have implemented in the Exciting-Plus code. Direct bands gap at the Γ point are found for all spinel-type nitrides. The calculated band gaps of Si3N4, Ge3N4 and Sn3N4 agree with experiment. We also find that for all systems studied, our GW calculations with and without the plasmon-pole approximation give very similar results, even when the system contains semi-core 3d electrons. These spinel-type nitrides are novel materials for potential optoelectronic applications. This work is supported by NSF/DMR-0804407 and DOE/BES-DE-FG02-02ER45995. Computations are performed using facilities at NERSC.

  15. GW quasiparticle energy study of ternary tetradymite Bi{sub 2}Te{sub 2}Se and Bi{sub 2}Te{sub 2}S thin films

    SciTech Connect

    Shuaibu, Alhassan; Rahman, Md. Mahmudur; Zainuddin, Hishamuddin; Talib, Zainal Abidin; Muhida, Rifki

    2015-04-24

    In this work, we have evaluated the quasiparticle energies of ternary tetradymite Bi{sub 2}Te{sub 2}Se and Bi{sub 2}Te{sub 2}S using first-principles calculation within the G{sub 0}W{sub 0} methods. We have also performed a broad convergence tests in order to investigate the quasiparticle corrections to the structural parameters and to the semi core d electrons in both of the compounds. For each case, we have calculated the many-body corrections within a one-shot GW method of the compounds. Our results have shown that for Bi{sub 2}Te{sub 2}Se the GW corrections increase the band gap to almost 10%, and for specific atomic positions, the band structure shows a close value to the experimental one. For Bi{sub 2}Te{sub 2}S, despite increase in the band gap due to the GW corrections, possibility of bulk resistivity that can be significant for photovoltaic applications was observed.

  16. BerkeleyGW: A massively parallel computer package for the calculation of the quasiparticle and optical properties of materials and nanostructures

    NASA Astrophysics Data System (ADS)

    Deslippe, Jack; Samsonidze, Georgy; Strubbe, David A.; Jain, Manish; Cohen, Marvin L.; Louie, Steven G.

    2012-06-01

    based on density-functional theory, fail to correctly capture this physics. Solution method: We construct and solve the Dyson's equation for the quasiparticle energies and wavefunctions within the GW approximation for the electron self-energy. We additionally construct and solve the Bethe-Salpeter equation for the correlated electron-hole (exciton) wavefunctions and excitation energies. Restrictions: The material size is limited in practice by the computational resources available. Materials with up to 500 atoms per periodic cell can be studied on large HPCs. Additional comments: The distribution file for this program is approximately 110 Mbytes and therefore is not delivered directly when download or E-mail is requested. Instead a html file giving details of how the program can be obtained is sent. Running time: 1-1000 minutes (depending greatly on system size and processor number).

  17. GW calculations on post-transition-metal oxides

    NASA Astrophysics Data System (ADS)

    Kang, Youngho; Kang, Gijae; Nahm, Ho-Hyun; Cho, Seong-Ho; Park, Young Soo; Han, Seungwu

    2014-04-01

    In order to establish the reliable GW scheme that can be consistently applied to post-transition-metal oxides (post-TMOs), we carry out comprehensive GW calculations on electronic structures of ZnO, Ga2O3,In2O3, and SnO2, the four representative post-TMOs. Various levels of self-consistency (G0W0, GW0, and QPGW0) and different starting functionals (GGA, GGA + U, and hybrid functional) are tested and their influence on the resulting electronic structure is closely analyzed. It is found that the GW0 scheme with GGA + U as the initial functional turns out to give the best agreement with experiment, implying that describing the position of metal-d level precisely in the ground state plays a critical role for the accurate dielectric property and quasiparticle band gap. Nevertheless, the computation on ZnO still suffers from the shallow Zn-d level and we propose a modified approach (GW0+Ud) that additionally considers an effective Hubbard U term during GW0 iterations and thereby significantly improves the band gap. It is also shown that a GGA + U-based GW0(+Ud) scheme produces an accurate energy gap of crystalline InGaZnO4, implying that this can serve as a standard scheme that can be applied to general structures of post-TMOs.

  18. First Principles Real-Space GW+BSE Calculations for Confined Systems

    NASA Astrophysics Data System (ADS)

    Hung, Linda; Ogut, Serdar; Souto, Jaime; Lee, Alex; Lena, Charles; Chelikowsky, James R.; Jornada, Felipe H. Da; Louie, Steven G.

    2014-03-01

    We investigate the performance of various levels of GW theories for electronic excitations as well as the resulting solutions of the Bethe-Salpeter-Equation (BSE) for optical excitations in a wide range of confined systems including atoms, ions, diatomic molecules, and organic molecules relevant for photovoltaic applications. Starting with solutions of the Kohn-Sham equations for ground state properties computed via the real-space ab initio pseudopotential code PARSEC, we perform the GW calculations in the space of single-particle transitions at various levels of theory, and compare the results with photoemission data. The levels of theory include such approximations as G0W0 with RPA screening, G0Wf that includes vertex corrections through the use of a dielectric screening within the time-dependent-local-density approximation (TDLDA), the GW0 , and the self-consistent GW. The resulting quasiparticle energies and wave functions from the GW calculations are used to solve the BSE for optical excitations, which are then compared with experiments and results from calculations performed within the TDLDA. The effects of the vertex corrections, self-consistency in GW, and core-valence partitioning are discussed. Supported by DOE Grant No. DE-SC0001853.

  19. GW-1000. GW Pharmaceuticals.

    PubMed

    Smith, Paul F

    2004-07-01

    GW Pharmaceuticals is developing GW-1000 (Sativex), a narrow ratio delta9-tetrahydrocannabinol:cannabidiol product for the potential treatment of multiple sclerosis, spinal cord injury, neurogenic pain and peripheral neuropathy. In March 2003, the company filed for approval for the treatment of MS with the UK Medicines Control Agency, and in May 2004, filed for new drug submission with Health Canada. PMID:15298072

  20. Improved quasiparticle wave functions and mean field for G0W0 calculations: Initialization with the COHSEX operator

    NASA Astrophysics Data System (ADS)

    Jain, Manish; Deslippe, Jack; Samsonidze, Georgy; Cohen, Marvin L.; Chelikowsky, James R.; Louie, Steven G.

    2014-09-01

    The GW approximation to the electron self-energy has become a standard method for ab initio calculation of excited-state properties of condensed-matter systems. In many calculations, the GW self-energy operator, Σ, is taken to be diagonal in the density functional theory (DFT) Kohn-Sham basis within the G0W0 scheme. However, there are known situations in which this diagonal G0W0 approximation starting from DFT is inadequate. We present two schemes to resolve such problems. The first, which we called sc -COHSEX+GW, involves construction of an improved mean field using the static limit of GW, known as COHSEX (Coulomb hole and screened exchange), which is significantly simpler to treat than GW. In this scheme, frequency-dependent self energy Σ (ω), is constructed and taken to be diagonal in the COHSEX orbitals after the system is solved self-consistently within this formalism. The second method is called off diagonal-COHSEX GW (od -COHSEX+GW). In this method, one does not self-consistently change the mean-field starting point but diagonalizes the COHSEX Hamiltonian within the Kohn-Sham basis to obtain quasiparticle wave functions and uses the resulting orbitals to construct the GW Σ in the diagonal form. We apply both methods to a molecular system, silane, and to two bulk systems, Si and Ge under pressure. For silane, both methods give good quasiparticle wave functions and energies. Both methods give good band gaps for bulk silicon and maintain good agreement with experiment. Further, the sc -COHSEX+GW method solves the qualitatively incorrect DFT mean-field starting point (having a band overlap) in bulk Ge under pressure.

  1. Self-consistent klystron simulations

    SciTech Connect

    Carlsten, B.E.; Tallerico, P.J.

    1985-01-01

    A numerical analysis of large-signal klystron behavior based on general wave-particle interaction theory is presented. The computer code presented is tailored for the minimum amount of complexity needed in klystron simulation. The code includes self-consistent electron motion, space-charge fields, and intermediate and output fields. It also includes use of time periodicity to simplify the problem, accurate representation of the space-charge fields, accurate representation of the cavity standing-wave fields, and a sophisticated particle-pushing routine. In the paper, examples are given that show the effects of cavity detunings, of varying the magnetic field profile, of electron beam asymmetries from the gun, and of variations in external load impedance. 4 refs., 7 figs.

  2. Self-consistent flattened isochrones

    NASA Astrophysics Data System (ADS)

    Binney, James

    2014-05-01

    We present a family of self-consistent axisymmetric stellar systems that have analytic distribution functions (DFs) of the form f(J), so they depend on three integrals of motion and have triaxial velocity ellipsoids. The models, which are generalizations of Hénon's isochrone sphere, have four dimensionless parameters, two determining the part of the DF that is even in Lz and two determining the odd part of the DF (which determines the azimuthal velocity distribution). Outside their cores, the velocity ellipsoids of all models tend to point to the model's centre, and we argue that this behaviour is generic, so near the symmetry axis of a flattened model, the long axis of the velocity ellipsoid is naturally aligned with the symmetry axis and not perpendicular to it as in many published dynamical models of well-studied galaxies. By varying one of the DF parameters, the intensity of rotation can be increased from zero up to a maximum value set by the requirement that the DF be non-negative. Since angle-action coordinates are easily computed for these models, they are ideally suited for perturbative treatments and stability analysis. They can also be used to choose initial conditions for an N-body model that starts in perfect equilibrium, and to model observations of early-type galaxies. The modelling technique introduced here is readily extended to different radial density profiles, more complex kinematics and multicomponent systems. A number of important technical issues surrounding the determination of the models' observable properties are explained in two appendices.

  3. Two-particle self-consistent approach to unconventional superconductivity

    NASA Astrophysics Data System (ADS)

    Otsuki, Junya

    2012-03-01

    A nonperturbative approach to unconventional superconductivity is developed based on the idea of two-particle self-consistent (TPSC) theory by Vilk and Tremblay. A sum rule, which the momentum-dependent pairing susceptibility satisfies, is derived. An effective pairing interaction between quasiparticles is determined so that the susceptibility should fulfill this exact sum rule, in which fluctuations belonging to different symmetries couple at finite momentum. It is demonstrated that the mode coupling between d-wave and s-wave pairing fluctuations leads to suppression of the d-wave fluctuation near the Mott insulator.

  4. GW and Bethe-Salpeter study of small water clusters

    NASA Astrophysics Data System (ADS)

    Blase, Xavier; Boulanger, Paul; Bruneval, Fabien; Fernandez-Serra, Marivi; Duchemin, Ivan

    2016-01-01

    We study within the GW and Bethe-Salpeter many-body perturbation theories the electronic and optical properties of small (H2O)n water clusters (n = 1-6). Comparison with high-level CCSD(T) Coupled-Cluster at the Single Double (Triple) levels and ADC(3) Green's function third order algebraic diagrammatic construction calculations indicates that the standard non-self-consistent G0W0@PBE or G0W0@PBE0 approaches significantly underestimate the ionization energy by about 1.1 eV and 0.5 eV, respectively. Consequently, the related Bethe-Salpeter lowest optical excitations are found to be located much too low in energy when building transitions from a non-self-consistent G0W0 description of the quasiparticle spectrum. Simple self-consistent schemes, with update of the eigenvalues only, are shown to provide a weak dependence on the Kohn-Sham starting point and a much better agreement with reference calculations. The present findings rationalize the theory to experiment possible discrepancies observed in previous G0W0 and Bethe-Salpeter studies of bulk water. The increase of the optical gap with increasing cluster size is consistent with the evolution from gas to dense ice or water phases and results from an enhanced screening of the electron-hole interaction.

  5. Self-Consistent Magnetosphere-Ionosphere Coupling

    NASA Technical Reports Server (NTRS)

    Six, N. Frank (Technical Monitor); Khazanov, G. V.; Newman, T. S.; Liemohn, M. W.; Fok, M. C.; Spiro, R. W.

    2002-01-01

    A self-consistent ring current (RC) model has been developed that couples electron and ion magnetospheric dynamics with the calculation of the electric field. Two new features were taken into account in order to close the self-consistent magnetosphere-ionosphere coupling loop. First, in addition to the RC ions, we have solved an electron kinetic equation in our model. Second, using the relation of Galand and Richmond, we have calculated the height integrated ionospheric conductances as a function of the precipitated high energy magnetospheric electrons and ions that are produced by our model. To validate the results of our model we simulate the magnetic storm of May 2, 1986, a storm that has been comprehensively studied by Fok et al., and have compared our results with different theoretical approaches. The self-consistent inclusion of the hot electrons and their effect on the conductance results in deeper penetration of the magnetospheric electric field. In addition, a slight westward rotation of the potential pattern (compared to previous self-consistent results) is evident in the inner magnetosphere. These effects change the hot plasma distribution, especially by allowing increased access of plasma sheet ions and electrons to low L shells.

  6. The ideas behind self-consistent expansion

    NASA Astrophysics Data System (ADS)

    Schwartz, Moshe; Katzav, Eytan

    2008-04-01

    In recent years we have witnessed a growing interest in various non-equilibrium systems described in terms of stochastic nonlinear field theories. In some of those systems, like KPZ and related models, the interesting behavior is in the strong coupling regime, which is inaccessible by traditional perturbative treatments such as dynamical renormalization group (DRG). A useful tool in the study of such systems is the self-consistent expansion (SCE), which might be said to generate its own 'small parameter'. The self-consistent expansion (SCE) has the advantage that its structure is just that of a regular expansion, the only difference is that the simple system around which the expansion is performed is adjustable. The purpose of this paper is to present the method in a simple and understandable way that hopefully will make it accessible to a wider public working on non-equilibrium statistical physics.

  7. Self-Consistent Random Phase Approximation

    NASA Astrophysics Data System (ADS)

    Rohr, Daniel; Hellgren, Maria; Gross, E. K. U.

    2012-02-01

    We report self-consistent Random Phase Approximation (RPA) calculations within the Density Functional Theory. The calculations are performed by the direct minimization scheme for the optimized effective potential method developed by Yang et al. [1]. We show results for the dissociation curve of H2^+, H2 and LiH with the RPA, where the exchange correlation kernel has been set to zero. For H2^+ and H2 we also show results for RPAX, where the exact exchange kernel has been included. The RPA, in general, over-correlates. At intermediate distances a maximum is obtained that lies above the exact energy. This is known from non-self-consistent calculations and is still present in the self-consistent results. The RPAX energies are higher than the RPA energies. At equilibrium distance they accurately reproduce the exact total energy. In the dissociation limit they improve upon RPA, but are still too low. For H2^+ the RPAX correlation energy is zero. Consequently, RPAX gives the exact dissociation curve. We also present the local potentials. They indicate that a peak at the bond midpoint builds up with increasing bond distance. This is expected for the exact KS potential.[4pt] [1] W. Yang, and Q. Wu, Phys. Rev. Lett., 89, 143002 (2002)

  8. Density neutron self-consistent caliper

    SciTech Connect

    Paske, W.C.; Rodney, P.F.; Roeder, R.A.

    1988-12-20

    This patent describes a system for determining the caliber of a borehole during drilling operations in an earth formation, comprising: first means adapted to make a first measurement of a first physical characteristic of an interior property of the formation; second means adapted to make a second measurement of a second physical characteristic of an interior property of the formation. The second physical characteristic being different from the first physical characteristic; means for determining the lithology of the formation; and means to compare the first and second measurements and to initiate an interation process based at least in part upon the determined lithology, to determine a self-consistent borehole caliber.

  9. Plasma Diffusion in Self-Consistent Fluctuations

    NASA Technical Reports Server (NTRS)

    Smets, R.; Belmont, G.; Aunai, N.

    2012-01-01

    The problem of particle diffusion in position space, as a consequence ofeleclromagnetic fluctuations is addressed. Numerical results obtained with a self-consistent hybrid code are presented, and a method to calculate diffusion coefficient in the direction perpendicular to the mean magnetic field is proposed. The diffusion is estimated for two different types of fluctuations. The first type (resuiting from an agyrotropic in itiai setting)is stationary, wide band white noise, and associated to Gaussian probability distribution function for the magnetic fluctuations. The second type (result ing from a Kelvin-Helmholtz instability) is non-stationary, with a power-law spectrum, and a non-Gaussian probabi lity distribution function. The results of the study allow revisiting the question of loading particles of solar wind origin in the Earth magnetosphere.

  10. Self-Consistent Scattering and Transport Calculations

    NASA Astrophysics Data System (ADS)

    Hansen, S. B.; Grabowski, P. E.

    2015-11-01

    An average-atom model with ion correlations provides a compact and complete description of atomic-scale physics in dense, finite-temperature plasmas. The self-consistent ionic and electronic distributions from the model enable calculation of x-ray scattering signals and conductivities for material across a wide range of temperatures and densities. We propose a definition for the bound electronic states that ensures smooth behavior of these measurable properties under pressure ionization and compare the predictions of this model with those of less consistent models for Be, C, Al, and Fe. SNL is a multi-program laboratory managed and operated by Sandia Corp., a wholly owned subsidiary of Lockheed Martin Corp, for the U.S. DoE NNSA under contract DE-AC04-94AL85000. This work was supported by DoE OFES Early Career grant FWP-14-017426.

  11. Plasma diffusion in self-consistent fluctuations

    SciTech Connect

    Smets, R.; Belmont, G.; Aunai, N.; Rezeau, L.

    2011-10-15

    The problem of particle diffusion in position space, as a consequence of electromagnetic fluctuations is addressed. Numerical results obtained with a self-consistent hybrid code are presented, and a method to calculate diffusion coefficient in the direction perpendicular to the mean magnetic field is proposed. The diffusion is estimated for two different types of fluctuations. The first type (resulting from an agyrotropic initial setting) is stationary, wide band white noise, and associated to Gaussian probability distribution function for the magnetic fluctuations. The second type (resulting from a Kelvin-Helmholtz instability) is non-stationary, with a power-law spectrum, and a non-Gaussian probability distribution function. The results of the study allow revisiting the question of loading particles of solar wind origin in the Earth magnetosphere.

  12. Quasiparticle properties of Ge(111)-2 times 1 surface

    SciTech Connect

    Zhu, X.; Louie, S.G.

    1992-08-01

    We have studied from first principles the quasiparticle properties of the 2 {times} 1 reconstructed (111) surface of Ge. Quasiparticle energies are calculated using the GW expansion of the electron self energy operator. The calculations explain a spectrum of experimental results obtained from photoemission, inverse photoemission, optical absorption, scanning tunneling microscopy, etc., for this surface. We also present a quasiparticle theory for the photoelectric threshold and examine the effects of many body corrections for this quantity.

  13. Self consistency grouping: a stringent clustering method

    PubMed Central

    2012-01-01

    Background Numerous types of clustering like single linkage and K-means have been widely studied and applied to a variety of scientific problems. However, the existing methods are not readily applicable for the problems that demand high stringency. Methods Our method, self consistency grouping, i.e. SCG, yields clusters whose members are closer in rank to each other than to any member outside the cluster. We do not define a distance metric; we use the best known distance metric and presume that it measures the correct distance. SCG does not impose any restriction on the size or the number of the clusters that it finds. The boundaries of clusters are determined by the inconsistencies in the ranks. In addition to the direct implementation that finds the complete structure of the (sub)clusters we implemented two faster versions. The fastest version is guaranteed to find only the clusters that are not subclusters of any other clusters and the other version yields the same output as the direct implementation but does so more efficiently. Results Our tests have demonstrated that SCG yields very few false positives. This was accomplished by introducing errors in the distance measurement. Clustering of protein domain representatives by structural similarity showed that SCG could recover homologous groups with high precision. Conclusions SCG has potential for finding biological relationships under stringent conditions. PMID:23320864

  14. Neutrinoless double {beta}-decay nuclear matrix elements within the SRQRPA with self-consistent short range correlations

    SciTech Connect

    Benes, Petr; Simkovic, Fedor

    2009-11-09

    The nuclear matrix elements M{sup 0v} of the neutrinoless double beta decay (0v{beta}{beta}-decay) are systematically evaluated using the self-consistent renormalized quasiparticle random phase approximation (SRQRPA). The residual interaction and the two-nucleon short-range correlations are derived from the charge-dependent Bonn (CD-Bonn) potential. The importance of further progress in the calculation of the 0v{beta}{beta}-decay nuclear matrix elements is stressed.

  15. Towards a unified description of ground and excited state properties: GW vs RPA and beyond

    NASA Astrophysics Data System (ADS)

    Rinke, Patrick

    2015-03-01

    In the quest for finding an ``optimal'' first principles electronic structure method, that combines accuracy and tractability with transferability across different chemical environments and dimensionalities (e.g. molecules, wires/tubes, surfaces, solids), the treatment of exchange and correlation in terms of ``exact-exchange plus correlation in the random-phase approximation (EX+cRPA)'' offers a promising avenue. Likewise one can express the same level of theory in the Green's function context through the GW approximation, which has the additional advantage that quasiparticle spectra as measured by direct and inverse photoemission become accessible. In this talk I will contrast both approaches and present the latest results from our continuous assessment. We find that self-consistent (sc) GW provides excellent charge densities, which is particularly important for charge transfer systems. Spectral properties for closed shell molecules are generally in good agreement with photoemission spectra, although a judicial choice of the starting point in perturbative G0W0 calculations can outperform scGW. Other ground state properties do not improve over EX+cRPA calculations. EX+cRPA, on the other hand, provides a good description of the ground state even for challenging cases like chemical reaction barrier heights and the f-electron metal cerium. The notorious underbinding of EX+cRPA can be corrected by going beyond RPA to renormalised second order perturbation theory (rPT2) that gives the overall most balanced performance. I will also discuss the associated rPT2 self-energy that goes beyond GW .

  16. Spin-Orbit Effects in the Quasiparticle Bandstructure of Noble Metals

    NASA Astrophysics Data System (ADS)

    Mustafa, Jamal; Louie, Steven

    2014-03-01

    Applications of the GW approximation to the electron self-energy have proven quite successful for calculating the quasiparticle properties of materials. We find that for the noble metals, in line with previous work in such calculations, the semicore states need to be taken into account. We show that, with these semicore states, a large cutoff must be used to describe the screening and, in turn, a large number of empty states must be included. Taking all of this into account, and carefully checking convergence, shows G0W0 can describe experimental results from angle-resolved photoemission spectroscopy quite well when the effects of spin-orbit coupling is also included. We compare our results to recent self-consistent GW calculations on gold. This work was supported by NSF grant No. DMR10-1006184 and U.S. DOE under Contract No. DE-AC02-05CH11231. Computational resources have been provided by DOE at LBNL's NERSC facility.

  17. Self-Consistent-Field Calculation on Lithium Hydride for Undergraduates.

    ERIC Educational Resources Information Center

    Rioux, Frank; Harriss, Donald K.

    1980-01-01

    Describes a self-consistent-field-linear combination of atomic orbitals-molecular orbital calculation on the valence electrons of lithium hydride using the method of Roothaan. This description is intended for undergraduate physics students.

  18. Self-consistent dynamics of a Josephson junction in the presence of an arbitrary environment.

    PubMed

    Joyez, Philippe

    2013-05-24

    We derive microscopically the dynamics associated with the dc Josephson effect in a superconducting tunnel junction interacting with an arbitrary electromagnetic environment. To do so, we extend to superconducting junctions the so-called P(E) theory (see, e.g., Ingold and Nazarov, arXiv:cond-mat/0508728) that accurately describes the interaction of a nonsuperconducting tunnel junction with its environment. We show the dynamics of this system is described by a small set of coupled correlation functions that take into account both Cooper pair and quasiparticle tunneling. When the phase fluctuations are small the problem is fully solved self-consistently, using and providing the exact linear admittance Y(ω) of the interacting junction. PMID:23745914

  19. GW100: Benchmarking G0W0 for Molecular Systems.

    PubMed

    van Setten, Michiel J; Caruso, Fabio; Sharifzadeh, Sahar; Ren, Xinguo; Scheffler, Matthias; Liu, Fang; Lischner, Johannes; Lin, Lin; Deslippe, Jack R; Louie, Steven G; Yang, Chao; Weigend, Florian; Neaton, Jeffrey B; Evers, Ferdinand; Rinke, Patrick

    2015-12-01

    We present the GW100 set. GW100 is a benchmark set of the ionization potentials and electron affinities of 100 molecules computed with the GW method using three independent GW codes and different GW methodologies. The quasi-particle energies of the highest-occupied molecular orbitals (HOMO) and lowest-unoccupied molecular orbitals (LUMO) are calculated for the GW100 set at the G0W0@PBE level using the software packages TURBOMOLE, FHI-aims, and BerkeleyGW. The use of these three codes allows for a quantitative comparison of the type of basis set (plane wave or local orbital) and handling of unoccupied states, the treatment of core and valence electrons (all electron or pseudopotentials), the treatment of the frequency dependence of the self-energy (full frequency or more approximate plasmon-pole models), and the algorithm for solving the quasi-particle equation. Primary results include reference values for future benchmarks, best practices for convergence within a particular approach, and average error bars for the most common approximations. PMID:26642984

  20. Quasiparticle energies and excitonic effects in dense solid hydrogen near metallization

    NASA Astrophysics Data System (ADS)

    Dvorak, Marc; Chen, Xiao-Jia; Wu, Zhigang

    2014-07-01

    We investigate the crucial metallization pressure of the Cmca-12 phase of solid hydrogen (H) using many-body perturbation theory within the GW approximation. We consider the effects of self-consistency, plasmon-pole models, and the vertex correction on the quasiparticle band gap (Eg). Our calculations show that self-consistency leads to an increase in Eg by 0.33 eV over the one-shot G0W0 approach. Because of error cancellation between the effects of self-consistency and the vertex correction, the simplest G0W0 method underestimates Eg by only 0.16 eV compared with the prediction of the more accurate GWΓ approach. Employing the plasmon-pole models underestimates Eg by 0.1-0.2 eV compared to the full-frequency numerical integration results. We thus predict a metallization pressure around 280 GPa, instead of 260 GPa predicted previously. Furthermore, we compute the optical absorption including the electron-hole interaction by solving the Bethe-Salpeter equation (BSE). The resulting absorption spectra demonstrate substantial redshifts and enhancement of absorption peaks compared to the calculated spectra neglecting excitonic effects. We find that the exciton binding energy decreases with increasing pressure from 66 meV at 100 GPa to 12 meV at 200 GPa due to the enhanced electronic screening as solid H approaches metallization. Because optical measurements are so important in identifying the structure of solid H, our BSE results should improve agreement between theory and experiment.

  1. Exact-exchange-based quasiparticle calculations

    SciTech Connect

    Aulbur, Wilfried G.; Staedele, Martin; Goerling, Andreas

    2000-09-15

    One-particle wave functions and energies from Kohn-Sham calculations with the exact local Kohn-Sham exchange and the local density approximation (LDA) correlation potential [EXX(c)] are used as input for quasiparticle calculations in the GW approximation (GWA) for eight semiconductors. Quasiparticle corrections to EXX(c) band gaps are small when EXX(c) band gaps are close to experiment. In the case of diamond, quasiparticle calculations are essential to remedy a 0.7 eV underestimate of the experimental band gap within EXX(c). The accuracy of EXX(c)-based GWA calculations for the determination of band gaps is as good as the accuracy of LDA-based GWA calculations. For the lowest valence band width a qualitatively different behavior is observed for medium- and wide-gap materials. The valence band width of medium- (wide-) gap materials is reduced (increased) in EXX(c) compared to the LDA. Quasiparticle corrections lead to a further reduction (increase). As a consequence, EXX(c)-based quasiparticle calculations give valence band widths that are generally 1-2 eV smaller (larger) than experiment for medium- (wide-) gap materials. (c) 2000 The American Physical Society.

  2. An approach to a self-consistent nuclear energy system

    SciTech Connect

    Fujii-e, Yoichi ); Arie, Kazuo; Endo, Hiroshi )

    1992-01-01

    A nuclear energy system should provide a stable supply of energy without endangering the environment or humans. If there is fear about exhausting world energy resources, accumulating radionuclides, and nuclear reactor safety, tension is created in human society. Nuclear energy systems of the future should be able to eliminate fear from people's minds. In other words, the whole system, including the nuclear fuel cycle, should be self-consistent. This is the ultimate goal of nuclear energy. If it can be realized, public acceptance of nuclear energy will increase significantly. In a self-consistent nuclear energy system, misunderstandings between experts on nuclear energy and the public should be minimized. The way to achieve this goal is to explain using simple logic. This paper proposes specific targets for self-consistent nuclear energy systems and shows that the fast breeder reactor (FBR) lies on the route to attaining the final goal.

  3. Analytic evaluation of the electronic self-energy in the GW approximation for two electrons on a sphere

    NASA Astrophysics Data System (ADS)

    Schindlmayr, Arno

    2013-02-01

    The GW approximation for the electronic self-energy is an important tool for the quantitative prediction of excited states in solids, but its mathematical exploration is hampered by the fact that it must, in general, be evaluated numerically even for very simple systems. In this paper I describe a nontrivial model consisting of two electrons on the surface of a sphere, interacting with the normal long-range Coulomb potential, and show that the GW self-energy, in the absence of self-consistency, can in fact be derived completely analytically in this case. The resulting expression is subsequently used to analyze the convergence of the energy gap between the highest occupied and the lowest unoccupied quasiparticle orbital with respect to the total number of states included in the spectral summations. The asymptotic formula for the truncation error obtained in this way, whose dominant contribution is proportional to the cutoff energy to the power -3/2, may be adapted to extrapolate energy gaps in other systems.

  4. Final Report Fermionic Symmetries and Self consistent Shell Model

    SciTech Connect

    Larry Zamick

    2008-11-07

    In this final report in the field of theoretical nuclear physics we note important accomplishments.We were confronted with "anomoulous" magnetic moments by the experimetalists and were able to expain them. We found unexpected partial dynamical symmetries--completely unknown before, and were able to a large extent to expain them.The importance of a self consistent shell model was emphasized.

  5. The Self-Consistency Model of Subjective Confidence

    ERIC Educational Resources Information Center

    Koriat, Asher

    2012-01-01

    How do people monitor the correctness of their answers? A self-consistency model is proposed for the process underlying confidence judgments and their accuracy. In answering a 2-alternative question, participants are assumed to retrieve a sample of representations of the question and base their confidence on the consistency with which the chosen…

  6. Efficient self-consistent quantum transport simulator for quantum devices

    SciTech Connect

    Gao, X. Mamaluy, D.; Nielsen, E.; Young, R. W.; Lilly, M. P.; Bishop, N. C.; Carroll, M. S.; Muller, R. P.; Shirkhorshidian, A.

    2014-04-07

    We present a self-consistent one-dimensional (1D) quantum transport simulator based on the Contact Block Reduction (CBR) method, aiming for very fast and robust transport simulation of 1D quantum devices. Applying the general CBR approach to 1D open systems results in a set of very simple equations that are derived and given in detail for the first time. The charge self-consistency of the coupled CBR-Poisson equations is achieved by using the predictor-corrector iteration scheme with the optional Anderson acceleration. In addition, we introduce a new way to convert an equilibrium electrostatic barrier potential calculated from an external simulator to an effective doping profile, which is then used by the CBR-Poisson code for transport simulation of the barrier under non-zero biases. The code has been applied to simulate the quantum transport in a double barrier structure and across a tunnel barrier in a silicon double quantum dot. Extremely fast self-consistent 1D simulations of the differential conductance across a tunnel barrier in the quantum dot show better qualitative agreement with experiment than non-self-consistent simulations.

  7. Efficient self-consistent quantum transport simulator for quantum devices

    NASA Astrophysics Data System (ADS)

    Gao, X.; Mamaluy, D.; Nielsen, E.; Young, R. W.; Shirkhorshidian, A.; Lilly, M. P.; Bishop, N. C.; Carroll, M. S.; Muller, R. P.

    2014-04-01

    We present a self-consistent one-dimensional (1D) quantum transport simulator based on the Contact Block Reduction (CBR) method, aiming for very fast and robust transport simulation of 1D quantum devices. Applying the general CBR approach to 1D open systems results in a set of very simple equations that are derived and given in detail for the first time. The charge self-consistency of the coupled CBR-Poisson equations is achieved by using the predictor-corrector iteration scheme with the optional Anderson acceleration. In addition, we introduce a new way to convert an equilibrium electrostatic barrier potential calculated from an external simulator to an effective doping profile, which is then used by the CBR-Poisson code for transport simulation of the barrier under non-zero biases. The code has been applied to simulate the quantum transport in a double barrier structure and across a tunnel barrier in a silicon double quantum dot. Extremely fast self-consistent 1D simulations of the differential conductance across a tunnel barrier in the quantum dot show better qualitative agreement with experiment than non-self-consistent simulations.

  8. SOCIAL COMPARISON, SELF-CONSISTENCY AND THE PRESENTATION OF SELF.

    ERIC Educational Resources Information Center

    MORSE, STANLEY J.; GERGEN, KENNETH J.

    TO DISCOVER HOW A PERSON'S (P) SELF-CONCEPT IS AFFECTED BY THE CHARACTERISTICS OF ANOTHER (O) WHO SUDDENLY APPEARS IN THE SAME SOCIAL ENVIRONMENT, SEVERAL QUESTIONNAIRES, INCLUDING THE GERGEN-MORSE (1967) SELF-CONSISTENCY SCALE AND HALF THE COOPERSMITH SELF-ESTEEM INVENTORY, WERE ADMINISTERED TO 78 UNDERGRADUATE MEN WHO HAD ANSWERED AN AD FOR WORK…

  9. Ishimori-I equation with self-consistent sources

    SciTech Connect

    Hu Juan; Hu Xingbiao; Tam, H.-W.

    2009-05-15

    In this paper, Grammian solutions of the Ishimori-I (Ish-I) equation are first obtained by Hirota's direct method. Utilizing the source generation procedure, this equation with self-consistent sources is then presented and the corresponding Grammian solutions are derived. Finally, as a simple case, the (1,1) dromion solution is examined.

  10. Self-Consistent Magnetosphere-Ionosphere Coupling: Theoretical Studies

    NASA Technical Reports Server (NTRS)

    Khazanov, G. V.; Newman, T. S.; Liemohn, M. W.; Fok, M.-C.; Spiro, R. W.; Six, N. Frank (Technical Monitor)

    2002-01-01

    A self-consistent ring current (RC) model has been developed that couples electron and ion magnetospheric dynamics with the calculation of the electric field. Two new features were taken into account in order to close the self-consistent magnetosphere- ionosphere coupling loop. First, in addition to the RC ions, we solve an electron kinetic equation in our model. Second, using the relation of Galand and Richmond [2001], we calculate the height integrated ionospheric conductances as a function of the precipitated high energy magnetospheric electrons and ions that are produced by our code. To validate the results of our model we simulate the magnetic storm of May 2, 1986, a storm that has has been comprehensively studied, and compare our results with different theoretical approaches. The self-consistent inclusion of the hot electrons and, their effect on the conductance results in deeper penetration of the magnetospheric electric field. In addition, a slight westward rotation of the potential pattern (compared to previous self-consistent results) is evident in the inner magnetosphere. This changes the hot plasma distribution, especially by allowing increased access of plasma sheet ions and electrons to low L shells.

  11. Evolution of inhomogeneous condensates: Self-consistent variational approach

    SciTech Connect

    Boyanovsky, D.; Cooper, F.; de Vega, H.J.; Sodano, P.

    1998-07-01

    We establish a self-consistent variational framework that allows us to study numerically the non-equilibrium evolution of non-perturbative inhomogeneous field configurations including quantum back reaction effects. After discussing the practical merits and disadvantages of different approaches we provide a closed set of local and renormalizable update equations that determine the dynamical evolution of inhomogeneous condensates and can be implemented numerically. These incorporate self-consistently the back reaction of quantum fluctuations and particle production. This program requires the solution of a self-consistent inhomogeneous problem to provide initial Cauchy data for the inhomogeneous condensates and Green`s functions. We provide a simple solvable ansatz for such an initial value problem for the sine-Gordon and {phi}{sup 4} quantum field theories in one spatial dimension. We compare exact known results of the sine-Gordon model to this simple ansatz. We also study the linear sigma model in the large N limit in three spatial dimensions as a microscopic model for pion production in ultrarelativistic collisions. We provide a solvable self-consistent ansatz for the initial value problem with cylindrical symmetry. For this case we also obtain a closed set of local and renormalized update equations that can be numerically implemented. A novel phenomenon of spinodal instabilities and pion production arises as a result of a Klein paradox for large amplitude inhomogeneous condensate configurations. {copyright} {ital 1998} {ital The American Physical Society}

  12. The development of an instrument to measure self-consistency.

    PubMed

    Zhan, L; Shen, C

    1994-09-01

    The maintenance of self-consistency is a task that engages elderly people, and it can be viewed as an indicator of how well a person can cope with stress in the ageing process. However, there is no reliable and valid instrument to date that measures this phenomenon. To help elderly people to accomplish the task warrants the necessity to develop an instrument to measure self-consistency. The purpose of this study was to develop an instrument to measure self-consistency in elderly people with chronic conditions. The Self-Consistency Scale (SCS) was designed and administered to hearing-impaired elderly people (n = 130) in the north-east part of the USA. Psychometric properties of the SCS were evaluated and the results indicated a promising psychometric integrity. The obtained alpha coefficient for the SCS total scale was 0.89, with a score range of 51-104 and a mean total score of 85.10 (SD = 11.04). Convergent validity of the SCS was established by correlating the SCS to a Visual Analogue Scale--A Sense of Self (VAS), r = 0.60 (P < 0.01). Divergent validity with the Geriatric Depression Scale (GDS) was established, r = -0.57 (P < 0.01). Maximum likelihood factor analysis with oblimin rotation resulted in a two-factor solution: Factor I, self-knowledge; and Factor II, stability of self-concept. PMID:7963057

  13. Self-consistent nonperturbative theory for classical systems.

    PubMed

    Mederos, L; Navascués, G; Velasco, E

    2002-01-01

    We construct a self-consistent nonperturbative theory for the structure and thermodynamics of a classical system of particles that goes beyond the usual approaches based on perturbation theory. Our theory, which gives accurate predictions for the phase diagram, is based on two ingredients: first, use is made of an exact expression for the free energy of a many-body system in terms of a reference system and a coupling integral connecting the latter to the final system; second, correlation functions may be very accurately approximated using a number of sum rules relating the radial distribution function with thermodynamic quantities. Consistency between the coupling integral expression and the sum rules may be achieved by means of a self-consistent process. PMID:11800760

  14. Self-consistent Castaing Distribution of Solar Wind Turbulent Fluctuations

    NASA Astrophysics Data System (ADS)

    Sorriso-Valvo, L.; Marino, R.; Lijoi, L.; Perri, S.; Carbone, V.

    2015-07-01

    The intermittent behavior of solar wind turbulent fluctuations has often been investigated through the modeling of their probability distribution functions (PDFs). Among others, the Castaing model has successfully been used in the past. In this paper, the energy dissipation field of solar wind turbulence has been studied for fast, slow, and polar wind samples recorded by Helios 2 and Ulysses spacecraft. The statistical description of the dissipation rate has then been used to remove intermittency through conditioning of the PDFs. Based on such observation, a self-consistent, parameter-free Castaing model is presented. The self-consistent model is tested against experimental PDFs, showing good agreement and supporting the picture of a multifractal energy cascade at the origin of solar wind intermittency.

  15. Self-consistent calculations of alpha-decay energies

    SciTech Connect

    Tolokonnikov, S. V.; Lutostansky, Yu. S.; Saperstein, E. E.

    2013-06-15

    On the basis of the self-consistent theory of finite Fermi systems, the energies of alphadecay chains were calculated for several new superheavy nuclei discovered recently in experiments of the Dubna-Livermore Collaboration headed by Yu.Ts. Oganessian. The approach in question is implemented on the basis of the generalized method of the density functional proposed by Fayans and his coauthors. The version used here relies on the functional DF3-a proposed recently for describing a wide array of nuclear data, including data on superheavy nuclei. A detailed comparison of the results obtained on this basis with the predictions of different approaches, including the self-consistent Skyrme-Hartree-Fock method, the micro-macro method in the version developed by Sobiczewski and his coauthors, and the phenomenological method of Liran and his coauthors, is performed. The resulting alpha-decay energies are used to calculate respective lifetimes with the aid of the phenomenological Parkhomenko-Sobiczewski formula.

  16. The Kelvin equation and self-consistent nucleation theory

    SciTech Connect

    Wilemski, G. |

    1995-07-15

    Issues of self-consistency are reviewed for several unary equilibrium size distributions based on the capillarity approximation. Some apparent difficulties of interpretation are resolved. In terms of the kinetic approach to nucleation theory, the influence of self-consistency on the nucleation rate is shown to arise entirely from differences in the dimer evaporation rates for nearly all versions of classical theory. The nucleation rate behavior of the Kelvin model is explored. In this model, the Kelvin equation is used to prescribe all cluster evaporation rates. Nucleation rates predicted by the Kelvin model are quantitatively similar to those of the self-consistent classical (SCC) theory, but not to other simple versions of the classical theory. This behavior arises entirely from the relatively close coincidence of the SCC and Kelvin dimer evaporation rates. This means that, for the distribution-based versions of classical theory, the SCC model is the closest analogue of the Kelvin model. Because the Kelvin equation is fundamentally inadequate for very small clusters, the close relationship between the Kelvin and SCC formulations indicates that both are equally lacking in fundamental justification. The Kelvin model may, however, have some pragmatic utility, and a simple analytical rate expression is also derived for it to simplify the calculation of nucleation rates for this model. {copyright} {ital 1995} {ital American} {ital Institute} {ital of} {ital Physics}.

  17. Superfluid 4He dynamics beyond quasiparticle excitations

    NASA Astrophysics Data System (ADS)

    Beauvois, K.; Campbell, C. E.; Dawidowski, J.; Fâk, B.; Godfrin, H.; Krotscheck, E.; Lauter, H.-J.; Lichtenegger, T.; Ollivier, J.; Sultan, A.

    2016-07-01

    The dynamics of superfluid 4He at and above the Landau quasiparticle regime is investigated by high-precision inelastic neutron scattering measurements of the dynamic structure factor. A highly structured response is observed above the familiar phonon-maxon-roton spectrum, characterized by sharp thresholds for phonon-phonon, maxon-roton, and roton-roton coupling processes. The experimental dynamic structure factor is compared to the calculation of the same physical quantity by a dynamic many-body theory including three-phonon processes self-consistently. The theory is found to provide a quantitative description of the dynamics of the correlated bosons for energies up to about three times that of the Landau quasiparticles.

  18. Quasiparticle band structure of HgSe

    SciTech Connect

    Rohlfing, M.; Louie, S.G.

    1998-04-01

    Motivated by a recent discussion about the existence of a fundamental gap in HgSe [Phys. Rev. Lett. {bold 78}, 3165 (1997)], we calculate the quasiparticle band structure of HgSe within the GW approximation for the electron self-energy. The band-structure results show that HgSe is a semimetal, which is in agreement with most experimental data. We observe a strong wave-vector dependence of the self-energy of the lowest conduction band, leading to an increased dispersion and a small effective mass. This may help to interpret recent photoemission spectroscopy measurements. {copyright} {ital 1998} {ital The American Physical Society}

  19. Stellar Turbulent Convection: A Self-consistent Model

    NASA Astrophysics Data System (ADS)

    Canuto, V. M.; Goldman, I.; Mazzitelli, I.

    1996-12-01

    We present a self-consistent model for stellar turbulent convection that is similar in spirit to the CM model (Canuto & Mazzitelli 1991) since it accounts for the full spectrum of the turbulent eddies rather than only one eddy, as done in the mixing length theory (MLT). The model differs from the CM model in the treatment of the rate of energy input nS(k) from the source that generates the turbulence. In the present model, nS(k) is controlled by both the source and the turbulence it ultimately generates, thus ensuring a self-consistent modeling of the turbulence. This improves the CM model in which nS(k) was taken to be equal to the growth rate of the linear unstable convective modes. However, since the formulation of a self-consistent treatment is far from simple, we were forced to use a representation of the nonlinear interactions less complete than the one in the CM model. The ensuing equations were solved numerically for a wide range of convective efficiencies. The results are the convective flux, the mean square turbulent velocity, the root mean squared turbulent pressure and the turbulent viscosity. We implemented the model in the ATON stellar structure code and computed the evolution of a solar model. The results are generally similar to those of the CM model and thus quite different from the MLT. The present model requires a smaller overshoot into the upper radiative zone than does the CM model, in accord with recent empirical estimates. Application to Population II stars and comparison with the very metal-poor globular cluster M68 yields an age in the range 11-12 Gyr. This is somewhat younger than the CM age, which in turn is younger than the corresponding MLT age, a result of possible cosmological interest

  20. a Self-Consistent Model of the Black Hole Evaporation

    NASA Astrophysics Data System (ADS)

    Kawai, Hikaru; Matsuo, Yoshinori; Yokokura, Yuki

    2013-06-01

    We construct a self-consistent model which describes a black hole from formation to evaporation including the backreaction from the Hawking radiation. In the case where a null shell collapses, at the beginning the evaporation occurs, but it stops eventually, and a horizon and singularity appear. On the other hand, in the generic collapse process of a continuously distributed null matter, the black hole evaporates completely without forming a macroscopically large horizon nor singularity. We also find a stationary solution in the heat bath, which can be regarded as a normal thermodynamic object.

  1. Self-Consistent Electromagnetic Modeling of Electron Sources

    SciTech Connect

    Hess, Mark

    2006-11-27

    The modeling of high-brightness electron sources, such as photoinjectors, requires a self-consistent technique for including the electromagnetics of tight electron bunches, as well as, the appropriate conductor boundary conditions of the source. A novel and effective technique for incorporating both of these effects utilizes time-dependent Green's functions. The advantages of this method are that Green's functions are generated by Delta function sources (making them ideal for tight bunch modeling) while simultaneously satisfying the conductor boundary condition. We demonstrate how these methods are used in a newly developed code called IRPSS (Indiana Rf Photocathode Source Simulator), and show initial simulations using IRPSS.

  2. Electrostatic Potentials from Self-Consistent Hirshfeld Atomic Charges.

    PubMed

    Van Damme, Sofie; Bultinck, Patrick; Fias, Stijn

    2009-02-10

    It is shown that molecular electrostatic potentials obtained from iterative or self-consistent Hirshfeld atomic point charges agree remarkably well with the ab initio computed electrostatic potentials. The iterative Hirshfeld scheme performs nearly as well as electrostatic potential derived atomic charges, having the advantage of allowing the definition of the atom in the molecule, rather than just yielding charges. The quality of the iterative Hirshfeld charges for computing electrostatic potentials is examined for a large set of molecules and compared to other commonly used techniques for population analysis. PMID:26610109

  3. Self-consistency based control scheme for magnetization dynamics

    SciTech Connect

    Albuquerque, G.; Miltat, J.; Thiaville, A.

    2001-06-01

    A numerical framework is presented for the solution of the Landau{endash}Lifshitz{endash}Gilbert equation of magnetization motion using a semi-implicit Crank{endash}Nicholson integration scheme. Along with the details of both space and time domain discretizations, we report on the development of a physically based self-consistency criterion that allows for a quantitative measurement of error in dynamic micromagnetic simulations. In essence, this criterion relies in recalculating from actual magnetization motion the imposed phenomenological damping constant. Test calculations were performed with special attention paid to the determination of suitable integration time steps. {copyright} 2001 American Institute of Physics.

  4. Applicability of self-consistent mean-field theory

    SciTech Connect

    Guo Lu; Sakata, Fumihiko; Zhao Enguang

    2005-02-01

    Within the constrained Hartree-Fock (CHF) theory, an analytic condition is derived to estimate whether a concept of the self-consistent mean field is realized in the level repulsive region. The derived condition states that an iterative calculation of the CHF equation does not converge when the quantum fluctuations coming from two-body residual interaction and quadrupole deformation become larger than a single-particle energy difference between two avoided crossing orbits. By means of numerical calculation, it is shown that the analytic condition works well for a realistic case.

  5. Wavelets in self-consistent electronic structure calculations

    SciTech Connect

    Wei, S.; Chou, M.Y.

    1996-04-01

    We report the first implementation of orthonormal wavelet bases in self-consistent electronic structure calculations within the local-density approximation. These local bases of different scales efficiently describe localized orbitals of interest. As an example, we studied two molecules, H{sub 2} and O{sub 2}, using pseudopotentials and supercells. Considerably fewer bases are needed compared with conventional plane-wave approaches, yet calculated binding properties are similar. Our implementation employs fast wavelet and Fourier transforms, avoiding evaluating any three-dimensional integral numerically. {copyright} {ital 1996 The American Physical Society.}

  6. Quasi-particle electronic band structure and alignment of the V-VI-VII semiconductors SbSI, SbSBr, and SbSeI for solar cells

    NASA Astrophysics Data System (ADS)

    Butler, Keith T.; McKechnie, Scott; Azarhoosh, Pooya; van Schilfgaarde, Mark; Scanlon, David O.; Walsh, Aron

    2016-03-01

    The ternary V-VI-VII chalcohalides consist of one cation and two anions. Trivalent antimony—with a distinctive 5s2 electronic configuration—can be combined with a chalcogen (e.g., S or Se) and halide (e.g., Br or I) to produce photoactive ferroelectric semiconductors with similarities to the Pb halide perovskites. We report—from relativistic quasi-particle self-consistent GW theory—that these materials have a multi-valley electronic structure with several electron and hole basins close to the band extrema. We predict ionisation potentials of 5.3-5.8 eV from first-principles for the three materials, and assess electrical contacts that will be suitable for achieving photovoltaic action from these unconventional compounds.

  7. Self-consistent chemical model of partially ionized plasmas

    SciTech Connect

    Arkhipov, Yu. V.; Baimbetov, F. B.; Davletov, A. E.

    2011-01-15

    A simple renormalization theory of plasma particle interactions is proposed. It primarily stems from generic properties of equilibrium distribution functions and allows one to obtain the so-called generalized Poisson-Boltzmann equation for an effective interaction potential of two chosen particles in the presence of a third one. The same equation is then strictly derived from the Bogolyubov-Born-Green-Kirkwood-Yvon (BBGKY) hierarchy for equilibrium distribution functions in the pair correlation approximation. This enables one to construct a self-consistent chemical model of partially ionized plasmas, correctly accounting for the close interrelation of charged and neutral components thereof. Minimization of the system free energy provides ionization equilibrium and, thus, permits one to study the plasma composition in a wide range of its parameters. Unlike standard chemical models, the proposed one allows one to study the system correlation functions and thereby to obtain an equation of state which agrees well with exact results of quantum-mechanical activity expansions. It is shown that the plasma and neutral components are strongly interrelated, which results in the short-range order formation in the corresponding subsystem. The mathematical form of the results obtained enables one to both firmly establish this fact and to determine a characteristic length of the structure formation. Since the cornerstone of the proposed self-consistent chemical model of partially ionized plasmas is an effective pairwise interaction potential, it immediately provides quite an efficient calculation scheme not only for thermodynamical functions but for transport coefficients as well.

  8. SELF-CONSISTENT SIZE AND VELOCITY DISTRIBUTIONS OF COLLISIONAL CASCADES

    SciTech Connect

    Pan, Margaret; Schlichting, Hilke E. E-mail: hilke@ucla.edu

    2012-03-10

    The standard theoretical treatment of collisional cascades derives a steady-state size distribution assuming a single constant velocity dispersion for all bodies regardless of size. Here we relax this assumption and solve self-consistently for the bodies' steady-state size and size-dependent velocity distributions. Specifically, we account for viscous stirring, dynamical friction, and collisional damping of the bodies' random velocities in addition to the mass conservation requirement typically applied to find the size distribution in a steady-state cascade. The resulting size distributions are significantly steeper than those derived without velocity evolution. For example, accounting self-consistently for the velocities can change the standard q = 3.5 power-law index of the Dohnanyi differential size spectrum to an index as large as q = 4. Similarly, for bodies held together by their own gravity, the corresponding power-law index range 2.88 < q < 3.14 of Pan and Sari can steepen to values as large as q = 3.26. Our velocity results allow quantitative predictions of the bodies' scale heights as a function of size. Together with our predictions, observations of the scale heights for different-sized bodies for the Kuiper belt, the asteroid belt, and extrasolar debris disks may constrain the mass and number of large bodies stirring the cascade as well as the colliding bodies' internal strengths.

  9. Self-consistent perturbation theory for two dimensional twisted bilayers

    NASA Astrophysics Data System (ADS)

    Shirodkar, Sharmila N.; Tritsaris, Georgios A.; Kaxiras, Efthimios

    Theoretical modeling and ab-initio simulations of two dimensional heterostructures with arbitrary angles of rotation between layers involve unrealistically large and expensive calculations. To overcome this shortcoming, we develop a methodology for weakly interacting heterostructures that treats the effect of one layer on the other as perturbation, and restricts the calculations to their primitive cells. Thus, avoiding computationally expensive supercells. We start by approximating the interaction potential between the twisted bilayers to that of a hypothetical configuration (viz. ideally stacked untwisted layers), which produces band structures in reasonable agreement with full-scale ab-initio calculations for commensurate and twisted bilayers of graphene (Gr) and Gr/hexagonal boron nitride (h-BN) heterostructures. We then self-consistently calculate the charge density and hence, interaction potential of the heterostructures. In this work, we test our model for bilayers of various combinations of Gr, h-BN and transition metal dichalcogenides, and discuss the advantages and shortcomings of the self-consistently calculated interaction potential. Department of Physics, Harvard University, Cambridge, Massachusetts 02138, USA.

  10. Self-Consistent Superthermal Electron Effects on Plasmaspheric Refilling

    NASA Technical Reports Server (NTRS)

    Liemohn, M. W.; Khazanov, G. V.; Moore, T. E.; Guiter, S. M.

    1997-01-01

    The effects of self-consistently including superthermal electrons in the definition of the ambipolar electric field are investigated for the case of plasmaspheric refilling after a geomagnetic storm. By using the total electron population in the hydrodynamic equations, a method for incorporating superthermal electron parameters in the electric field and electron temperature calculation is developed. Also, the ambipolar electric field is included in the kinetic equation for the superthermal electrons through a change of variables using the total energy and the first adiabatic invariant. Calculations based on these changes are performed by coupling time-dependent models of the thermal plasma and superthermal electrons. Results from this treatment of the electric field and the self-consistent development of the solution are discussed in detail. Specifically, there is a decreased thermal electron density in the plasmasphere during the first few minutes of refilling, a slightly accelerated proton shock front, and a decreased superthermal electron flux due to the deceleration by the electric field. The timescales of plasmaspheric refilling are discussed and determined to be somewhat shorter than previously calculated for the thermal plasma and superthermal electron population due to the effects of the field-aligned potential.

  11. Efficient self-consistency for magnetic tight binding

    NASA Astrophysics Data System (ADS)

    Soin, Preetma; Horsfield, A. P.; Nguyen-Manh, D.

    2011-06-01

    Tight binding can be extended to magnetic systems by including an exchange interaction on an atomic site that favours net spin polarisation. We have used a published model, extended to include long-ranged Coulomb interactions, to study defects in iron. We have found that achieving self-consistency using conventional techniques was either unstable or very slow. By formulating the problem of achieving charge and spin self-consistency as a search for stationary points of a Harris-Foulkes functional, extended to include spin, we have derived a much more efficient scheme based on a Newton-Raphson procedure. We demonstrate the capabilities of our method by looking at vacancies and self-interstitials in iron. Self-consistency can indeed be achieved in a more efficient and stable manner, but care needs to be taken to manage this. The algorithm is implemented in the code PLATO. Program summaryProgram title:PLATO Catalogue identifier: AEFC_v2_0 Program summary URL:http://cpc.cs.qub.ac.uk/summaries/AEFC_v2_0.html Program obtainable from: CPC Program Library, Queen's University, Belfast, N. Ireland Licensing provisions: Standard CPC licence, http://cpc.cs.qub.ac.uk/licence/licence.html No. of lines in distributed program, including test data, etc.: 228 747 No. of bytes in distributed program, including test data, etc.: 1 880 369 Distribution format: tar.gz Programming language: C and PERL Computer: Apple Macintosh, PC, Unix machines Operating system: Unix, Linux, Mac OS X, Windows XP Has the code been vectorised or parallelised?: Yes. Up to 256 processors tested RAM: Up to 2 Gbytes per processor Classification: 7.3 External routines: LAPACK, BLAS and optionally ScaLAPACK, BLACS, PBLAS, FFTW Catalogue identifier of previous version: AEFC_v1_0 Journal reference of previous version: Comput. Phys. Comm. 180 (2009) 2616 Does the new version supersede the previous version?: Yes Nature of problem: Achieving charge and spin self-consistency in magnetic tight binding can be very

  12. Isoscalar and isovector giant resonances in a self-consistent phonon coupling approach

    NASA Astrophysics Data System (ADS)

    Lyutorovich, N.; Tselyaev, V.; Speth, J.; Krewald, S.; Grümmer, F.; Reinhard, P.-G.

    2015-10-01

    We present fully self-consistent calculations of isoscalar giant monopole and quadrupole as well as isovector giant dipole resonances in heavy and light nuclei. The description is based on Skyrme energy-density functionals determining the static Hartree-Fock ground state and the excitation spectra within random-phase approximation (RPA) and RPA extended by including the quasiparticle-phonon coupling at the level of the time-blocking approximation (TBA). All matrix elements were derived consistently from the given energy-density functional and calculated without any approximation. As a new feature in these calculations, the single-particle continuum was included thus avoiding the artificial discretization usually implied in RPA and TBA. The step to include phonon coupling in TBA leads to small, but systematic, down shifts of the centroid energies of the giant resonances. These shifts are similar in size for all Skyrme parametrizations investigated here. After all, we demonstrate that one can find Skyrme parametrizations which deliver a good simultaneous reproduction of all three giant resonances within TBA.

  13. Electronic structure and metallization of cubic GdH3 under pressure: Ab initio many-body GW calculations

    NASA Astrophysics Data System (ADS)

    Kong, Bo; Zhang, Yachao

    2016-07-01

    The electronic structures of the cubic GdH3 are extensively investigated using the ab initio many-body GW calculations treating the Gd 4f electrons either in the core (4f-core) or in the valence states (4f-val). Different degrees of quasiparticle (QP) self-consistent calculations with the different starting points are used to correct the failures of the GGA/GGA + U/HSE03 calculations. In the 4f-core case, GGA + G0W0 calculations give a fundamental band gap of 1.72 eV, while GGA+ GW0 or GGA + GW calculations present a larger band gap. In the 4f-val case, the nonlocal exchange-correlation (xc) functional HSE03 can account much better for the strong localization of the 4f states than the semilocal or Hubbard U corrected xc functional in the Kohn-Sham equation. We show that the fundamental gap of the antiferromagnetic (AFM) or ferromagnetic (FM) GdH3 can be opened up by solving the QP equation with improved starting point of eigenvalues and wave functions given by HSE03. The HSE03 + G0W0 calculations present a fundamental band gap of 2.73 eV in the AFM configuration, and the results of the corresponding GW0 and GW calculations are 2.89 and 3.03 eV, respectively. In general, for the cubic structure, the fundamental gap from G0W0 calculations in the 4f-core case is the closest to the real result. By G0W0 calculations in the 4f-core case, we find that H or Gd defects can strongly affect the band structure, especially the H defects. We explain the mechanism in terms of the possible electron correlation on the hydrogen site. Under compression, the insulator-to-metal transition in the cubic GdH3 occurs around 40 GPa, which might be a satisfied prediction.

  14. Self-consistent quantum kinetic theory of diatomic molecule formation

    SciTech Connect

    Forrey, Robert C.

    2015-07-14

    A quantum kinetic theory of molecule formation is presented which includes three-body recombination and radiative association for a thermodynamically closed system which may or may not exchange energy with its surrounding at a constant temperature. The theory uses a Sturmian representation of a two-body continuum to achieve a steady-state solution of a governing master equation which is self-consistent in the sense that detailed balance between all bound and unbound states is rigorously enforced. The role of quasibound states in catalyzing the molecule formation is analyzed in complete detail. The theory is used to make three predictions which differ from conventional kinetic models. These predictions suggest significant modifications may be needed to phenomenological rate constants which are currently in wide use. Implications for models of low and high density systems are discussed.

  15. Linear Multigrid Techniques in Self-consistent Electronic Structure Calculations

    SciTech Connect

    Fattebert, J-L

    2000-05-23

    Ab initio DFT electronic structure calculations involve an iterative process to solve the Kohn-Sham equations for an Hamiltonian depending on the electronic density. We discretize these equations on a grid by finite differences. Trial eigenfunctions are improved at each step of the algorithm using multigrid techniques to efficiently reduce the error at all length scale, until self-consistency is achieved. In this paper we focus on an iterative eigensolver based on the idea of inexact inverse iteration, using multigrid as a preconditioner. We also discuss how this technique can be used for electrons described by general non-orthogonal wave functions, and how that leads to a linear scaling with the system size for the computational cost of the most expensive parts of the algorithm.

  16. Wakes in complex plasmas: A self-consistent kinetic theory

    NASA Astrophysics Data System (ADS)

    Kompaneets, Roman; Morfill, Gregor E.; Ivlev, Alexei V.

    2016-06-01

    In ground-based experiments with complex (dusty) plasmas, charged microparticles are levitated against gravity by an electric field, which also drives ion flow in the parent gas. Existing analytical approaches to describe the electrostatic interaction between microparticles in such conditions generally ignore the field and ion-neutral collisions, assuming free ion flow with a certain approximation for the ion velocity distribution function (usually a shifted Maxwellian). We provide a comprehensive analysis of our previously proposed self-consistent kinetic theory including the field, ion-neutral collisions, and the corresponding ion velocity distribution. We focus on various limiting cases and demonstrate how the interplay of these factors results in different forms of the shielding potential.

  17. Self-Consistent Ornstein-Zernike Approximation for Lattice Gases

    SciTech Connect

    Dickman, R.; Stell, G. |

    1996-08-01

    A self-consistent approximation for the structure factor of three-dimensional lattice gases yields remarkably accurate predictions (less than 3{percent} error over most of the temperature range) for the correlation length, isothermal compressibility, specific heat, and the coexistence curve. Critical temperatures agree to within 0.2{percent}, and other critical properties to within (1{endash}2){percent}, of the best numerical estimates. Until temperature and density are within 1{percent} of their critical values, the approximate {ital effective} critical exponents do not differ appreciably from their estimated exact form; they attain their limiting spherical-model values only much closer to critical. The method should prove useful for a variety of three-dimensional lattice-gas and fluid problems; it is inappropriate to two dimensions, where it predicts criticality at zero temperature. {copyright} {ital 1996 The American Physical Society.}

  18. Using Self Consistent Field Theory on Polymeric Mixtures

    NASA Astrophysics Data System (ADS)

    von Konigslow, Kier; Park, Chul; Thompson, Russell

    The ability to predict the solubility of a particular solvent in a polymer fluid is essential to the production of polymer foams. For the past 40 years, the primary model employed to this end has been an expansion of Flory-Huggins lattice fluid theory developed by Sanchez and Lacombe (S-L theory). S-L theory, while useful in the uniform limit, is limited to homogeneous systems. Self-Consistent Field Theory (SCFT), which has long been in use in polymer physics, is a mean-field theory capable of modeling the equilibrium behaviour of both homogeneous and inhomogeneous systems. We are investigating whether SCFT, applied to polymer-solvent mixtures, is in agreement with SL-theory in the homogeneous limit. Should this prove successful, we hope to use SCFT to model more general mixtures, including inhomogeneous nanocellular polymer foam systems.

  19. Soft core thermodynamics from self-consistent hard core fluids.

    PubMed

    Schöll-Paschinger, Elisabeth; Reiner, Albert

    2006-10-28

    In an effort to generalize the self-consistent Ornstein-Zernike approximation (SCOZA)-an accurate liquid state theory that has been restricted so far to hard core systems-to arbitrary soft core systems we study a combination of SCOZA with a recently developed perturbation theory. The latter was constructed by Ben-Amotz and Stell [J. Phys. Chem. B 108, 6877 (2004)] as a reformulation of the Weeks-Chandler-Andersen [J. Chem. Phys. 54, 5237 (1971)] perturbation theory directly in terms of an arbitrary hard sphere reference system. We investigate the accuracy of the combined approach for the Lennard-Jones fluid in comparison with simulation data and pure perturbation theory predictions and determine the dependence of the thermodynamic properties and the phase behavior on the choice of the effective hard core diameter of the reference system. PMID:17092101

  20. The Brittle-Ductile Transition - A Self-Consistent Approach.

    NASA Astrophysics Data System (ADS)

    Hobbs, B.; Regenauer-Lieb, K.; Ord, A.; Yuen, D. A.

    2006-12-01

    The brittle-ductile transition (BDT) in the Earth is commonly viewed as a switch between two different constitutive behaviors, plastic and viscous, and is represented in models by various formulations. We show that thermal-mechanical coupling leads to a self consistent view where the BDT emerges naturally within one constitutive framework once a critical temperature is attained. Viscous folding occurs above this temperature and brittle fracturing below. Seismic activity is maximised at the BDT. Orogenesis emerges as a thermal-mechanical decoupling near the BDT during flexing of the lithosphere with the development of "crocodile" -like structures, fold-nappe systems and far-travelled thrust sheets. For quartz- feldspar composite materials this transition lies in a critical range of 500 K to 580 K.

  1. Wakes in complex plasmas: A self-consistent kinetic theory.

    PubMed

    Kompaneets, Roman; Morfill, Gregor E; Ivlev, Alexei V

    2016-06-01

    In ground-based experiments with complex (dusty) plasmas, charged microparticles are levitated against gravity by an electric field, which also drives ion flow in the parent gas. Existing analytical approaches to describe the electrostatic interaction between microparticles in such conditions generally ignore the field and ion-neutral collisions, assuming free ion flow with a certain approximation for the ion velocity distribution function (usually a shifted Maxwellian). We provide a comprehensive analysis of our previously proposed self-consistent kinetic theory including the field, ion-neutral collisions, and the corresponding ion velocity distribution. We focus on various limiting cases and demonstrate how the interplay of these factors results in different forms of the shielding potential. PMID:27415371

  2. Self-consistent simulation of cyclotron autoresonance maser amplifiers

    SciTech Connect

    Pendergast, K.D.; Danly, B.G.; Temkin, R.J.; Wurtele, J.S.

    1988-04-01

    A self-consistent, one-dimensional model of the cyclotron autoresonance maser (CARM) amplifier is developed, and numerical simulations based on this model are described. Detailed studies of the CARM gain and efficiency for a wide range of initial energy and velocity spreads are presented. The interaction efficiency is found to be substantially increased when the axial magnetic field is tapered. For example, efficiencies of greater than 41 percent are obtained for a 140-GHz CARM amplifier with a tapered axial magnetic field and a 700-kV 4.5-A electron beam with parallel velocity spreads of less than 1 percent. A discussion of the nonlinear bandwidth and interaction sensitivity to axial field inhomogeneities is presented.

  3. Binary nucleation kinetics. I. Self-consistent size distribution

    SciTech Connect

    Wilemski, G.; Wyslouzil, B.E. ||

    1995-07-15

    Using the principle of detailed balance, we derive a new self-consistency requirement, termed the kinetic product rule, relating the evaporation coefficients and equilibrium cluster distribution for a binary system. We use this result to demonstrate and resolve an inconsistency for an idealized Kelvin model of nucleation in a simple binary mixture. We next examine several common forms for the equilibrium distribution of binary clusters based on the capillarity approximation and ideal vapor behavior. We point out fundamental deficiencies for each expression. We also show that each distribution yields evaporation coefficients that formally satisfy the new kinetic product rule but are physically unsatisfactory because they depend on the monomer vapor concentrations. We then propose a new form of the binary distribution function that is free of the deficiencies of the previous functions except for its reliance on the capillarity approximation. This new self-consistent classical (SCC) size distribution for binary clusters has the following properties: It satisfies the law of mass action; it reduces to an SCC unary distribution for clusters of a single component; and it produces physically acceptable evaporation rate coefficients that also satisfy the new kinetic product rule. Since it is possible to construct other examples of similarly well-behaved distributions, our result is not unique in this respect, but it does give reasonable predictions. As an illustration, we calculate binary nucleation rates and vapor activities for the ethanol--hexanol system at 260 K using the new SCC distribution and compare them to experimental results. The theoretical rates are uniformly higher than the experimental values over the entire vapor composition range. Although the predicted activities are lower, we find good agreement between the measured and theoretical slope of the critical vapor activity curve at a constant nucleation rate of 10{sup 7} cm{sup {minus}3} s{sup {minus}2}.

  4. Size-extensive vibrational self-consistent field method

    NASA Astrophysics Data System (ADS)

    Keçeli, Murat; Hirata, So

    2011-10-01

    The vibrational self-consistent field (VSCF) method is a mean-field approach to solve the vibrational Schrödinger equation and serves as a basis of vibrational perturbation and coupled-cluster methods. Together they account for anharmonic effects on vibrational transition frequencies and vibrationally averaged properties. This article reports the definition, programmable equations, and corresponding initial implementation of a diagrammatically size-extensive modification of VSCF, from which numerous terms with nonphysical size dependence in the original VSCF equations have been eliminated. When combined with a quartic force field (QFF), this compact and strictly size-extensive VSCF (XVSCF) method requires only quartic force constants of the partial ^4 V / partial Q_i^2 partial Q_j^2 type, where V is the electronic energy and Qi is the ith normal coordinate. Consequently, the cost of a XVSCF calculation with a QFF increases only quadratically with the number of modes, while that of a VSCF calculation grows quartically. The effective (mean-field) potential of XVSCF felt by each mode is shown to be harmonic, making the XVSCF equations subject to a self-consistent analytical solution without matrix diagonalization or a basis-set expansion, which are necessary in VSCF. Even when the same set of force constants is used, XVSCF is nearly three orders of magnitude faster than VSCF implemented similarly. Yet, the results of XVSCF and VSCF are shown to approach each other as the molecular size is increased, implicating the inclusion of unnecessary, nonphysical terms in VSCF. The diagrams of the XVSCF energy expression and their evaluation rules are also proposed, underscoring their connected structures.

  5. Modeling self-consistent multi-class dynamic traffic flow

    NASA Astrophysics Data System (ADS)

    Cho, Hsun-Jung; Lo, Shih-Ching

    2002-09-01

    In this study, we present a systematic self-consistent multiclass multilane traffic model derived from the vehicular Boltzmann equation and the traffic dispersion model. The multilane domain is considered as a two-dimensional space and the interaction among vehicles in the domain is described by a dispersion model. The reason we consider a multilane domain as a two-dimensional space is that the driving behavior of road users may not be restricted by lanes, especially motorcyclists. The dispersion model, which is a nonlinear Poisson equation, is derived from the car-following theory and the equilibrium assumption. Under the concept that all kinds of users share the finite section, the density is distributed on a road by the dispersion model. In addition, the dynamic evolution of the traffic flow is determined by the systematic gas-kinetic model derived from the Boltzmann equation. Multiplying Boltzmann equation by the zeroth, first- and second-order moment functions, integrating both side of the equation and using chain rules, we can derive continuity, motion and variance equation, respectively. However, the second-order moment function, which is the square of the individual velocity, is employed by previous researches does not have physical meaning in traffic flow. Although the second-order expansion results in the velocity variance equation, additional terms may be generated. The velocity variance equation we propose is derived from multiplying Boltzmann equation by the individual velocity variance. It modifies the previous model and presents a new gas-kinetic traffic flow model. By coupling the gas-kinetic model and the dispersion model, a self-consistent system is presented.

  6. Behavior of one-quasiparticle levels in odd isotonic chains of heavy nuclei

    SciTech Connect

    Adamian, G. G.; Antonenko, N. V.; Kuklin, S. N.; Malov, L. A.; Lu, B. N.; Zhou, S. G.

    2011-08-15

    The low-lying one-quasiparticle states are studied in the isotonic chains with N=147, 149, 151, 153, and 155 within the microscopic-macroscopic and self-consistent approaches. The energies of one-quasiparticle states change rather smoothly in the isotonic chains if there is no cross of the proton subshell. The {alpha}-decay schemes of several nuclei are suggested. The isomeric states in the odd isotopes of Fm and No are discussed.

  7. Quasiparticle energy studies of bulk semiconductors, surfaces and nanotubes

    SciTech Connect

    Blase, X.F.

    1994-12-01

    Effects of many-body effects on electronic excitation energies (quasiparticle band structure) of these materials are explored. GW approximation, including local field effects, for self-energy operator is used to calculate quasi-particle energies. The newly discovered carbon nanotubes are studied; structural stability and band structures are calculated. BN nanotubes are also studied, and their stability is predicted. Unexpected electronic features are predicted for both systems. Filling of carbon nanotubes with metal atoms and the doping of BN nanotubes by carbon and other impurites is also studied. The occupied surface states at H/Si(111)-(1{times}1) surface are studied; it is shown that the electronic structure requires a full quasiparticle calculation even for this simple chemisorption system. The core level shift of the Si 2p levels for atoms near the H/Si(111)-(1{times}1) surface is calculated; a simple first order perturbation theory using pseudopotential and the local density approximation gives good results for the photoemission spectra of the core electrons. The quasiparticle energies of bulk hexagonal BN and those of an isolated BN sheet are studied; this provides an understanding of the quasiparticle band structure of BN nanotubes. A nearly free electron state with a wavefunction in the interlayer or vacuum region composes the bottom of the conduction bands. A mixed-space formalism is presented for calculating the dynamical screening effects and electron self-energy operator in solids; this provides an efficient algorithm to calculate quasiparticle energies for large systems.

  8. Self-consistent dynamical and thermodynamical evolutions of protoplanetary disks.

    NASA Astrophysics Data System (ADS)

    Baillie, K.; Charnoz, S.; Taillifet, E.; Piau, L.

    2012-09-01

    Astronomical observations reveal the diversity of protoplanetary disk evolutions. In order to understand the global evolution of these disks from their birth, during the collapse of the molecular cloud, to their evaporation because of the stellar radiation, many processes with different timescales must be coupled: stellar evolution, thermodynamical evolution, photoevaporation, cloud collapse, viscous spreading... Simulating all these processes simultaneously is beyond the capacity of modern computers. However, by modeling the results of large scale simulations and coupling them with models of viscous evolution, we have designed a one dimension full model of disk evolution. In order to generate the most realistic protoplanetary disk, we minimize the number of input parameters and try to calculate most of them from self-consistent processes, as early as possible in the history of the disk; starting with the collapse of the molecular cloud that feeds the disk in gas. We start from the Hueso and Guillot, 2005 [2] model of disk evolution and couple the radiative transfer description of Calvet et al, 1991 [1] allowing us to handle a non-isothermal disk which midplane temperature is defined by an irradiation term form the central star and a viscous heating term depending on the optical depth of the disk. Our new model of the disk photosphere profile allows us to estimate self-consistent photosphere heights and midplane temperatures at the same time. We then follow the disk evolution using an upgrade of the viscous spreading equation from Lynden-Bell and Pringle, 1981 [3]. In particular, the molecular cloud collapse adds a time varying term to the temporal variation of the surface mass density of the disk, in the same manner that photo-evaporation introduces a density loss term. The central star itself is modeled using recent stellar evolution code described in Piau et al, 2011 [4]. Using the same temperature model in the vertical direction, we estimate 2D thermal maps of

  9. Self-consistent discharge growing model of helicon plasma

    NASA Astrophysics Data System (ADS)

    Isayama, Shogo; Hada, Tohru; Shinohara, Shunjiro; Tanikawa, Takao

    2015-11-01

    Helicon plasma is a high-density and low-temperature plasma generated by the electromagnetic (Helicon) wave excited in the plasma. It is thought to be useful for various applications including electric thrusters. Physics of helicon plasma production involves such fundamental processes as the wave propagation (dispersion relation), collisional and non-collisional wave damping, plasma heating, ionization/recombination of neutral particles, and modification of the dispersion relation by newly ionized plasma. There remain a number of unsolved physical issues such as, how the Helicon and the TG modes influence the plasma density, electron temperature and their spatial profiles. While the Helicon mode is absorbed in the bulk plasma, the TG mode is mostly absorbed near the edge of the plasma. The local power deposition in the helicon plasma is mostly balanced by collisional loss. This local power balance can give rise to the inhomogeneous electron temperature profile that leads to time evolution of density profile and dispersion relation. In our study, we construct a self-consistent model of the discharge evolution that includes the wave excitation, electron heat transfer, and diffusion of charged particles.

  10. Unifying Self-Consistent Field Theory for Weak Polyelectrolytes

    NASA Astrophysics Data System (ADS)

    Witte, Kevin; Won, You-Yeon

    2008-03-01

    A self-consistent field (SCF) theory for weak polyelectrolytes has been derived from a grand canonical partition function. The formalism accounts for the location and mixing of the charged and uncharged polymer species, treating the local (spatially dependent) charge fraction as a field variable with which to minimize the total free energy. This method of the derivation gives the resulting equations, especially those governing the local charge fraction, that are identical to the results obtained by Szleifer and coworkers (J. Polym. Sci. B Polym. Phys., 2006) who built upon the mean-field ``annealed'' free energy expression proposed by Raphael and Joanny (Europhys. Lett., 1990). However, we show that these results are further identical to the ``two-state'' model of Borukhov, Andelman and Orland (Eur. Phys. J. B, 1998), namely, the potential field due to the polymer charges with which the chains interact and the local charge fraction are shown to be exactly equal. This annealed model is derived by averaging the partition function with regard to the monomer charges. The charged and uncharged states are weighted by their probabilities which is, in our notation, the bulk charge fraction and one minus the bulk charge fraction, respectively. The utility of this theory is demonstrated by comparing its predictions against various experimental results from bulk potentiometric measurements and also from polyelectrolyte brush compression studies.

  11. Self-Consistent and Time-Dependent Solar Wind Models

    NASA Technical Reports Server (NTRS)

    Ong, K. K.; Musielak, Z. E.; Rosner, R.; Suess, S. T.; Sulkanen, M. E.

    1997-01-01

    We describe the first results from a self-consistent study of Alfven waves for the time-dependent, single-fluid magnetohydrodynamic (MHD) solar wind equations, using a modified version of the ZEUS MHD code. The wind models we examine are radially symmetrical and magnetized; the initial outflow is described by the standard Parker wind solution. Our study focuses on the effects of Alfven waves on the outflow and is based on solving the full set of the ideal nonlinear MHD equations. In contrast to previous studies, no assumptions regarding wave linearity, wave damping, and wave-flow interaction are made; thus, the models naturally account for the back-reaction of the wind on the waves, as well as for the nonlinear interaction between different types of MHD waves. Our results clearly demonstrate when momentum deposition by Alfven waves in the solar wind can be sufficient to explain the origin of fast streams in solar coronal holes; we discuss the range of wave amplitudes required to obtained such fast stream solutions.

  12. Self-consistent MPI-IO performance requirements and expectations.

    SciTech Connect

    Gropp, W. D.; Kimpe, D.; Ross, R.; Thakur, R.; Traff, J. L.; Mathematics and Computer Science; Univ. of Illinois; Katholieke Univ. Leuven; NEC Laboratories Europe

    2008-01-01

    We recently introduced the idea of self-consistent performance requirements for MPI communication. Such requirements provide a means to ensure consistent behavior of an MPI library, thereby ensuring a degree of performance portability by making it unnecessary for a user to perform implementation-dependent optimizations by hand. For the collective operations in particular, a large number of such rules could sensibly be formulated, without making hidden assumptions about the underlying communication system or otherwise constraining the MPI implementation. In this paper, we extend this idea to the realm of parallel I/O (MPI-IO), where the issues are far more subtle. In particular, it is not always possible to specify performance requirements without making assumptions about the implementation or without a priori knowledge of the I/O access pattern. For such cases, we introduce the notion of performance expectations, which specify the desired behavior for good implementations of MPI-IO. I/O performance requirements as well as expectations could be automatically checked by an appropriate benchmarking tool.

  13. Two regimes of self-consistent heating of charged particles

    SciTech Connect

    Doveil, Fabrice; Macor, Alessandro

    2011-10-15

    An experimental illustration of the transition between two basic transport regimes for wave-particle interaction is reported. A striking feature is that chaos, although present in both regimes, does not need to be invoked to explain the observed behavior. This experimental realization for a theoretical paradigm opens the possibility to check the validity of basic models, as is normally required in physics. Indeed, among seemingly simple problems exhibiting complex behavior is the classical interaction of an electron with electrostatic waves. Launching a low-intensity electron beam in a Traveling Wave Tube (TWT) recently allowed observing the real-world consequences of the richness of the electron trajectory. Here we show that self-consistent effects are nevertheless acting on a modulated test electron beam through the collective excitation of beam modes. The transition between two different particle transport regimes (stochastic diffusion in a set of waves and slow chaos associated to a pulsating separatrix) is directly measured by increasing the amplitude of the excitation.

  14. Two regimes of self-consistent heating of charged particles.

    PubMed

    Doveil, Fabrice; Macor, Alessandro

    2011-10-01

    An experimental illustration of the transition between two basic transport regimes for wave-particle interaction is reported. A striking feature is that chaos, although present in both regimes, does not need to be invoked to explain the observed behavior. This experimental realization for a theoretical paradigm opens the possibility to check the validity of basic models, as is normally required in physics. Indeed, among seemingly simple problems exhibiting complex behavior is the classical interaction of an electron with electrostatic waves. Launching a low-intensity electron beam in a Traveling Wave Tube (TWT) recently allowed observing the real-world consequences of the richness of the electron trajectory. Here we show that self-consistent effects are nevertheless acting on a modulated test electron beam through the collective excitation of beam modes. The transition between two different particle transport regimes (stochastic diffusion in a set of waves and slow chaos associated to a pulsating separatrix) is directly measured by increasing the amplitude of the excitation. PMID:22181220

  15. A new mixed self-consistent field procedure

    NASA Astrophysics Data System (ADS)

    Alvarez-Ibarra, A.; Köster, A. M.

    2015-10-01

    A new approach for the calculation of three-centre electronic repulsion integrals (ERIs) is developed, implemented and benchmarked in the framework of auxiliary density functional theory (ADFT). The so-called mixed self-consistent field (mixed SCF) divides the computationally costly ERIs in two sets: far-field and near-field. Far-field ERIs are calculated using the newly developed double asymptotic expansion as in the direct SCF scheme. Near-field ERIs are calculated only once prior to the SCF procedure and stored in memory, as in the conventional SCF scheme. Hence the name, mixed SCF. The implementation is particularly powerful when used in parallel architectures, since all RAM available are used for near-field ERI storage. In addition, the efficient distribution algorithm performs minimal intercommunication operations between processors, avoiding a potential bottleneck. One-, two- and three-dimensional systems are used for benchmarking, showing substantial time reduction in the ERI calculation for all of them. A Born-Oppenheimer molecular dynamics calculation for the Na+55 cluster is also shown in order to demonstrate the speed-up for small systems achievable with the mixed SCF. Dedicated to Sourav Pal on the occasion of his 60th birthday.

  16. Self-Consistent Monte Carlo Simulations of Positive Column Discharges

    NASA Astrophysics Data System (ADS)

    Lawler, J. E.; Kortshagen, U.

    1998-10-01

    In recent years it has become widely recognized that electron distribution functions in atomic gas positive column discharges are best described as non local over most of the range of R× N (column radius × gas density) where positive columns are stable. The use of an efficient Monte Carlo code with a radial potential expansion in powers of r^2 and with judiciously chosen constraints on the potential near the axis and wall now provides fully self-consistent kinetic solutions using only small computers. A set of solutions at smaller R× N and lower currents are presented which exhibit the classic negative dynamic resistance of the positive column at low currents. The negative dynamic resistance is due to a non-negligible Debye length and is sometimes described as a transition from free to ambipolar diffusion. This phenomenon is sensitive to radial variations of key parameters in the positive column and thus kinetic theory simulations are likely to provide a more realistic description than classic isothermal fluid models of the positive column. Comparisons of kinetic theory simulations to various fluid models of the positive column continue to provide new insight on this `corner stone' problem of Gaseous Electronics.

  17. First principles molecular dynamics without self-consistent field optimization

    SciTech Connect

    Souvatzis, Petros; Niklasson, Anders M. N.

    2014-01-28

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

  18. Self-consistent gyroviscous fluid model of rotational discontinuities

    SciTech Connect

    Hau, L.N.; Sonnerup, B.U.O. )

    1991-09-01

    One-dimensional steady state equilibrium structures of rotational discontinuities (RDs) have been constructed by use of a nondissipative gyroviscous two-fluid model in which electron inertia is neglected. The gyroviscous two-fluid RDs presented in this paper satisfy the MHD jump conditions for RDs exactly, i.e., all plasma properties as well as the magnetic field magnitude are identical and uniform on both sides of the layer, and v{sub x}, the normal flow speed there, is equal to the Alfven speed based on the normal magnetic field component B{sub x}. A necessary condition for the existence of RD solutions is derived by performing fixed-point analysis at the possible upstream and downstream states of RDs. This analysis shows that in the B{sub y}B{sub z} plane the magnetic structure near the fixed point is either a center or a saddle point. For those upstream conditions that give rise to a center point, no RD structure exists; when the fixed point examined is a saddle point, numerical integration of the one-dimension, steady state, nonlinear gyroviscous two-fluid equations indicates that self-consistent complete RD structures may sometimes (but not always) exist. These structures involve a rotation of the tangential magnetic field by an angle {phi} that is determined by the angle {theta} between the upstream (or downstream) magnetic field and the vector normal to the discontinuity, by the plasma beta value, {beta} = 2{mu}{sub o}p/B{sup 2}, and by the electron to ion temperature ratio T{sub e}/T{sub i}. Only electron-polarized structures have been found to date: within these discontinuities the field magnitude {vert bar}B{vert bar} always displays a maximum. If ion-polarized RDs exist, they would contain a minimum in {vert bar}B{vert bar}. However, no such structures have been found.

  19. Self-consistent theory for systems with mesoscopic fluctuations.

    PubMed

    Ciach, A; Góźdź, W T

    2016-10-19

    We have developed a theory for inhomogeneous systems that allows for the incorporation of the effects of mesoscopic fluctuations. A hierarchy of equations relating the correlation and direct correlation functions for the local excess [Formula: see text] of the volume fraction of particles ζ has been obtained, and an approximation leading to a closed set of equations for the two-point functions has been introduced for the disordered inhomogeneous phase. We have numerically solved the self-consistent equations for one-dimensional (1D) and three-dimensional (3D) models with short-range attraction and long-range repulsion. Predictions for all of the qualitative properties of the 1D model agree with the exact results, but only semi-quantitative agreement is obtained in the simplest version of the theory. The effects of fluctuations in the two 3D models considered are significantly different, despite the very similar properties of these models in the mean-field approximation. In both cases we obtain the sequence of large-small-large compressibility for increasing ζ. The very small compressibility is accompanied by the oscillatory decay of correlations with correlation lengths that are orders of magnitude larger than the size of particles. In one of the two models considered, the small compressibility becomes very small and the large compressibility becomes very large with decreasing temperature, and eventually van der Waals loops appear. Further studies are necessary in order to determine the nature of the strongly inhomogeneous phase present for intermediate volume fractions in 3D. PMID:27545343

  20. Self-consistent Dynamical And Thermodynamical Evolutions Of Protoplanetary Disks.

    NASA Astrophysics Data System (ADS)

    Baillie, Kevin; Charnoz, S.; Taillifet, E.; Piau, L.

    2012-10-01

    Astronomical observations reveal the diversity of protoplanetary disk evolutions. In order to understand the global evolution of these disks from their birth, during the collapse of the molecular cloud, to their evaporation because of the stellar radiation, many processes with different timescales must be coupled: stellar evolution, thermodynamical evolution, photoevaporation, cloud collapse, viscous spreading...By modeling the results of large scale simulations and coupling them with models of viscous evolution, we have designed a one dimension full model of disk evolution for which most parameters are self-consistently calculated at each time step. We integrate the Hueso and Guillot, 2005 model of disk evolution and couple the radiative transfer description of Calvet et al, 1991 allowing us to handle a non-isothermal disk. We also take into account the collapse of the molecular cloud that feeds the disk. Using the same temperature model in the vertical direction, we estimate 2D thermal maps of the disk. The central star itself is modeled using recent stellar evolution code described in Piau et al, 2011. We first test our model in the case of an already formed Minimum Mass solar Nebula, trying to match the observational constraints on the radial surface density gradients and photosphere height profiles of the Taurus-Auriga or Ophiucus disks for instance. We then follow the full long-term evolution of a disk fed by the collapse of the molecular cloud. We estimate disk temperatures and accretion rates and try to constrain the favourable zone for the formation of the first solids. This will help targeting future JWST observations.

  1. Self-consistent dynamical and thermodynamical evolutions of protoplanetary disks

    NASA Astrophysics Data System (ADS)

    Baillié, K.; Charnoz, S.; Taillifet, E.; Piau, L.

    2012-12-01

    Astronomical observations reveal the diversity of protoplanetary disk evolutions. In order to understand the global evolution of these disks from their birth, during the collapse of the molecular cloud, to their evaporation because of the stellar radiation, many processes with different timescales must be coupled: stellar evolution, thermodynamical evolution, photoevaporation, cloud collapse, viscous spreading... By modeling the results of large scale simulations and coupling them with models of viscous evolution, we have designed a one dimension full model of disk evolution for which most parameters are self-consistently calculated at each time step. We integrate the Hueso and Guillot, 2005 model of disk evolution and couple the radiative transfer description of Calvet et al, 1991 allowing us to handle a non-isothermal disk. We also take into account the collapse of the molecular cloud that feeds the disk. Using the same temperature model in the vertical direction, we estimate 2D thermal maps of the disk. The central star itself is modeled using recent stellar evolution code described in Piau et al, 2011. We first test our model in the case of an already formed Minimum Mass solar Nebula, trying to match the observational constraints on the radial surface density gradients and photosphere height profiles of the Taurus-Auriga or Ophiucus disks for instance. We then follow the full long-term evolution of a disk fed by the collapse of the molecular cloud. We estimate disk temperatures and accretion rates and try to constrain the favourable zone for the formation of the first solids. This will help targeting future JWST observations.

  2. Accurate Ionization Potentials and Electron Affinities of Acceptor Molecules III: A Benchmark of GW Methods.

    PubMed

    Knight, Joseph W; Wang, Xiaopeng; Gallandi, Lukas; Dolgounitcheva, Olga; Ren, Xinguo; Ortiz, J Vincent; Rinke, Patrick; Körzdörfer, Thomas; Marom, Noa

    2016-02-01

    The performance of different GW methods is assessed for a set of 24 organic acceptors. Errors are evaluated with respect to coupled cluster singles, doubles, and perturbative triples [CCSD(T)] reference data for the vertical ionization potentials (IPs) and electron affinities (EAs), extrapolated to the complete basis set limit. Additional comparisons are made to experimental data, where available. We consider fully self-consistent GW (scGW), partial self-consistency in the Green's function (scGW0), non-self-consistent G0W0 based on several mean-field starting points, and a "beyond GW" second-order screened exchange (SOSEX) correction to G0W0. We also describe the implementation of the self-consistent Coulomb hole with screened exchange method (COHSEX), which serves as one of the mean-field starting points. The best performers overall are G0W0+SOSEX and G0W0 based on an IP-tuned long-range corrected hybrid functional with the former being more accurate for EAs and the latter for IPs. Both provide a balanced treatment of localized vs delocalized states and valence spectra in good agreement with photoemission spectroscopy (PES) experiments. PMID:26731609

  3. Self-consistent Modeling of Elastic Anisotropy in Shale

    NASA Astrophysics Data System (ADS)

    Kanitpanyacharoen, W.; Wenk, H.; Matthies, S.; Vasin, R.

    2012-12-01

    Elastic anisotropy in clay-rich sedimentary rocks has increasingly received attention because of significance for prospecting of petroleum deposits, as well as seals in the context of nuclear waste and CO2 sequestration. The orientation of component minerals and pores/fractures is a critical factor that influences elastic anisotropy. In this study, we investigate lattice and shape preferred orientation (LPO and SPO) of three shales from the North Sea in UK, the Qusaiba Formation in Saudi Arabia, and the Officer Basin in Australia (referred to as N1, Qu3, and L1905, respectively) to calculate elastic properties and compare them with experimental results. Synchrotron hard X-ray diffraction and microtomography experiments were performed to quantify LPO, weight proportions, and three-dimensional SPO of constituent minerals and pores. Our preliminary results show that the degree of LPO and total amount of clays are highest in Qu3 (3.3-6.5 m.r.d and 74vol%), moderately high in N1 (2.4-5.6 m.r.d. and 70vol%), and lowest in L1905 (2.3-2.5 m.r.d. and 42vol%). In addition, porosity in Qu3 is as low as 2% while it is up to 6% in L1605 and 8% in N1, respectively. Based on this information and single crystal elastic properties of mineral components, we apply a self-consistent averaging method to calculate macroscopic elastic properties and corresponding seismic velocities for different shales. The elastic model is then compared with measured acoustic velocities on the same samples. The P-wave velocities measured from Qu3 (4.1-5.3 km/s, 26.3%Ani.) are faster than those obtained from L1905 (3.9-4.7 km/s, 18.6%Ani.) and N1 (3.6-4.3 km/s, 17.7%Ani.). By making adjustments for pore structure (aspect ratio) and single crystal elastic properties of clay minerals, a good agreement between our calculation and the ultrasonic measurement is obtained.

  4. Fast RPA and GW calculations: cubic system size scaling

    NASA Astrophysics Data System (ADS)

    Kresse, Georg

    The random phase approximation (RPA) to the correlation energy and the related GW approximation are among the most promising methods to obtain accurate correlation energy differences and QP energies from diagrammatic perturbation theory at reasonable computational cost. The calculations are, however, usually one to two orders of magnitude more demanding than conventional density functional theory calculations. Here, we show that a cubic system size scaling can be readily obtained reducing the computation time by one to two orders of magnitude for large systems. Furthermore, the scaling with respect to the number of k points used to sample the Brillouin zone can be reduced to linear order. In combination, this allows accurate and very well-converged single-point RPA and GW calculations, with a time complexity that is roughly on par or better than for self-consistent Hartree-Fock and hybrid-functional calculations. Furthermore, the talk discusses the relation between the RPA correlation energy and the GW approximation: the self-energy is the derivative of the RPA correlation energy with respect to the Green's function. The calculated self-energy can be used to compute QP-energies in the GW approximation, any first derivative of the total energy, as well as corrections to the correlation energy from the changes of the charge density when switching from DFT to a many-body body description (GW singles energy contribution).

  5. Relativistic quasiparticle time blocking approximation: Dipole response of open-shell nuclei

    NASA Astrophysics Data System (ADS)

    Litvinova, E.; Ring, P.; Tselyaev, V.

    2008-07-01

    The self-consistent relativistic quasiparticle random-phase approximation (RQRPA) is extended by the quasiparticle-phonon coupling (QPC) model using the quasiparticle time blocking approximation (QTBA). The method is formulated in terms of the Bethe-Salpeter equation (BSE) in the two-quasiparticle space with an energy-dependent two-quasiparticle residual interaction. This equation is solved either in the basis of Dirac states forming the self-consistent solution of the ground state or in the momentum representation. Pairing correlations are treated within the Bardeen-Cooper-Schrieffer (BCS) model with a monopole-monopole interaction. The same NL3 set of the coupling constants generates the Dirac-Hartree-BCS single-quasiparticle spectrum, the static part of the residual two-quasiparticle interaction and the quasiparticle-phonon coupling amplitudes. A quantitative description of electric dipole excitations in the chain of tin isotopes (Z=50) with the mass numbers A=100,106,114,116,120, and 130 and in the chain of isotones with (N=50) Sr88, Zr90, Mo92 is performed within this framework. The RQRPA extended by the coupling to collective vibrations generates spectra with a multitude of 2q⊗phonon (two quasiparticles plus phonon) states providing a noticeable fragmentation of the giant dipole resonance as well as of the soft dipole mode (pygmy resonance) in the nuclei under investigation. The results obtained for the photo absorption cross sections and for the integrated contributions of the low-lying strength to the calculated dipole spectra agree very well with the available experimental data.

  6. Self-consistent Pauli corrections in Brueckner-Hartree-Fock calculations

    NASA Technical Reports Server (NTRS)

    Braley, R. C.; Ford, W. F.

    1972-01-01

    A scheme is introduced which makes it feasible to make completely self-consistent Brueckner-Hartree-Fock (BHF) and renormalized BHF calculations for spherical, closed-shell and axially-symmetric deformed nuclei. The usual requirement or orbital self-consistency has been imposed, as well as self-consistency in the starting energies and occupation probabilities. Previously, only approximate forms were used for the Pauli operator. This approximation is removed and a method for making the necessary Pauli corrections to the reaction matrix during the approach to self-consistency is presented. A discussion of the symmetries which reduce the problem to one of manageable proportions is included.

  7. Quasiparticle theory of transport coefficients for hadronic matter at finite temperature and baryon density

    NASA Astrophysics Data System (ADS)

    Albright, M.; Kapusta, J. I.

    2016-01-01

    We develop a flexible quasiparticle theory of transport coefficients of hot hadronic matter at finite baryon density. We begin with a hadronic quasiparticle model which includes a scalar and a vector mean field. Quasiparticle energies and the mean fields depend on temperature and baryon chemical potential. Starting with the quasiparticle dispersion relation, we derive the Boltzmann equation and use the Chapman-Enskog expansion to derive formulas for the shear and bulk viscosities and thermal conductivity. We obtain both relaxation-time approximation formulas and more general integral equations. Throughout the work, we explicitly enforce the Landau-Lifshitz conditions of fit and ensure the theory is thermodynamically self-consistent. The derived formulas should be useful for predicting the transport coefficients of the hadronic phase of matter produced in heavy-ion collisions at the Relativistic Heavy Ion Collider and at other accelerators.

  8. Excitation spectra of aromatic molecules within a real-space G W -BSE formalism: Role of self-consistency and vertex corrections

    NASA Astrophysics Data System (ADS)

    Hung, Linda; da Jornada, Felipe H.; Souto-Casares, Jaime; Chelikowsky, James R.; Louie, Steven G.; Ã-ǧüt, Serdar

    2016-08-01

    We present first-principles calculations on the vertical ionization potentials (IPs), electron affinities (EAs), and singlet excitation energies on an aromatic-molecule test set (benzene, thiophene, 1,2,5-thiadiazole, naphthalene, benzothiazole, and tetrathiafulvalene) within the G W and Bethe-Salpeter equation (BSE) formalisms. Our computational framework, which employs a real-space basis for ground-state and a transition-space basis for excited-state calculations, is well suited for high-accuracy calculations on molecules, as we show by comparing against G0W0 calculations within a plane-wave-basis formalism. We then generalize our framework to test variants of the G W approximation that include a local density approximation (LDA)-derived vertex function (ΓLDA) and quasiparticle-self-consistent (QS) iterations. We find that ΓLDA and quasiparticle self-consistency shift IPs and EAs by roughly the same magnitude, but with opposite sign for IPs and the same sign for EAs. G0W0 and QS G W ΓLDA are more accurate for IPs, while G0W0ΓLDA and QS G W are best for EAs. For optical excitations, we find that perturbative G W -BSE underestimates the singlet excitation energy, while self-consistent G W -BSE results in good agreement with previous best-estimate values for both valence and Rydberg excitations. Finally, our work suggests that a hybrid approach, in which G0W0 energies are used for occupied orbitals and G0W0ΓLDA for unoccupied orbitals, also yields optical excitation energies in good agreement with experiment but at a smaller computational cost.

  9. Evaluating the GW Approximation with CCSD(T) for Charged Excitations Across the Oligoacenes.

    PubMed

    Rangel, Tonatiuh; Hamed, Samia M; Bruneval, Fabien; Neaton, Jeffrey B

    2016-06-14

    Charged excitations of the oligoacene family of molecules, relevant for astrophysics and technological applications, are widely studied and therefore provide an excellent system for benchmarking theoretical methods. In this work, we evaluate the performance of many-body perturbation theory within the GW approximation relative to new high-quality CCSD(T) reference data for charged excitations of the acenes. We compare GW calculations with a number of hybrid density functional theory starting points and with eigenvalue self-consistency. Special focus is given to elucidating the trend of GW-predicted excitations with molecule length increasing from benzene to hexacene. We find that GW calculations with starting points based on an optimally tuned range-separated hybrid (OTRSH) density functional and eigenvalue self-consistency can yield quantitative ionization potentials for the acenes. However, for larger acenes, the predicted electron affinities can deviate considerably from reference values. Our work paves the way for predictive and cost-effective GW calculations of charged excitations of molecules and identifies certain limitations of current GW methods used in practice for larger molecules. PMID:27123935

  10. Finite amplitude method for the quasiparticle random-phase approximation

    SciTech Connect

    Avogadro, Paolo; Nakatsukasa, Takashi

    2011-07-15

    We present the finite amplitude method (FAM), originally proposed in Ref. [17], for superfluid systems. A Hartree-Fock-Bogoliubov code may be transformed into a code of the quasiparticle-random-phase approximation (QRPA) with simple modifications. This technique has advantages over the conventional QRPA calculations, such as coding feasibility and computational cost. We perform the fully self-consistent linear-response calculation for the spherical neutron-rich nucleus {sup 174}Sn, modifying the hfbrad code, to demonstrate the accuracy, feasibility, and usefulness of the FAM.

  11. 0{nu}{beta}{beta}-decay nuclear matrix elements with self-consistent short-range correlations

    SciTech Connect

    Simkovic, Fedor; Faessler, Amand; Muether, Herbert; Rodin, Vadim; Stauf, Markus

    2009-05-15

    A self-consistent calculation of nuclear matrix elements of the neutrinoless double-beta decays (0{nu}{beta}{beta}) of {sup 76}Ge, {sup 82}Se, {sup 96}Zr, {sup 100}Mo, {sup 116}Cd, {sup 128}Te, {sup 130}Te, and {sup 136}Xe is presented in the framework of the renormalized quasiparticle random phase approximation (RQRPA) and the standard QRPA. The pairing and residual interactions as well as the two-nucleon short-range correlations are for the first time derived from the same modern realistic nucleon-nucleon potentials, namely, from the charge-dependent Bonn potential (CD-Bonn) and the Argonne V18 potential. In a comparison with the traditional approach of using the Miller-Spencer Jastrow correlations, matrix elements for the 0{nu}{beta}{beta} decay are obtained that are larger in magnitude. We analyze the differences among various two-nucleon correlations including those of the unitary correlation operator method (UCOM) and quantify the uncertainties in the calculated 0{nu}{beta}{beta}-decay matrix elements.

  12. Quasiparticle and optical properties of polythiophene-derived polymers

    NASA Astrophysics Data System (ADS)

    Samsonidze, Georgy; Ribeiro, Filipe J.; Cohen, Marvin L.; Louie, Steven G.

    2014-07-01

    Electron donor conjugated polymers blended with electron acceptor fullerene derivatives is one of the promising technologies for organic photovoltaics. However, with the energy conversion efficiency of only 9% in a single bulk heterojunction device structure, these solar cells are not yet competitive with conventional inorganic semiconductor technology. Some of the limitations are large optical band gaps and small electron affinities of polymers preventing the absorption of infrared radiation and leading to energy losses during charge separation at the donor-acceptor interface, respectively. In this work, we compute from first principles the quasiparticle and optical spectra of several different thiophene-, ethyne-, and vinylene-based copolymers using the GW method and the GW plus Bethe-Salpeter equation approach. One of the polymers is identified which has a preferential alignment of the energy levels at the interface with fullerene molecule compared to the reference case of polythiophene.

  13. Self-consistent Green's function embedding for advanced electronic structure methods based on a dynamical mean-field concept

    NASA Astrophysics Data System (ADS)

    Chibani, Wael; Ren, Xinguo; Scheffler, Matthias; Rinke, Patrick

    2016-04-01

    We present an embedding scheme for periodic systems that facilitates the treatment of the physically important part (here a unit cell or a supercell) with advanced electronic structure methods, that are computationally too expensive for periodic systems. The rest of the periodic system is treated with computationally less demanding approaches, e.g., Kohn-Sham density-functional theory, in a self-consistent manner. Our scheme is based on the concept of dynamical mean-field theory formulated in terms of Green's functions. Our real-space dynamical mean-field embedding scheme features two nested Dyson equations, one for the embedded cluster and another for the periodic surrounding. The total energy is computed from the resulting Green's functions. The performance of our scheme is demonstrated by treating the embedded region with hybrid functionals and many-body perturbation theory in the GW approach for simple bulk systems. The total energy and the density of states converge rapidly with respect to the computational parameters and approach their bulk limit with increasing cluster (i.e., computational supercell) size.

  14. On the periodic Toda lattice with a self-consistent source

    NASA Astrophysics Data System (ADS)

    Babajanov, Bazar; Fečkan, Michal; Urazboev, Gayrat

    2015-05-01

    This work is devoted to the application of inverse spectral problem for integration of the periodic Toda lattice with self-consistent source. The effective method of solution of the inverse spectral problem for the discrete Hill's equation is presented.

  15. Self-consistent Debye-Waller factors of the electron solid on liquid helium

    NASA Astrophysics Data System (ADS)

    Namaizawa, H.

    1980-05-01

    Based on the self-consistent field formalism we propose a shear-mode self-consistency for the high-frequency Debye-Waller factors (HFDWF) of the electron solid bound on a free surface of liquid helium. Our results are qualitatively in agreement with the empirical DW factor determined by Fisher, Halperin, and Platzman with the experiment of Grimes and Adams. We also report the analysis of the HFDWF according to the Lindemann law.

  16. Microscopic Model of Quasiparticle Wave Packets in Superfluids, Superconductors, and Paired Hall States

    NASA Astrophysics Data System (ADS)

    Parameswaran, S. A.; Kivelson, S. A.; Shankar, R.; Sondhi, S. L.; Spivak, B. Z.

    2012-12-01

    We study the structure of Bogoliubov quasiparticles, bogolons, the fermionic excitations of paired superfluids that arise from fermion (BCS) pairing, including neutral superfluids, superconductors, and paired quantum Hall states. The naive construction of a stationary quasiparticle in which the deformation of the pair field is neglected leads to a contradiction: it carries a net electrical current even though it does not move. However, treating the pair field self-consistently resolves this problem: in a neutral superfluid, a dipolar current pattern is associated with the quasiparticle for which the total current vanishes. When Maxwell electrodynamics is included, as appropriate to a superconductor, this pattern is confined over a penetration depth. For paired quantum Hall states of composite fermions, the Maxwell term is replaced by a Chern-Simons term, which leads to a dipolar charge distribution and consequently to a dipolar current pattern.

  17. Self-consistent calculations of the strength function and radiative neutron capture cross section for stable and unstable tin isotopes

    SciTech Connect

    Avdeenkov, A.; Goriely, S.; Kamerdzhiev, S.; Krewald, S.

    2011-06-15

    The E1 strength function for 15 stable and unstable Sn even-even isotopes from A=100 to A=176 are calculated using a self-consistent microscopic theory which, in addition to the standard (quasiparticle) random-phase approximation [(Q)RPA] approach, takes into account phonon coupling and the single-particle continuum (by means of the discretization procedure) with a cutoff of 100 MeV. Our analysis shows two distinct regions for which the integral characteristics of both the giant and pygmy resonances behave rather differently. For neutron-rich nuclei, starting from {sup 132}Sn, we obtain a giant E1 resonance which significantly deviates from the widely used systematics extrapolated from experimental data in the {beta}-stability valley. We show that the inclusion of phonon coupling is necessary for a proper description of the low-energy pygmy resonances and the corresponding transition densities for A<132 nuclei, while in the A>132 region the influence of phonon coupling is significantly smaller. The radiative neutron capture cross sections leading to the stable {sup 124}Sn and unstable {sup 132}Sn and {sup 150}Sn nuclei are calculated with both the (Q)RPA and the beyond-(Q)RPA strength functions and shown to be sensitive to both the predicted low-lying strength and the phonon-coupling contribution. The comparison with the widely used phenomenological generalized Lorentzian approach shows considerable differences both for the strength function and the radiative neutron capture cross section. In particular, for the neutron-rich {sup 150}Sn, the reaction cross section is found to be increased by a factor greater than 20. We conclude that the present approach may provide a complete and coherent description of the {gamma}-ray-strength function for astrophysics applications. In particular, such calculations are highly recommended for a reliable estimate of the electromagnetic properties of exotic nuclei.

  18. Doubly self-consistent field theory of grafted polymers under simple shear in steady state

    SciTech Connect

    Suo, Tongchuan; Whitmore, Mark D.

    2014-03-21

    We present a generalization of the numerical self-consistent mean-field theory of polymers to the case of grafted polymers under simple shear. The general theoretical framework is presented, and then applied to three different chain models: rods, Gaussian chains, and finitely extensible nonlinear elastic (FENE) chains. The approach is self-consistent at two levels. First, for any flow field, the polymer density profile and effective potential are calculated self-consistently in a manner similar to the usual self-consistent field theory of polymers, except that the calculation is inherently two-dimensional even for a laterally homogeneous system. Second, through the use of a modified Brinkman equation, the flow field and the polymer profile are made self-consistent with respect to each other. For all chain models, we find that reasonable levels of shear cause the chains to tilt, but it has very little effect on the overall thickness of the polymer layer, causing a small decrease for rods, and an increase of no more than a few percent for the Gaussian and FENE chains. Using the FENE model, we also probe the individual bond lengths, bond correlations, and bond angles along the chains, the effects of the shear on them, and the solvent and bonded stress profiles. We find that the approximations needed within the theory for the Brinkman equation affect the bonded stress, but none of the other quantities.

  19. An overview of self-consistent methods for fiber-reinforced composites

    NASA Technical Reports Server (NTRS)

    Gramoll, Kurt C.; Freed, Alan D.; Walker, Kevin P.

    1991-01-01

    The Walker et al. (1989) self-consistent method to predict both the elastic and the inelastic effective material properties of composites is examined and compared with the results of other self-consistent and elastically based solutions. The elastic part of their method is shown to be identical to other self-consistent methods for non-dilute reinforced composite materials; they are the Hill (1965), Budiansky (1965), and Nemat-Nasser et al. (1982) derivations. A simplified form of the non-dilute self-consistent method is also derived. The predicted, elastic, effective material properties for fiber reinforced material using the Walker method was found to deviate from the elasticity solution for the v sub 31, K sub 12, and mu sub 31 material properties (fiber is in the 3 direction) especially at the larger volume fractions. Also, the prediction for the transverse shear modulus, mu sub 12, exceeds one of the accepted Hashin bounds. Only the longitudinal elastic modulus E sub 33 agrees with the elasticity solution. The differences between the Walker and the elasticity solutions are primarily due to the assumption used in the derivation of the self-consistent method, i.e., the strain fields in the inclusions and the matrix are assumed to remain constant, which is not a correct assumption for a high concentration of inclusions.

  20. Quasiparticle electronic structure of bulk and slab Bi2Se3 and Bi2Te3

    NASA Astrophysics Data System (ADS)

    Barker, Bradford; Deslippe, Jack; Yazyev, Oleg; Louie, Steven G.

    2014-03-01

    We present ab initio calculations of the quasiparticle electronic band structure of three-dimensional topological insulator materials Bi2Se3 and Bi2Te3. The mean-field DFT calculation is performed with fully relativistic pseudopotentials, generating spinor wavefunctions in a plane-wave basis. Quasiparticle properties are computed with a one-shot ab initio GW calculation. We use both bulk and slab forms of the materials to better understand the quasiparticle band gaps and Fermi velocities of the topological surface states of these materials. This work was supported by NSF grant No. DMR10-1006184 and U.S. DOE under Contract No. DE-AC02-05CH11231. Computational resources have been provided by DOE at LBNL's NERSC facility and the NSF through XSEDE resources at NICS.

  1. Effective dipole moment for the mode coupling instability: Mapping of self-consistent wake models

    SciTech Connect

    Roecker, T. B.; Zhdanov, S. K.; Ivlev, A. V.; Morfill, G. E.; Lampe, M.; Joyce, G.

    2012-07-15

    The theory of the mode coupling instability operating in two-dimensional plasma crystals is generalized, by employing the linear plasma response formalism to describe the interparticle interactions self-consistently. In this approach, the underlying ion distribution function is calculated from first principles. Subthermal and suprathermal regimes of the ion flow are considered. A mapping procedure is proposed, which relates the self-consistent coupling coefficients to the effective dipole moment of the wake-the parameter which characterizes the mode coupling in the framework of the conventionally used Yukawa/point-wake model. The importance of the self-consistent approach is demonstrated by comparing the theoretically obtained dipole moments with the values deduced from experiments.

  2. Self-consistent dynamical models for early-type galaxies in the CALIFA Survey

    NASA Astrophysics Data System (ADS)

    Posti, L.; van de Ven, G.; Binney, J.; Nipoti, C.; Ciotti, L.

    2016-06-01

    We present the first application of self-consistent, continuous models with distribution functions (DFs) depending on the action integrals, to a sample of nearby early-type galaxies in the CALIFA Survey. Each model is axisymmetric, flattened, anisotropic and rotating and the total gravitational potential is self-consistently generated by the density distribution. The spatially-resolved kinematics of the CALIFA Survey gives solid constraints to the models' parameters: we fit the galaxies' surface brightness and the galaxies' spatially resolved kinematics and we estimate dynamical masses in agreement with other dynamical modelling approaches. For each galaxy, the best model provides an analytic DF which fully characterizes the velocity distribution of the stars. The fact that the DF depends on the action integrals makes it easy to extend the present models to have multiple components, such as bulge, stellar disc and dark and stellar halo, in equilibrium with their self-consistent gravitational potential.

  3. MultiSIMNRA: A computational tool for self-consistent ion beam analysis using SIMNRA

    NASA Astrophysics Data System (ADS)

    Silva, T. F.; Rodrigues, C. L.; Mayer, M.; Moro, M. V.; Trindade, G. F.; Aguirre, F. R.; Added, N.; Rizzutto, M. A.; Tabacniks, M. H.

    2016-03-01

    SIMNRA is widely adopted by the scientific community of ion beam analysis for the simulation and interpretation of nuclear scattering techniques for material characterization. Taking advantage of its recognized reliability and quality of the simulations, we developed a computer program that uses multiple parallel sessions of SIMNRA to perform self-consistent analysis of data obtained by different ion beam techniques or in different experimental conditions of a given sample. In this paper, we present a result using MultiSIMNRA for a self-consistent multi-elemental analysis of a thin film produced by magnetron sputtering. The results demonstrate the potentialities of the self-consistent analysis and its feasibility using MultiSIMNRA.

  4. The generalized active space concept in multiconfigurational self-consistent field methods.

    PubMed

    Ma, Dongxia; Li Manni, Giovanni; Gagliardi, Laura

    2011-07-28

    A multiconfigurational self-consistent field method based on the concept of generalized active space (GAS) is presented. GAS wave functions are obtained by defining an arbitrary number of active spaces with arbitrary occupation constraints. By a suitable choice of the GAS spaces, numerous ineffective configurations present in a large complete active space (CAS) can be removed, while keeping the important ones in the CI space. As a consequence, the GAS self-consistent field approach retains the accuracy of the CAS self-consistent field (CASSCF) ansatz and, at the same time, can deal with larger active spaces, which would be unaffordable at the CASSCF level. Test calculations on the Gd atom, Gd(2) molecule, and oxoMn(salen) complex are presented. They show that GAS wave functions achieve the same accuracy as CAS wave functions on systems that would be prohibitive at the CAS level. PMID:21806111

  5. Statistical dynamics of classical systems: A self-consistent field approach

    SciTech Connect

    Grzetic, Douglas J. Wickham, Robert A.; Shi, An-Chang

    2014-06-28

    We develop a self-consistent field theory for particle dynamics by extremizing the functional integral representation of a microscopic Langevin equation with respect to the collective fields. Although our approach is general, here we formulate it in the context of polymer dynamics to highlight satisfying formal analogies with equilibrium self-consistent field theory. An exact treatment of the dynamics of a single chain in a mean force field emerges naturally via a functional Smoluchowski equation, while the time-dependent monomer density and mean force field are determined self-consistently. As a simple initial demonstration of the theory, leaving an application to polymer dynamics for future work, we examine the dynamics of trapped interacting Brownian particles. For binary particle mixtures, we observe the kinetics of phase separation.

  6. Ab initio quasiparticle bandstructure of ABA and ABC-stacked graphene trilayers

    NASA Astrophysics Data System (ADS)

    Menezes, Marcos; Capaz, Rodrigo; Louie, Steven

    2013-03-01

    We obtain the quasiparticle band structure of ABA and ABC-stacked graphene trilayers through ab initio density functional theory (DFT) and many-body quasiparticle calculations within the GW approximation. To interpret our results, we fit the DFT and GW π bands to a low energy tight-binding model, which is found to reproduce very well the observed features near the K point. The values of the extracted hopping parameters are reported and compared with available theoretical and experimental data. For both stackings, the quasiparticle corrections lead to a renormalization of the Fermi velocity, an effect also observed in previous calculations on monolayer graphene. They also increase the separation between the higher energy bands, which is proportional to the nearest neighbor interlayer hopping parameter γ1. Both features are brought to closer agreement with experiment through the quasiparticle corrections. Finally, other effects, such as trigonal warping, electron-hole assymetry and energy gaps are discussed in terms of the associated parameters. This work was supported by the Brazilian funding agencies: CAPES, CNPq, FAPERJ and INCT-Nanomateriais de Carbono. It was also supported by NSF grant No. DMR10-1006184 and U.S. DOE under Contract No. DE-AC02-05CH11231.

  7. Thermodynamically self-consistent non-stochastic micromagnetic model for the ferromagnetic state

    SciTech Connect

    Dvornik, Mykola Vansteenkiste, Arne; Van Waeyenberge, Bartel

    2014-10-20

    In this work, a self-consistent thermodynamic approach to micromagnetism is presented. The magnetic degrees of freedom are modeled using the Landau-Lifshitz-Baryakhtar theory, which separates the different contributions to the magnetic damping, and thereby allows them to be coupled to the electron and phonon systems in a self-consistent way. We show that this model can quantitatively reproduce ultrafast magnetization dynamics in Nickel suggesting that in ferromagnetic metals the ultrafast angular momentum transfer happens via the relativistic spin-electron scattering.

  8. A self-consistent theory of collective alpha particle losses induced by Alfvenic turbulence

    SciTech Connect

    Biglari, H.; Diamond, P.H.

    1992-01-01

    The nonlinear dynamics of kinetic Alfven waves, resonantly excited by energetic ions/alpha particles, is investigated. It is shown that {alpha}-particles govern both linear instability and nonlinear saturation dynamics, while the background MHD turbulence results only in a nonlinear real frequency shift. The most efficient saturation mechanism is found to be self-induced profile modification. Expressions for the fluctuation amplitudes and the {alpha}-particle radial flux are self-consistently derived. The work represents the first self-consistent, turbulent treatment of collective {alpha}-particle losses by Alfvenic fluctuations.

  9. Self-consistent RPA based on a many-body vacuum

    SciTech Connect

    Jemaie, M.; Schuck, P.

    2011-08-15

    Self-Consistent RPA is extended in a way so that it is compatible with a variational ansatz for the ground-state wave function as a fermionic many-body vacuum. Employing the usual equation-of-motion technique, we arrive at extended RPA equations of the Self-Consistent RPA structure. In principle the Pauli principle is, therefore, fully respected. However, the correlation functions entering the RPA matrix can only be obtained from a systematic expansion in powers of some combinations of RPA amplitudes. We demonstrate for a model case that this expansion may converge rapidly.

  10. Self-consistent descriptions of vector mesons in hot matter reexamined

    SciTech Connect

    Riek, Felix; Knoll, Joern

    2010-08-15

    Technical concepts are presented that improve the self-consistent treatment of vector mesons in a hot and dense medium. First applications concern an interacting gas of pions and {rho} mesons. As an extension of earlier studies, we thereby include random-phase-approximation-type vertex corrections and further use dispersion relations to calculate the real part of the vector-meson self-energy. An improved projection method preserves the four transversality of the vector-meson polarization tensor throughout the self-consistent calculations, thereby keeping the scheme void of kinematical singularities.

  11. Efficient self-consistent DFT calculation of nondynamic correlation based on the B05 method

    PubMed Central

    Proynov, Emil; Shao, Yihan; Kong, Jing

    2010-01-01

    Becke’s B05 method of describing nondynamic electron correlation in Density Functional Theory is implemented self-consistently with computational efficiency. Important modifications of the method are proposed in order to make the self-consistency feasible. Resolution-of-identity technique is used to reduce dramatically the cost of the required exact-exchange energy density. The method is briefly validated on a variety of properties. It describes accurately for the first time the subtle energetics of the NO dimer, an exemplary system of strong nondynamic correlation. The efficient algorithm for the exact-exchange energy density can be applied to other functionals that use this quantity. PMID:20640046

  12. Characterisation of gunshot residue particles using self-consistent ion beam analysis

    NASA Astrophysics Data System (ADS)

    Bailey, M. J.; Jeynes, C.

    2009-06-01

    Individual particles of gunshot residue were studied with particle-induced X-ray emission and backscattering spectrometry using a 2.5 MeV H + beam focussed to ˜4 μm and self-consistent fitting of the data. The geometry of these spherical particles was considered in order to accurately fit the corresponding particle spectrum and therefore to quantify the trace element composition of these particles. The demonstrable self-consistency of this method allows the compositions of most residue particles to be determined unambiguously and with a higher sensitivity to trace elements than conventional methods.

  13. The Progenitor of GW150914

    NASA Astrophysics Data System (ADS)

    Woosley, S. E.

    2016-06-01

    The spectacular detection of gravitational waves (GWs) from GW150914 and its reported association with a gamma-ray burst (GRB) offer new insights into the evolution of massive stars. Here, it is shown that no single star of any mass and credible metallicity is likely to produce the observed GW signal. Stars with helium cores in the mass range 35–133 M ⊙ encounter the pair instability and either explode or pulse until the core mass is less than 45 M ⊙, smaller than the combined mass of the observed black holes. The rotation of more massive helium cores is either braked by interaction with a slowly rotating hydrogen envelope, if one is present, or by mass loss, if one is not. The very short interval between the GW signal and the observed onset of the putative GRB in GW150914 is also too short to have come from a single star. A more probable model for making the gravitational radiation is the delayed merger of two black holes made by 70 and 90 M ⊙ stars in a binary system. The more massive component was a pulsational-pair instability supernova before making the first black hole.

  14. Importance of semicore states in GW calculations for simulating accurately the photoemission spectra of metal phthalocyanine molecules

    NASA Astrophysics Data System (ADS)

    Umari, P.; Fabris, S.

    2012-05-01

    The quasi-particle energy levels of the Zn-Phthalocyanine (ZnPc) molecule calculated with the GW approximation are shown to depend sensitively on the explicit description of the metal-center semicore states. We find that the calculated GW energy levels are in good agreement with the measured experimental photoemission spectra only when explicitly including the Zn 3s and 3p semicore states in the valence. The main origin of this effect is traced back to the exchange term in the self-energy GW approximation. Based on this finding, we propose a simplified approach for correcting GW calculations of metal phthalocyanine molecules that avoids the time-consuming explicit treatment of the metal semicore states. Our method allows for speeding up the calculations without compromising the accuracy of the computed spectra.

  15. Condensing Non-Abelian Quasiparticles

    SciTech Connect

    Hermanns, M.

    2010-02-05

    A most interesting feature of certain fractional quantum Hall states is that their quasiparticles obey non-Abelian fractional statistics. So far, candidate non-Abelian wave functions have been constructed from conformal blocks in cleverly chosen conformal field theories. In this work we present a hierarchy scheme by which we can construct daughter states by condensing non-Abelian quasiparticles (as opposed to quasiholes) in a parent state, and show that the daughters have a non-Abelian statistics that differs from the parent. In particular, we discuss the daughter of the bosonic, spin-polarized Moore-Read state at nu=4/3 as an explicit example.

  16. Self-consistent Gödel cosmology with spin-density in Riemann-Cartan spacetime

    NASA Astrophysics Data System (ADS)

    Smalley, Larry L.

    1986-01-01

    We show that the Gödel metric for a rotating cosmology is compatible with the self-consistent formulation of the Einstein-Cartan metric-torsion theory for a spinning fluid. Fruitful discussions with J.R. Ray and F.W. Hehl are greatfully acknowledged. I wish to also thank J. Bradas who brought a reference to my attention.

  17. Subjective Confidence in Perceptual Judgments: A Test of the Self-Consistency Model

    ERIC Educational Resources Information Center

    Koriat, Asher

    2011-01-01

    Two questions about subjective confidence in perceptual judgments are examined: the bases for these judgments and the reasons for their accuracy. Confidence in perceptual judgments has been claimed to rest on qualitatively different processes than confidence in memory tasks. However, predictions from a self-consistency model (SCM), which had been…

  18. Self-consistency tests of large-scale dynamics parameterizations for single-column modeling

    SciTech Connect

    Edman, Jacob P.; Romps, David M.

    2015-03-18

    Large-scale dynamics parameterizations are tested numerically in cloud-resolving simulations, including a new version of the weak-pressure-gradient approximation (WPG) introduced by Edman and Romps (2014), the weak-temperature-gradient approximation (WTG), and a prior implementation of WPG. We perform a series of self-consistency tests with each large-scale dynamics parameterization, in which we compare the result of a cloud-resolving simulation coupled to WTG or WPG with an otherwise identical simulation with prescribed large-scale convergence. In self-consistency tests based on radiative-convective equilibrium (RCE; i.e., no large-scale convergence), we find that simulations either weakly coupled or strongly coupled to either WPG or WTG are self-consistent, but WPG-coupled simulations exhibit a nonmonotonic behavior as the strength of the coupling to WPG is varied. We also perform self-consistency tests based on observed forcings from two observational campaigns: the Tropical Warm Pool International Cloud Experiment (TWP-ICE) and the ARM Southern Great Plains (SGP) Summer 1995 IOP. In these tests, we show that the new version of WPG improves upon prior versions of WPG by eliminating a potentially troublesome gravity-wave resonance.

  19. Enriching Elementary Quantum Mechanics with the Computer: Self-Consistent Field Problems in One Dimension

    ERIC Educational Resources Information Center

    Bolemon, Jay S.; Etzold, David J.

    1974-01-01

    Discusses the use of a small computer to solve self-consistent field problems of one-dimensional systems of two or more interacting particles in an elementary quantum mechanics course. Indicates that the calculation can serve as a useful introduction to the iterative technique. (CC)

  20. Self-Consistent Non-Stationary Theory of Multipactor in DLA Structures

    SciTech Connect

    Sinitsyn, O. V.; Nusinovich, G. S.; Antonsen, T. M.; Kishek, R.

    2009-01-22

    In this paper a non-stationary self-consistent theoretical model of multipactor in dielectric loaded accelerator structures is proposed. In comparison with our previous work, the effects of the cylindricity are included. The corresponding numerical implementation of the model is described and some simulation results are shown.

  1. Remembering and telling self-consistent and self-discrepant memories.

    PubMed

    Mutlutürk, Aysu; Tekcan, Ali I

    2016-04-01

    It has been argued that memories that are inconsistent with one's self would differ from those that are consistent with the self. The present study addresses retrieval, phenomenology, rehearsal and narrative characteristics of autobiographical memories that are consistent versus discrepant with one's self. One hundred participants were asked to recall one self-consistent and one self-discrepant memory as well as an episode of telling these memories to others. They also filled out the Autobiographical Memory Questionnaire and the Centrality of Event Scale for each memory. Results showed no difference between self-consistent and self-discrepant memories in retrieval time, specificity or phenomenology. However, self-discrepant memory narratives contained more meaning-making statements and less autonomy than self-consistent memories. Compared to self-consistent memories, self-discrepant memories were told to fewer people, and listener responses were more negative when they were told. Results are discussed in relation to the functions these memories serve. PMID:25785773

  2. Self-consistently improved finite temperature effective potential for gauge theories

    SciTech Connect

    Amelino-Camelia, G. )

    1994-03-15

    The finite temperature effective potential of the Abelian Higgs model is studied using the self-consistent composite operator method, which can be used to sum up the contributions of daisy and superdaisy diagrams. The effect of the momentum dependence of the effective masses is estimated by using a Rayleigh-Ritz variational approximation.

  3. Self-consistency tests of large-scale dynamics parameterizations for single-column modeling

    DOE PAGESBeta

    Edman, Jacob P.; Romps, David M.

    2015-03-18

    Large-scale dynamics parameterizations are tested numerically in cloud-resolving simulations, including a new version of the weak-pressure-gradient approximation (WPG) introduced by Edman and Romps (2014), the weak-temperature-gradient approximation (WTG), and a prior implementation of WPG. We perform a series of self-consistency tests with each large-scale dynamics parameterization, in which we compare the result of a cloud-resolving simulation coupled to WTG or WPG with an otherwise identical simulation with prescribed large-scale convergence. In self-consistency tests based on radiative-convective equilibrium (RCE; i.e., no large-scale convergence), we find that simulations either weakly coupled or strongly coupled to either WPG or WTG are self-consistent, butmore » WPG-coupled simulations exhibit a nonmonotonic behavior as the strength of the coupling to WPG is varied. We also perform self-consistency tests based on observed forcings from two observational campaigns: the Tropical Warm Pool International Cloud Experiment (TWP-ICE) and the ARM Southern Great Plains (SGP) Summer 1995 IOP. In these tests, we show that the new version of WPG improves upon prior versions of WPG by eliminating a potentially troublesome gravity-wave resonance.« less

  4. Moment-based, self-consistent linear analysis of gyrotron oscillators

    NASA Astrophysics Data System (ADS)

    Braunmueller, F.; Tran, T. M.; Alberti, S.; Hogge, J.-Ph.; Tran, M. Q.

    2014-04-01

    A new model for simulating gyrotron oscillators in the monomode time-dependent linear self-consistent regime is presented. Starting from a nonlinear time-dependent monomode model, the linearization and the following simplification of the model, based on a moment approach, are described. This simplified model represents a numerically efficient model and allows to have a deeper physical insight, in particular, for regimes dominated by self-consistent effects such as for the gyro-backward wave instability. One specific case of a gyrotron cavity is studied in detail and compared with experimental results, with special attention to self-consistent effects and to the differences with a model using a fixed field profile. Self-consistent linear simulations are, amongst other applications, important for the design of frequency-tunable gyrotrons or high-power gyrotrons with cavities having a relatively low quality factor, but also for studies of parasitic oscillations as they may occur in beam ducts and/or in the launcher section following the interaction cavity.

  5. Moment-based, self-consistent linear analysis of gyrotron oscillators

    SciTech Connect

    Braunmueller, F. Tran, T. M.; Alberti, S.; Hogge, J.-Ph.; Tran, M. Q.

    2014-04-15

    A new model for simulating gyrotron oscillators in the monomode time-dependent linear self-consistent regime is presented. Starting from a nonlinear time-dependent monomode model, the linearization and the following simplification of the model, based on a moment approach, are described. This simplified model represents a numerically efficient model and allows to have a deeper physical insight, in particular, for regimes dominated by self-consistent effects such as for the gyro-backward wave instability. One specific case of a gyrotron cavity is studied in detail and compared with experimental results, with special attention to self-consistent effects and to the differences with a model using a fixed field profile. Self-consistent linear simulations are, amongst other applications, important for the design of frequency-tunable gyrotrons or high-power gyrotrons with cavities having a relatively low quality factor, but also for studies of parasitic oscillations as they may occur in beam ducts and/or in the launcher section following the interaction cavity.

  6. Self-consistent chaos in a mean-field Hamiltonian model of fluids and plasmas

    NASA Astrophysics Data System (ADS)

    del-Castillo-Negrete, D.; Firpo, Marie-Christine

    2002-11-01

    We present a mean-field Hamiltonian model that describes the collective dynamics of marginally stable fluids and plasmas. In plasmas, the model describes the self-consistent evolution of electron holes and clumps in phase space. In fluids, the model describes the dynamics of vortices with negative and positive circulation in shear flows. The mean-field nature of the system makes it a tractable model to study the dynamics of large degrees-of-freedom, coupled Hamiltonian systems. Here we focus in the role of self-consistent chaos in the formation and destruction of phase space coherent structures. Numerical simulations in the finite N and in the Narrow kinetic limit (where N is the number of particles) show the existence of coherent, rotating dipole states. We approximate the dipole as two macroparticles, and show that the N = 2 limit has a family of rotating integrable solutions described by a one degree-of-freedom nontwist Hamiltonian. The coherence of the dipole is explained in terms of a parametric resonance between the rotation frequency of the macroparticles and the oscillation frequency of the self-consistent mean field. For a class of initial conditions, the mean field exhibits a self-consistent, elliptic-hyperbolic bifurcation that leads to the destruction of the dipole and violent mixing of the phase space.

  7. The GW electronic structure of cubic RbMF3 perovskites (M = Be, Mg, Ca, Sr, Ba)

    NASA Astrophysics Data System (ADS)

    Syrotyuk, Stepan V.; Shved, Vira M.

    2015-09-01

    The electronic energy band spectra of cubic RbMF3 perovskites (M = Be, Mg, Ca, Sr, Ba) have been evaluated within the projector augmented waves (PAW) approach by means of the ABINIT code. The Kohn-Sham single-particle states have been found in the LDA framework. Our parameters of the electron energy bands obtained in the LDA are in good agreement with the published results of other authors. The calculated GW band gaps and dielectric constants are well compared with the available experimental data. The quasiparticle energies and density of states of electrons as well as the dielectric constants were obtained in the approximation GW for the first time.

  8. Using Hartree-Fock pseudopotentials in GW calculations

    NASA Astrophysics Data System (ADS)

    Hamann, D. R.; Vanderbilt, David

    2010-03-01

    The issue of including shallow ``semi-core'' states as valence has recently resurfaced in the context of self-consistent GW calculations.footnotetextF. Bruneval et al., Phys. Rev. Lett. 97, 267601 (2006). Supposing that semi-core-valence exchange is the dominant process necessitating the inclusion of semi-cores, we have investigated whether the use Hartree-Fock pseudopotentialsfootnotetextW. A. Al-Saidi, E. J. Walter, and A. M. Rappe, Phys. Rev. B 77, 075122 (2008). instead of density-functional psp's might obviate the need for semi-cores. The answers to this question appear to be ``yes'' for the case of CuCl (filled d shell), and ``semi-cores don't matter anyway'' for ScN (empty d shell). Opportunity permitting, additional examples will be discussed.

  9. Predictive GW calculations using plane waves and pseudopotentials

    NASA Astrophysics Data System (ADS)

    Klimeš, Jiří; Kaltak, Merzuk; Kresse, Georg

    2014-08-01

    We derive a finite-basis-set correction for quasiparticle (QP) energies in the GW approximation and many-body correlation energies in the random phase approximation. Since the correction requires only knowledge of the ground-state density distribution, it is straightforward to implement in any plane-wave code and significantly improves convergence at negligible computational cost. The expression also indicates that QP energies might converge to the wrong value using the projector augmented wave (PAW) method since the overlap densities of occupied orbitals and high-energy, plane-wave-like orbitals are inaccurately described. The error is shown to be related to the incompleteness of the partial waves inside the atomic spheres. It can be avoided by adopting norm-conserving partial waves. G0W0 and GW0 results based on such norm-conserving PAW potentials are presented for a large set of semiconductors and insulators. Accurate extrapolation procedures to the infinite-basis-set limit and infinite-k-point limit are discussed in detail.

  10. Quasiparticle Interactions in Neutron Matter for Applications in Neutron Stars

    NASA Technical Reports Server (NTRS)

    Wambach, J; Ainsworth, T. L.; Pines, D.

    1993-01-01

    A microscopic model for the quasiparticle interaction in neutron matter is presented. Both-particle (pp) and particle-hole (ph) correlations are included. The pp correlations are treated in semi-empirical way, while ph correlations are incorporated by solving coupled two-body equations for particle-hole interaction and the scattering amplitude of the Fermi sphere. The resulting integral equations self-consistently sum the ph reducible diagrams. Antisymmetry is kept at all stages and hence the forward-scattering sum rules for the scattering amplitude are obeyed. Results for Landau parameters and transport coefficients in a density regime representing the crust of a neutron star are presented. We also estimate the (1)S(sub 0) gap parameter for neutron superfluidity and comment briefly on neutron-star implications.

  11. Quasiparticle and optical properties of rutile and anatase TiO2

    NASA Astrophysics Data System (ADS)

    Kang, Wei; Hybertsen, Mark S.

    2010-08-01

    Quasiparticle excitation energies and optical properties of TiO2 in the rutile and anatase structures are calculated using many-body perturbation-theory methods. Calculations are performed for a frozen crystal lattice; electron-phonon coupling is not explicitly considered. In the GW method, several approximations are compared and it is found that inclusion of the full frequency dependence as well as explicit treatment of the Ti semicore states are essential for accurate calculation of the quasiparticle energy-band gap. The calculated quasiparticle energies are in good agreement with available photoemission and inverse photoemission experiments. The results of the GW calculations, together with the calculated static screened Coulomb interaction, are utilized in the Bethe-Salpeter equation to calculate the dielectric function γ2(ω) for both the rutile and anatase structures. The results are in good agreement with experimental observations, particularly the onset of the main absorption features around 4 eV. For comparison to low-temperature optical-absorption measurements that resolve individual excitonic transitions in rutile, the low-lying discrete excitonic energy levels are calculated with electronic screening only. The lowest energy exciton found in the energy gap of rutile has a binding energy of 0.13 eV. In agreement with experiment, it is not dipole allowed but the calculated exciton energy exceeds that measured in absorption experiments by about 0.22 eV and the scale of the exciton binding energy is also too large. The quasiparticle energy alignment of rutile is calculated for nonpolar (110) surfaces. In the GW approximation, the valence-band maximum is 7.8 eV below the vacuum level, showing a small shift from density-functional theory results.

  12. Self-consistent solution for proximity effect and Josephson current in ballistic graphene SNS Josephson junctions

    SciTech Connect

    Black-Schaffer, Annica M.

    2010-04-06

    We use a tight-binding Bogoliubov-de Gennes (BdG) formalism to self-consistently calculate the proximity effect, Josephson current, and local density of states in ballistic graphene SNS Josephson junctions. Both short and long junctions, with respect to the superconducting coherence length, are considered, as well as different doping levels of the graphene. We show that self-consistency does not notably change the current-phase relationship derived earlier for short junctions using the non-selfconsistent Dirac-BdG formalism but predict a significantly increased critical current with a stronger junction length dependence. In addition, we show that in junctions with no Fermi level mismatch between the N and S regions superconductivity persists even in the longest junctions we can investigate, indicating a diverging Ginzburg-Landau superconducting coherence length in the normal region.

  13. Self-consistent dual boson approach to single-particle and collective excitations in correlated systems

    NASA Astrophysics Data System (ADS)

    Stepanov, E. A.; van Loon, E. G. C. P.; Katanin, A. A.; Lichtenstein, A. I.; Katsnelson, M. I.; Rubtsov, A. N.

    2016-01-01

    We propose an efficient dual boson scheme, which extends the dynamical mean-field theory paradigm to collective excitations in correlated systems. The theory is fully self-consistent both on the one- and on the two-particle level, thus describing the formation of collective modes as well as the renormalization of electronic and bosonic spectra on equal footing. The method employs an effective impurity model comprising both fermionic and bosonic hybridization functions. Only single- and two-electron Green's functions of the reference problem enter the theory, due to the optimal choice of the self-consistency condition for the effective bosonic bath. We show that the theory is naturally described by a dual Luttinger-Ward functional and obeys the relevant conservation laws.

  14. Self-Consistent Sources for Integrable Equations Via Deformations of Binary Darboux Transformations

    NASA Astrophysics Data System (ADS)

    Chvartatskyi, Oleksandr; Dimakis, Aristophanes; Müller-Hoissen, Folkert

    2016-06-01

    We reveal the origin and structure of self-consistent source extensions of integrable equations from the perspective of binary Darboux transformations. They arise via a deformation of the potential that is central in this method. As examples, we obtain in particular matrix versions of self-consistent source extensions of the KdV, Boussinesq, sine-Gordon, nonlinear Schrödinger, KP, Davey-Stewartson, two-dimensional Toda lattice and discrete KP equation. We also recover a (2+1)-dimensional version of the Yajima-Oikawa system from a deformation of the pKP hierarchy. By construction, these systems are accompanied by a hetero binary Darboux transformation, which generates solutions of such a system from a solution of the source-free system and additionally solutions of an associated linear system and its adjoint. The essence of all this is encoded in universal equations in the framework of bidifferential calculus.

  15. Communication: A difference density picture for the self-consistent field ansatz

    NASA Astrophysics Data System (ADS)

    Parrish, Robert M.; Liu, Fang; Martínez, Todd J.

    2016-04-01

    We formulate self-consistent field (SCF) theory in terms of an interaction picture where the working variable is the difference density matrix between the true system and a corresponding superposition of atomic densities. As the difference density matrix directly represents the electronic deformations inherent in chemical bonding, this "difference self-consistent field (dSCF)" picture provides a number of significant conceptual and computational advantages. We show that this allows for a stable and efficient dSCF iterative procedure with wholly single-precision Coulomb and exchange matrix builds. We also show that the dSCF iterative procedure can be performed with aggressive screening of the pair space. These approximations are tested and found to be accurate for systems with up to 1860 atoms and >10 000 basis functions, providing for immediate overall speedups of up to 70% in the heavily optimized TeraChem SCF implementation.

  16. Self-consistent simulation studies of periodically focused intense charged-particle beams

    NASA Astrophysics Data System (ADS)

    Chen, C.; Jameson, R. A.

    1995-09-01

    A self-consistent two-dimensional model is used to investigate intense charged-particle beam propagation through a periodic solenoidal focusing channel, particularly in the regime in which there is a mismatch between the beam and the focusing channel. The present self-consistent studies confirm that mismatched beams exhibit nonlinear resonances and chaotic behavior in the envelope evolution, as predicted by an earlier envelope analysis [C. Chen and R. C. Davidson, Phys. Rev. Lett. 72, 2195 (1994)]. Transient effects due to emittance growth are studied, and halo formation is investigated. The halo size is estimated. The halo characteristics for a periodic focusing channel are found to be qualitatively the same as those for a uniform focusing channel. A threshold condition is obtained numerically for halo formation in mismatched beams in a uniform focusing channel, which indicates that relative envelope mismatch must be kept well below 20% to prevent space-charge-dominated beams from developing halos.

  17. Self-Consistent Sources for Integrable Equations Via Deformations of Binary Darboux Transformations

    NASA Astrophysics Data System (ADS)

    Chvartatskyi, Oleksandr; Dimakis, Aristophanes; Müller-Hoissen, Folkert

    2016-08-01

    We reveal the origin and structure of self-consistent source extensions of integrable equations from the perspective of binary Darboux transformations. They arise via a deformation of the potential that is central in this method. As examples, we obtain in particular matrix versions of self-consistent source extensions of the KdV, Boussinesq, sine-Gordon, nonlinear Schrödinger, KP, Davey-Stewartson, two-dimensional Toda lattice and discrete KP equation. We also recover a (2+1)-dimensional version of the Yajima-Oikawa system from a deformation of the pKP hierarchy. By construction, these systems are accompanied by a hetero binary Darboux transformation, which generates solutions of such a system from a solution of the source-free system and additionally solutions of an associated linear system and its adjoint. The essence of all this is encoded in universal equations in the framework of bidifferential calculus.

  18. Thermodynamic approach to the interpretation of self-consistent pressure profiles in a tokamak

    SciTech Connect

    Dyabilin, K. S.; Razumova, K. A.

    2015-09-15

    The phenomenon of invariable pressure profiles in tokamaks is interpreted in the framework of the thermodynamic approach suggesting that invariable self-consistent states correspond to the minimum of free energy. Solutions qualitatively consistent with the experiment are obtained under the assumption that the mechanism for the formation of self-consistent profiles is directly related to equilibrium diamagnetic currents. The dynamics of the system and specific transport phenomena, such as energy and particle pinching and a decrease in the local density under auxiliary electron cyclotron resonance heating (density pump-out), are analyzed in the vicinity of an equilibrium state characterized by a stable pressure profile. The scaling for the energy confinement time deduced from the transport model agrees qualitatively with the ITER scaling based on the analysis of experimental data obtained in many tokamaks. The possibility of using generalized Tsallis statistics to analyze pressure profiles is considered.

  19. Electron Pumping under Non-Markovian Dissipation: The Role of the Self-Consistent Field

    NASA Astrophysics Data System (ADS)

    Grossmann, Frank; Sakurai, Atsunori; Tanimura, Yoshitaka

    2016-03-01

    Focusing on electron transport through a periodically driven resonant tunneling diode, we study the generation of a non-vanishing dc-current by applying symmetry breaking external ac fields with phase difference φ in a statically unbiased system. The effect of an environment is investigated using the system-bath Hamiltonian represented by the electron system coupled to harmonic oscillator modes with a Drude-Lorentz spectral density. To carry out simulations, we use the hierarchal equations of motion approach in the Wigner representation including a self-consistently constructed electric field that is determined from the electron distribution using the Poisson equation. We show that the maximal pumping current at a phase difference near φ = π/2 is strongly influenced by the system-bath coupling strength. The effect of dissipation is diminished if the self-consistent part of the potential is ignored.

  20. Integrated fusion simulation with self-consistent core-pedestal coupling

    NASA Astrophysics Data System (ADS)

    Meneghini, O.; Snyder, P. B.; Smith, S. P.; Candy, J.; Staebler, G. M.; Belli, E. A.; Lao, L. L.; Park, J. M.; Green, D. L.; Elwasif, W.; Grierson, B. A.; Holland, C.

    2016-04-01

    Accurate prediction of fusion performance in present and future tokamaks requires taking into account the strong interplay between core transport, pedestal structure, current profile, and plasma equilibrium. An integrated modeling workflow capable of calculating the steady-state self-consistent solution to this strongly coupled problem has been developed. The workflow leverages state-of-the-art components for collisional and turbulent core transport, equilibrium and pedestal stability. Testing against a DIII-D discharge shows that the workflow is capable of robustly predicting the kinetic profiles (electron and ion temperature and electron density) from the axis to the separatrix in a good agreement with the experiments. An example application is presented, showing self-consistent optimization for the fusion performance of the 15 MA D-T ITER baseline scenario as functions of the pedestal density and ion effective charge Zeff .

  1. Self-consistent axial modeling of surface-wave-produced discharges at low and intermediate pressures.

    PubMed

    Petrova, T; Benova, E; Petrov, G; Zhelyazkov, I

    1999-07-01

    A model for description of the axial structure of a surface-wave-produced and -sustained plasma based on numerical calculation of a complete set of electrodynamic and kinetic equations is presented. The model includes a self-consistent solution to the electron Boltzmann equation, a set of particle balance equations for electrons, excited atoms, atomic and molecular ions, as well as Maxwell's equations with appropriate boundary conditions. A gas thermal balance equation is used to predict the neutral gas temperature self-consistently. Precise calculations of discharge characteristics of an argon plasma column sustained by an azimuthally symmetric surface wave at low and intermediate gas pressures have been performed. A comparison with available experimental data is done in order to test the validity of the model. PMID:11969832

  2. Self-consistent calculation of the electron structure and x-ray spectra of chromium nitride

    SciTech Connect

    Bekenev, V.L.; Lisenko, A.A.; Zhurakovskii, E.A.

    1986-02-01

    The authors calculate the energy band structure of cubic chromium nitride by the self-consistent method of associated plane waves for a broad energy range. Self-consistency led to overlapping of the p-band of nitrogen and the d-band of chromium and to the appearance of an energy discontinuity in the region of unbounded states. The total and local partial densities of the states are calculated. With allowance for the probability of transition, the KB/sub 5/ and L/sub 111/ -emission bands of chromium, the Ka -band of nitrogen, and the K-edge of absorption of chromium in chromium nitride are calculated in a dipole approximation. The possibility of calculating the absorption edge with allowance for the effect of shell holes is discussed. Satisfactory agreement is obtained with experimental data.

  3. Globular structures of a helix-coil copolymer: self-consistent treatment.

    PubMed

    Nowak, C; Rostiashvili, V G; Vilgis, T A

    2007-01-21

    A self-consistent-field theory was developed in the grand canonical ensemble formulation to study transitions in a helix-coil multiblock globule. Helical and coil parts are treated as stiff rods and self-avoiding walks of variable lengths correspondingly. The resulting field theory takes, in addition to the conventional Zimm-Bragg, [J. Chem. Phys. 31, 526 (1959)] parameters, also three-dimensional interaction terms into account. The appropriate differential equations which determine the self-consistent fields were solved numerically with finite element method. Three different phase states are found: open chain, amorphous globule, and nematic liquid-crystalline (LC) globule. The LC-globule formation is driven by the interplay between the hydrophobic helical segment attraction and the anisotropic globule surface energy of an entropic nature. The full phase diagram of the helix-coil copolymer was calculated and thoroughly discussed. The suggested theory shows a clear interplay between secondary and tertiary structures in globular homopolypeptides. PMID:17249898

  4. Dynamics and Self-consistent Chaos in a Mean Field Hamiltonian Model

    NASA Astrophysics Data System (ADS)

    del-Castillo-Negrete, Diego

    We study a mean field Hamiltonian model that describes the collective dynamics of marginally stable fluids and plasmas in the finite N and N-> infty kinetic limit (where N is the number of particles). The linear stability of equilibria in the kinetic model is studied as well as the initial value problem including Landau damping . Numerical simulations show the existence of coherent, rotating dipole states. We approximate the dipole as two macroparticles and show that the N=2 limit has a family of rotating integrable solutions that provide an accurate description of the dynamics. We discuss the role of self-consistent Hamiltonian chaos in the formation of coherent structures, and discuss a mechanism of "violent" mixing caused by a self-consistent elliptic-hyperbolic bifurcation in phase space.

  5. Self-consistent simulation of radio frequency multipactor on micro-grooved dielectric surface

    SciTech Connect

    Cai, Libing; Wang, Jianguo; Cheng, Guoxin; Zhu, Xiangqin; Xia, Hongfu

    2015-02-07

    The multipactor plays a key role in the surface breakdown on the feed dielectric window irradiated by high power microwave. To study the suppression of multipactor, a 2D electrostatic PIC-MCC simulation code was developed. The space charge field, including surface deposited charge and multipactor electron charge field, is obtained by solving 2D Poisson's equation in time. Therefore, the simulation is self-consistent and does not require presetting a fixed space charge field. By using this code, the self-consistent simulation of the RF multipactor on the periodic micro-grooved dielectric surface is realized. The 2D space distributions of the multipactor electrons and space charge field are presented. From the simulation results, it can be found that only half slopes have multipactor discharge when the slope angle exceeds a certain value, and the groove presents a pronounced suppression effect on the multipactor.

  6. Integrated fusion simulation with self-consistent core-pedestal coupling

    SciTech Connect

    Meneghini, Orso; Snyder, P. B.; Smith, S. P.; Candy, J.; Staebler, G. M.; Belli, E. A.; Lao, L. L.; Park, J. M.; Green, David L; Elwasif, Wael R; Grierson, Brian A.; Holland, C.

    2016-01-01

    Accurate prediction of fusion performance in present and future tokamaks requires taking into account the strong interplay between core transport, pedestal structure, current profile and plasma equilibrium. An integrated modeling workflow capable of calculating the steady-state self- consistent solution to this strongly-coupled problem has been developed. The workflow leverages state-of-the-art components for collisional and turbulent core transport, equilibrium and pedestal stability. Validation against DIII-D discharges shows that the workflow is capable of robustly pre- dicting the kinetic profiles (electron and ion temperature and electron density) from the axis to the separatrix in good agreement with the experiments. An example application is presented, showing self-consistent optimization for the fusion performance of the 15 MA D-T ITER baseline scenario as functions of the pedestal density and ion effective charge Z eff.

  7. Criterion for DNA melting in the mean-field modified self-consistent phonon theory

    NASA Astrophysics Data System (ADS)

    Feng, Y.; Prohofsky, E. W.

    1991-04-01

    We have examined the validity of the first-order-perturbation method in calculating eigenfunctions and the criterion for helix melting of mean-field polymers in the modified self-consistent phonon approach (MSPA) theory. It is found that the instability in the self-consistent solution is due to the breakdown of the first-order perturbation. The instability as a criterion for helix melting is therefore techniquely inappropriate. However, the breakdown of the perturbation is due to facts that are directly related to the onset of softening. Previously predicted melting temperatures for various sequence DNA polymers may still represent good estimates to the actual melting temperatures. An alternative criterion is required to define the melting temperature of the polymer DNA double helix in the MSPA theory.

  8. Thermodynamic approach to the interpretation of self-consistent pressure profiles in a tokamak

    NASA Astrophysics Data System (ADS)

    Dyabilin, K. S.; Razumova, K. A.

    2015-09-01

    The phenomenon of invariable pressure profiles in tokamaks is interpreted in the framework of the thermodynamic approach suggesting that invariable self-consistent states correspond to the minimum of free energy. Solutions qualitatively consistent with the experiment are obtained under the assumption that the mechanism for the formation of self-consistent profiles is directly related to equilibrium diamagnetic currents. The dynamics of the system and specific transport phenomena, such as energy and particle pinching and a decrease in the local density under auxiliary electron cyclotron resonance heating (density pump-out), are analyzed in the vicinity of an equilibrium state characterized by a stable pressure profile. The scaling for the energy confinement time deduced from the transport model agrees qualitatively with the ITER scaling based on the analysis of experimental data obtained in many tokamaks. The possibility of using generalized Tsallis statistics to analyze pressure profiles is considered.

  9. Integrated fusion simulation with self-consistent core-pedestal coupling

    DOE PAGESBeta

    Meneghini, O.; Snyder, P. B.; Smith, S. P.; Candy, J.; Staebler, G. M.; Belli, E. A.; Lao, L. L.; Park, J. M.; Green, D. L.; Elwasif, W.; et al

    2016-04-20

    In this study, accurate prediction of fusion performance in present and future tokamaks requires taking into account the strong interplay between core transport, pedestal structure, current profile and plasma equilibrium. An integrated modeling workflow capable of calculating the steady-state self- consistent solution to this strongly-coupled problem has been developed. The workflow leverages state-of-the-art components for collisional and turbulent core transport, equilibrium and pedestal stability. Validation against DIII-D discharges shows that the workflow is capable of robustly pre- dicting the kinetic profiles (electron and ion temperature and electron density) from the axis to the separatrix in good agreement with the experiments.more » An example application is presented, showing self-consistent optimization for the fusion performance of the 15 MA D-T ITER baseline scenario as functions of the pedestal density and ion effective charge Zeff.« less

  10. Self-consistent Ornstein-Zernike approximation for molecules with soft cores.

    PubMed

    Høye, J S; Reiner, A

    2006-09-14

    The self-consistent Ornstein-Zernike approximation (SCOZA) is an accurate liquid state theory. So far it has been tied to interactions composed of hard core repulsion and long-range attraction, whereas real molecules have soft core repulsion at short distances. In the present work, this is taken into account through the introduction of an effective hard core with a diameter that depends upon temperature only. It is found that the contribution to the configurational internal energy due to the repulsive reference fluid is of prime importance and must be included in the thermodynamic self-consistency requirement on which SCOZA is based. An approximate but accurate evaluation of this contribution relies on the virial theorem to gauge the amplitude of the pair distribution function close to the molecular surface. Finally, the SCOZA equation is transformed by which the problem is reformulated in terms of the usual SCOZA with fixed hard core reference system and temperature-dependent interaction. PMID:16999537

  11. Self-consistent solution of the Dyson equation for atoms and molecules within a conserving approximation.

    PubMed

    Dahlen, Nils Erik; van Leeuwen, Robert

    2005-04-22

    We have calculated the self-consistent Green's function for a number of atoms and diatomic molecules. This Green's function is obtained from a conserving self-energy approximation, which implies that the observables calculated from the Green's functions agree with the macroscopic conservation laws for particle number, momentum, and energy. As a further consequence, the kinetic and potential energies agree with the virial theorem, and the many possible methods for calculating the total energy all give the same result. In these calculations we use the finite temperature formalism and calculate the Green's function on the imaginary time axis. This allows for a simple extension to nonequilibrium systems. We have compared the energies from self-consistent Green's functions to those of nonselfconsistent schemes and also calculated ionization potentials from the Green's functions by using the extended Koopmans' theorem. PMID:15945667

  12. Self-consistent field theory based molecular dynamics with linear system-size scaling

    SciTech Connect

    Richters, Dorothee; Kühne, Thomas D.

    2014-04-07

    We present an improved field-theoretic approach to the grand-canonical potential suitable for linear scaling molecular dynamics simulations using forces from self-consistent electronic structure calculations. It is based on an exact decomposition of the grand canonical potential for independent fermions and does neither rely on the ability to localize the orbitals nor that the Hamilton operator is well-conditioned. Hence, this scheme enables highly accurate all-electron linear scaling calculations even for metallic systems. The inherent energy drift of Born-Oppenheimer molecular dynamics simulations, arising from an incomplete convergence of the self-consistent field cycle, is circumvented by means of a properly modified Langevin equation. The predictive power of the present approach is illustrated using the example of liquid methane under extreme conditions.

  13. Self-consistent models for Coulomb heated X-ray pulsar atmospheres

    NASA Technical Reports Server (NTRS)

    Harding, A.; Meszaros, S. P.; Kirk, J.; Galloway, D.

    1983-01-01

    Calculations of accreting magnetized neutron star atmospheres heated by the gradual deceleration of protons via Coulomb collisions are presented. Self consistent determinations of the temperature and density structure for different accretion rates are made by assuming hydrostatic equilibrium and energy balance, coupled with radiative transfer. The full radiative transfer in two polarizations, using magnetic cross sections but with cyclotron resonance effects treated approximately, is carried out in the inhomogeneous atmospheres.

  14. Steady-state self-consistent analysis of diode-laser arrays

    SciTech Connect

    Kumar, T.

    1987-04-06

    A diode-laser array model, solving the two-dimensional Laplace equation in the passive layer and the carrier diffusion equation in the active layer in a self-consistent manner, is presented. The model is general and does not make use of the usual one-dimensional approximation to the current spreading problem. The beam-propagation method used in solving the wave equation eases the inclusion of any axial variations.

  15. Self-consistent modeling of the electrohydrodynamics of a conductive meniscus

    SciTech Connect

    Wright, G.S.; Krein, P.T.; Chato, J.C.

    1995-07-01

    A complete self-consistent model has been developed for the motion of a conductive liquid surface on a capillary orifice, under an applied electric field. The model uses a quasi-one-dimensional hydrodynamic formulation and a two-dimensional axisymmetric boundary element solution for electric field. The model permits simulation of meniscus behavior with time-varying electric excitation or pressure. Results for resonance behavior and drop emission are presented, with experimental results for comparison.

  16. Accuracy of the Generalized Self-Consistent Method in Modelling the Elastic Behaviour of Periodic Composites

    NASA Technical Reports Server (NTRS)

    Walker, Kevin P.; Freed, Alan D.; Jordan, Eric H.

    1993-01-01

    Local stress and strain fields in the unit cell of an infinite, two-dimensional, periodic fibrous lattice have been determined by an integral equation approach. The effect of the fibres is assimilated to an infinite two-dimensional array of fictitious body forces in the matrix constituent phase of the unit cell. By subtracting a volume averaged strain polarization term from the integral equation we effectively embed a finite number of unit cells in a homogenized medium in which the overall stress and strain correspond to the volume averaged stress and strain of the constrained unit cell. This paper demonstrates that the zeroth term in the governing integral equation expansion, which embeds one unit cell in the homogenized medium, corresponds to the generalized self-consistent approximation. By comparing the zeroth term approximation with higher order approximations to the integral equation summation, both the accuracy of the generalized self-consistent composite model and the rate of convergence of the integral summation can be assessed. Two example composites are studied. For a tungsten/copper elastic fibrous composite the generalized self-consistent model is shown to provide accurate, effective, elastic moduli and local field representations. The local elastic transverse stress field within the representative volume element of the generalized self-consistent method is shown to be in error by much larger amounts for a composite with periodically distributed voids, but homogenization leads to a cancelling of errors, and the effective transverse Young's modulus of the voided composite is shown to be in error by only 23% at a void volume fraction of 75%.

  17. Self-consistent Symmetries in the Proton-Neutron Hartree-Fock-Bogoliubov Approach

    SciTech Connect

    Rohozinski, S.; Dobaczewski, J.; Nazarewicz, Witold

    2010-01-01

    Symmetry properties of densities and mean fields appearing in the nuclear Density Functional Theory with pairing are studied. We consider energy functionals that depend only on local densities and their derivatives. The most important self-consistent symmetries ae discussed: spherical, axial, space-inverson, and mirror symmetries. In each case, the consequences of breaking or conserving the time-reversal and/or proton-neutron symmetries are discussed and summarized in a tabulated form, useful in practical applications.

  18. Magnetic field generation from Self-Consistent collective neutrino-plasma interactions

    SciTech Connect

    Brizard, A.J.; Murayama H.; Wurtele, J.S.

    1999-11-24

    A new Lagrangian formalism for self-consistent collective neutrino-plasma interactions is presented in which each neutrino species is described as a classical ideal fluid. The neutrino-plasma fluid equations are derived from a covariant relativistic variational principle in which finite-temperature effects are retained. This new formalism is then used to investigate the generation of magnetic fields and the production of magnetic helicity as a result of collective neutrino-plasma interactions.

  19. Convergence Properties of the Harris Density Functional and the Self-Consistent Atom Fragment Approximation

    SciTech Connect

    Averill, Frank; Painter, Gayle S

    2006-01-01

    Describing materials properties and behavior over increasing scales of dimension and complexity requires an optimal balance of completeness and accuracy in solving the local density equations. In this study, the convergence properties of a set of schemes that aim to achieve increasing accuracy are systematically examined according to the hierarchical approximations upon which they are based. Specifically, the Harris density functional (HDF) and related schemes that express the total energy in terms of atomic densities and limited self-consistency are compared within a single consistent framework. Convergence of the HDF energy relative to input density is first tested by carrying out calculations within the non-self-consistent atom fragment and self-consistent atom fragment (SCAF) approximations and then by supplementing the SCAF density by increasing numbers of partial waves about each atomic site using the self-consistent partial wave (SCPW) method. The construct of the SCPW method, that solves the local density equations with controlled precision according to the number of partial waves in the site density expansions, enables this study. The rapid convergence of structural properties with an increasing number of partial waves on each site, sometimes even with only L=0 partial waves, provides additional justification for HDF-based tight-binding and molecular dynamics methods where the interatomic potentials are obtained from the superposition of atomic-like densities. The convergence of ground state structural properties is demonstrated by application to the set of molecules: carbon monoxide, water, orthosilicic acid (H{sub 4}SiO{sub 4}) , formamide (HCONH{sub 2}) , iron pentacarbonyl [Fe(CO){sub 5}] , and dimanganese decacarbonyl [Mn{sub 2}(CO){sub 10}] .

  20. Two-way self-consistent coupling of HEIDI in SWMF

    NASA Astrophysics Data System (ADS)

    Ilie, R.; Liemohn, M. W.; Toth, G.

    2013-12-01

    In this study we present results from the two-way coupling between the kinetic Hot Electron and Ion Drift Integrator (HEIDI) model and the Space Weather Modeling Framework (SWMF). HEIDI solves the time dependent, gyration and bounced averaged kinetic equation for the phase space density of different ring current species and computes full pitch angle distributions for all local times and radial distances. This model was generalized to accommodate an arbitrary magnetic field and, through the coupling with SWMF, it obtains the magnetic field description along with the plasma distribution at the model boundaries from the Block Adaptive Tree Solar Wind Roe Upwind Scheme (BATS-R-US) magnetohydrodynamics (MHD) model within the SWMF. Electric field self-consistency is assured by the passing of convection potentials from the Ridley Ionosphere Model (RIM) within SWMF. Our study tests the various levels of coupling between the 3 models, highlighting the roles that the magnetic field, plasma sheet conditions and the cross polar cap potential play in the formation and evolution of the ring current. The results of the self-consistent coupling between HEIDI, BATSRUS and RIM during disturbed conditions emphasize the importance of a kinetic self-consistent approach to the description of the geospace.

  1. Self-Consistent Ring Current/Electromagnetic Ion Cyclotron Waves Modeling

    NASA Technical Reports Server (NTRS)

    Khazanov, G. V.; Gamayunov, K. V.; Gallagher, D. L.

    2006-01-01

    The self-consistent treatment of the RC ion dynamics and EMIC waves, which are thought to exert important influences on the ion dynamical evolution, is an important missing element in our understanding of the storm-and recovery-time ring current evolution. For example, the EMlC waves cause the RC decay on a time scale of about one hour or less during the main phase of storms. The oblique EMIC waves damp due to Landau resonance with the thermal plasmaspheric electrons, and subsequent transport of the dissipating wave energy into the ionosphere below causes an ionosphere temperature enhancement. Under certain conditions, relativistic electrons, with energies 21 MeV, can be removed from the outer radiation belt by EMIC wave scattering during a magnetic storm. That is why the modeling of EMIC waves is critical and timely issue in magnetospheric physics. This study will generalize the self-consistent theoretical description of RC ions and EMIC waves that has been developed by Khazanov et al. [2002, 2003] and include the heavy ions and propagation effects of EMIC waves in the global dynamic of self-consistent RC - EMIC waves coupling. The results of our newly developed model that will be presented at the meeting, focusing mainly on the dynamic of EMIC waves and comparison of these results with the previous global RC modeling studies devoted to EMIC waves formation. We also discuss RC ion precipitations and wave induced thermal electron fluxes into the ionosphere.

  2. Self-consistent Equilibrium Model of Low-aspect-ratio Toroidal Plasma with Energetic Beam Ions

    SciTech Connect

    E.V. Belova; N.N. Gorelenkov; C.Z. Cheng

    2003-04-09

    A theoretical model is developed which allows the self-consistent inclusion of the effects of energetic beam ions in equilibrium calculations of low-aspect-ratio toroidal devices. A two-component plasma is considered, where the energetic ions are treated using a kinetic Vlasov description, while a one-fluid magnetohydrodynamic description is used to represent the thermal plasma. The model allows for an anisotropic distribution function and a large Larmor radius of the beam ions. Numerical results are obtained for neutral-beam-heated plasmas in the National Spherical Torus Experiment (NSTX). Self-consistent equilibria with an anisotropic fast-ion distribution have been calculated for NSTX. It is shown for typical experimental parameters that the contribution of the energetic neutral-beam ions to the total current can be comparable to that of the background plasma, and that the kinetic modifications of the equilibrium can be significant. The range of validity of the finite-Larmor-radius expansion and of the reduced kinetic descriptions for the beam ions in NSTX is discussed. The calculated kinetic equilibria can be used for self-consistent numerical studies of beam-ion-driven instabilities in NSTX.

  3. Two-particle irreducible effective actions versus resummation: Analytic properties and self-consistency

    NASA Astrophysics Data System (ADS)

    Brown, Michael; Whittingham, Ian

    2015-11-01

    Approximations based on two-particle irreducible (2PI) effective actions (also known as Φ-derivable, Cornwall-Jackiw-Tomboulis or Luttinger-Ward functionals depending on context) have been widely used in condensed matter and non-equilibrium quantum/statistical field theory because this formalism gives a robust, self-consistent, non-perturbative and systematically improvable approach which avoids problems with secular time evolution. The strengths of 2PI approximations are often described in terms of a selective resummation of Feynman diagrams to infinite order. However, the Feynman diagram series is asymptotic and summation is at best a dangerous procedure. Here we show that, at least in the context of a toy model where exact results are available, the true strength of 2PI approximations derives from their self-consistency rather than any resummation. This self-consistency allows truncated 2PI approximations to capture the branch points of physical amplitudes where adjustments of coupling constants can trigger an instability of the vacuum. This, in effect, turns Dyson's argument for the failure of perturbation theory on its head. As a result we find that 2PI approximations perform better than Padé approximation and are competitive with Borel-Padé resummation. Finally, we introduce a hybrid 2PI-Padé method.

  4. A second order self-consistent IMEX method for radiation hydrodynamics

    SciTech Connect

    Kadioglu, Samet Y.; Knoll, Dana A.; Lowrie, Robert B.; Rauenzahn, Rick M.

    2010-11-01

    We present a second order self-consistent implicit/explicit (methods that use the combination of implicit and explicit discretizations are often referred to as IMEX (implicit/explicit) methods ) time integration technique for solving radiation hydrodynamics problems. The operators of the radiation hydrodynamics are splitted as such that the hydrodynamics equations are solved explicitly making use of the capability of well-understood explicit schemes. On the other hand, the radiation diffusion part is solved implicitly. The idea of the self-consistent IMEX method is to hybridize the implicit and explicit time discretizations in a nonlinearly consistent way to achieve second order time convergent calculations. In our self-consistent IMEX method, we solve the hydrodynamics equations inside the implicit block as part of the nonlinear function evaluation making use of the Jacobian-free Newton Krylov (JFNK) method . This is done to avoid order reductions in time convergence due to the operator splitting. We present results from several test calculations in order to validate the numerical order of our scheme. For each test, we have established second order time convergence.

  5. Calculations of molecular ionization energies using a self-consistent-charge Hartree-Fock-Slater method

    NASA Technical Reports Server (NTRS)

    Rosen, A.; Ellis, D. E.; Adachi, H.; Averill, F. W.

    1976-01-01

    A numerical-variational method for performing self-consistent molecular calculations in the Hartree-Fock-Slater (HFS) model is presented. Molecular wavefunctions are expanded in terms of basis sets constructed from numerical HFS solutions of selected one-center atomlike problems. Binding energies and wavefunctions for the molecules are generated using a discrete variational method for a given molecular potential. In the self-consistent-charge (SCC) approximation to the complete self-consistent-field (SCF) method, results of a Mulliken population analysis of the molecular eigenfunctions are used in each iteration to produce 'atomic' occupation numbers. The simplest SCC potential is then obtained from overlapping spherical atomlike charge distributions. Molecular ionization energies are calculated using the transition-state procedure; results are given for CO, H2O, H2S, AlCl, InCl, and the Ni5O surface complex. Agreement between experimental and theoretical ionization energies for the free-molecule valence levels is generally within 1 eV. The simple SCC procedure gives a reasonably good approximation to the molecular potential, as shown by comparison with experiment, and with complete SCF calculations for CO, H2O, and H2S.

  6. Recent Progress in GW-based Methods for Excited-State Calculations of Reduced Dimensional Systems

    NASA Astrophysics Data System (ADS)

    da Jornada, Felipe H.

    2015-03-01

    Ab initio calculations of excited-state phenomena within the GW and GW-Bethe-Salpeter equation (GW-BSE) approaches allow one to accurately study the electronic and optical properties of various materials, including systems with reduced dimensionality. However, several challenges arise when dealing with complicated nanostructures where the electronic screening is strongly spatially and directionally dependent. In this talk, we discuss some recent developments to address these issues. First, we turn to the slow convergence of quasiparticle energies and exciton binding energies with respect to k-point sampling. This is very effectively dealt with using a new hybrid sampling scheme, which results in savings of several orders of magnitude in computation time. A new ab initio method is also developed to incorporate substrate screening into GW and GW-BSE calculations. These two methods have been applied to mono- and few-layer MoSe2, and yielded strong environmental dependent behaviors in good agreement with experiment. Other issues that arise in confined systems and materials with reduced dimensionality, such as the effect of the Tamm-Dancoff approximation to GW-BSE, and the calculation of non-radiative exciton lifetime, are also addressed. These developments have been efficiently implemented and successfully applied to real systems in an ab initio framework using the BerkeleyGW package. I would like to acknowledge collaborations with Diana Y. Qiu, Steven G. Louie, Meiyue Shao, Chao Yang, and the experimental groups of M. Crommie and F. Wang. This work was supported by Department of Energy under Contract No. DE-AC02-05CH11231 and by National Science Foundation under Grant No. DMR10-1006184.

  7. Self-consistent calculations of optical properties of type I and type II quantum heterostructures

    NASA Astrophysics Data System (ADS)

    Shuvayev, Vladimir A.

    In this Thesis the self-consistent computational methods are applied to the study of the optical properties of semiconductor nanostructures with one- and two-dimensional quantum confinements. At first, the self-consistent Schrodinger-Poisson system of equations is applied to the cylindrical core-shell structure with type II band alignment without direct Coulomb interaction between carriers. The electron and hole states and confining potential are obtained from a numerical solution of this system. The photoluminescence kinetics is theoretically analyzed, with the nanostructure size dispersion taken into account. The results are applied to the radiative recombination in the system of ZnTe/ZnSe stacked quantum dots. A good agreement with both continuous wave and time-resolved experimental observations is found. It is shown that size distribution results in the photoluminescence decay that has essentially non-exponential behavior even at the tail of the decay where the carrier lifetime is almost the same due to slowly changing overlap of the electron and hole wavefunctions. Also, a model situation applicable to colloidal core-shell nanowires is investigated and discussed. With respect to the excitons in type I quantum wells, a new computationally efficient and flexible approach of calculating the characteristics of excitons, based on a self-consistent variational treatment of the electron-hole Coulomb interaction, is developed. In this approach, a system of self-consistent equations describing the motion of an electron-hole pair is derived. The motion in the growth direction of the quantum well is separated from the in-plane motion, but each of them occurs in modified potentials found self-consistently. This approach is applied to a shallow quantum well with the delta-potential profile, for which analytical expressions for the exciton binding energy and the ground state eigenfunctions are obtained, and to the quantum well with the square potential profile with several

  8. Beyond the GW approximation: a second-order screened exchange correction

    NASA Astrophysics Data System (ADS)

    Rinke, Patrick; Caruso, Fabio; Ren, Xinguo; Scheffler, Matthias; Marom, Noa

    2013-03-01

    Despite the success of the GW method in describing the photoemission spectra of solids, molecules and clusters, challenges remain. For aromatic molecules for example absolute as well as relative positions of ionisation energies and affinities are not well reproduced in perturbative G0W0 schemes with different starting points as well as in self-consistent GW, sometimes even giving the wrong orbital order. Motivated by renormalized second-order perturbation theory for the ground-state energy, we propose a second-order screened exchange correction (SOSEX) to the GW self-energy. This correction follows the spirit of the SOSEX correction to the random-phase approximation for the electron correlation energy and reduces the self-correlation error. The performance of the GW +SOSEX scheme has been benchmarked for a set of molecular systems, including the G2 set, commonly used acceptor molecules, benzene and the azabenzene molecules. We find that the SOSEX correction improves the description of the spectral properties including the orbital order with respect to the different GW schemes, highlighting the importance of reducing the self-correlation error.

  9. GW in the Gaussian and Plane Waves Scheme with Application to Linear Acenes.

    PubMed

    Wilhelm, Jan; Del Ben, Mauro; Hutter, Jürg

    2016-08-01

    We present an implementation of G0W0 and eigenvalue-self-consistent GW (evGW) in the Gaussian and plane waves scheme for molecules. We calculate the correlation self-energy for imaginary frequencies employing the resolution of the identity. The correlation self-energy for real frequencies is then evaluated by analytic continuation. This technique allows an efficient parallel implementation and application to systems with several hundreds of atoms. Various benchmark calculations are presented. In particular, the convergence with respect to the most important numerical parameters is assessed for the benzene molecule. Comparisons with respect to other G0W0 implementations are reported for a set of molecules, while the performance of the method has been measured for water clusters containing up to 480 atoms in a cc-TZVP basis. Additionally, G0W0 has been applied for studying the influence of the ligands on the gap of small CdSe nanoparticles. evGW has been employed to calculate the HOMO-LUMO gaps of linear acenes, linear chains formed of connected benzene rings. Distinct differences between the closed and the open-shell (broken-symmetry) evGW HOMO-LUMO gaps for long acenes are found. In future experiments, a comparison of measured HOMO-LUMO gaps and our calculated evGW values may be helpful to determine the electronic ground state of long acenes. PMID:27348184

  10. Some exploitations of the self-consistent QRPA approach with the Gogny force

    SciTech Connect

    Peru, S.; Martini, M.; Dupuis, M.

    2012-10-20

    Fully consistent axially-symmetric-deformed quasiparticle random phase approximation calculations have been performed with the D1S Gogny force. Giant resonances in exotic nuclei as well as in deformed Mg and Si isotopes have been studied. Dipole responses have been calculated in Ne isotopes and N=16 isotones to study the existence of soft dipole modes in exotic nuclei. The same formalism has been used to describe multipole responses up to octupole in the deformed and heavy nucleus {sup 238}U. Low energy spectroscopy of nickel isotopes has been studied, revealing 0{sup +} states which display a particular structure.

  11. Self-Consistent Green Function Calculations for Isospin Asymmetric Nuclear Matter

    NASA Astrophysics Data System (ADS)

    Mansour, H.; Gad, K.; Hassaneen, K. S. A.

    2010-04-01

    The one-body potentials for protons and neutrons are obtained from the self-consistent Green-function calculations of asymmetric nuclear matter, in particular their dependence on the degree of proton/neutron asymmetry. Results of the binding energy per nucleon as a function of the density and asymmetry parameter are presented for the self-consistent Green function approach using the CD-Bonn potential. For the sake of comparison, the same calculations are performed using the Brueckner-Hartree-Fock approximation. The contribution of the hole-hole terms leads to a repulsive contribution to the energy per nucleon which increases with the nuclear density. The incompressibility for asymmetric nuclear matter has been also investigated in the framework of the self-consistent Green-function approach using the CD-Bonn potential. The behavior of the incompressibility is studied for different values of the nuclear density and the neutron excess parameter. The nuclear symmetry potentia l at fixed nuclear density is also calculated and its value decreases with increasing the nucleon energy. In particular, the nuclear symmetry potential at saturation density changes from positive to negative values at nucleon kinetic energy of about 200 MeV. For the sake of comparison, the same calculations are performed using the Brueckner-Hartree-Fock approximation. The proton/neutron effective mass splitting in neutron-rich matter has been studied. The predicted isospin splitting of the proton/neutron effective mass splitting in neutron-rich matter is such that m(*_n) ≥ m(*_p) .

  12. Chiral and magnetic rotation in atomic nuclei studied within self-consistent mean-field methods

    NASA Astrophysics Data System (ADS)

    Olbratowski, P.

    2004-07-01

    Currently, one application of the mean-field methods in nuclear physics is the investigation of exotic nuclear symmetries. This is related, in particular, to the study of nuclear rotation about an axis tilted with respect to the principal axes of the mass distribution in the Tilted-Axis Cranking (TAC) model. The present work presents one of the first TAC calculations performed within fully self-consistent methods. The Hartree-Fock method with the Skyrme effective two-body interaction has been used. A computer code has been developed that allows for the breaking of all spatial symmetries of the solution. As a first application, calculations for the magnetic bands in 142Gd and for the chiral bands in 130Cs, 132La, 134Pr, and 136Pm have been carried out. The appearance of those bands is due to a new mechanism of breaking the spherical symmetry and to the spontaneous breaking of the chiral symmetry, respectively. The self-consistent solutions for 142Gd confirm the important role of the shears mechanism in generating the total angular momentum. However, the agreement with experimental data is not satisfactory, probably due to the lack of the pairing correlations in the calculations or to the possibly overestimated deformation. The results obtained for 132La constitute the first fully self-consistent proof that the nuclear rotation can attain a chiral character. It has been shown that the chiral rotation can only exist above a certain critical angular frequency. It has also been checked that the terms of the Skyrme mean field odd under the time reversal have no qualitative influence on the results.

  13. Self-consistent generalized Langevin-equation theory for liquids of nonspherically interacting particles.

    PubMed

    Elizondo-Aguilera, L F; Zubieta Rico, P F; Ruiz-Estrada, H; Alarcón-Waess, O

    2014-11-01

    A self-consistent generalized Langevin-equation theory is proposed to describe the self- and collective dynamics of a liquid of linear Brownian particles. The equations of motion for the spherical harmonics projections of the collective and self-intermediate-scattering functions, F_{lm,lm}(k,t) and F_{lm,lm}^{S}(k,t), are derived as a contraction of the description involving the stochastic equations of the corresponding tensorial one-particle density n_{lm}(k,t) and the translational (α=T) and rotational (α=R) current densities j_{lm}^{α}(k,t). Similar to the spherical case, these dynamic equations require as an external input the equilibrium structural properties of the system contained in the projections of the static structure factor, denoted by S_{lm,lm}(k). Complementing these exact equations with simple (Vineyard-like) approximate relations for the collective and the self-memory functions we propose a closed self-consistent set of equations for the dynamic properties involved. In the long-time asymptotic limit, these equations become the so-called bifurcation equations, whose solutions (the nonergodicity parameters) can be written, extending the spherical case, in terms of one translational and one orientational scalar dynamic order parameter, γ_{T} and γ_{R}, which characterize the possible dynamical arrest transitions of the system. As a concrete illustrative application of this theory we determine the dynamic arrest diagram of the dipolar hard-sphere fluid. In qualitative agreement with mode coupling theory, the present self-consistent equations also predict three different regions in the state space spanned by the macroscopic control parameters η (volume fraction) and T* (scaled temperature): a region of fully ergodic states, a region of mixed states, in which the translational degrees of freedom become arrested while the orientational degrees of freedom remain ergodic, and a region of fully nonergodic states. PMID:25493790

  14. Quasiparticle spectra from a nonempirical optimally tuned range-separated hybrid density functional

    SciTech Connect

    Refaely-Abramson, Sivan; Sharifzadeh, Sahar; Govind, Niranjan; Autschbach, Jochen; Neaton, Jeffrey B.; Baer, Roi; Kronik, Leeor

    2012-11-28

    We present a method for obtaining quasiparticle excitation energies from a DFT-based calculation, but with accuracy that is comparable to that of many-body perturbation theory within the GW approximation. The approach uses a range-separated hybrid density functional, with asymptotically exact and short-range fractional Fock exchange. The functional contains two parameters - the range separation and the short-range Fock fraction. Both are determined non-empirically, per system, based on satisfaction of exact physical constraints for the ionization potential and many-electron self-interaction, respectively. The accuracy of the method is demonstrated on the important benchmark molecule, 3,4,9,10-perylene-tetracarboxylic-dianydride (PTCDA), where it is shown to be the only non-empirical DFT-based method comparable to GW calculations. For any finite system, we envision that the approach could be useful directly as an inexpensive alternative to GW that offers good accuracy for both frontier and non-frontier quasiparticle excitation energies, opening the door to the studyof presently out of reach large-scale systems.

  15. Solvent effects in time-dependent self-consistent field methods. I. Optical response calculations.

    PubMed

    Bjorgaard, J A; Kuzmenko, V; Velizhanin, K A; Tretiak, S

    2015-01-28

    We implement and examine three excited state solvent models in time-dependent self-consistent field methods using a consistent formalism which unambiguously shows their relationship. These are the linear response, state specific, and vertical excitation solvent models. Their effects on energies calculated with the equivalent of COSMO/CIS/AM1 are given for a set of test molecules with varying excited state charge transfer character. The resulting solvent effects are explained qualitatively using a dipole approximation. It is shown that the fundamental differences between these solvent models are reflected by the character of the calculated excitations. PMID:25637965

  16. The concept of coupling impedance in the self-consistent plasma wake field excitation

    NASA Astrophysics Data System (ADS)

    Fedele, R.; Akhter, T.; De Nicola, S.; Migliorati, M.; Marocchino, A.; Massimo, F.; Palumbo, L.

    2016-09-01

    Within the framework of the Vlasov-Maxwell system of equations, we describe the self-consistent interaction of a relativistic charged-particle beam with the surroundings while propagating through a plasma-based acceleration device. This is done in terms of the concept of coupling (longitudinal) impedance in full analogy with the conventional accelerators. It is shown that also here the coupling impedance is a very useful tool for the Nyquist-type stability analysis. Examples of specific physical situations are finally illustrated.

  17. Self-consistent calculation of dephasing in quantum cascade structures within a density matrix method

    NASA Astrophysics Data System (ADS)

    Freeman, Will

    2016-05-01

    Dephasing in terahertz quantum cascade structures is studied within a density matrix formalism. We self-consistently calculate the pure dephasing time from the intrasubband interactions within the upper and lower lasing states. Interface roughness and ionized impurity scattering interactions are included in the calculation. Dephasing times are shown to be consistent with measured spontaneous emission spectra, and the lattice temperature dependence of the device output power is consistent with experiment. The importance of including multiple optical transitions when a lower miniband continuum is present and the resulting multi-longitudinal modes within the waveguide resonant cavity are also shown.

  18. Self-consistent relativistic random-phase approximation with vacuum polarization

    SciTech Connect

    Haga, A.; Toki, H.; Tamenaga, S.; Horikawa, Y.; Yadav, H.L.

    2005-09-01

    We present a theoretical formulation for the description of nuclear excitations within the framework of a relativistic random-phase approximation whereby the vacuum polarization arising from nucleon-antinucleon fields is duly accounted for. The vacuum contribution to the Lagrangian is explicitly described as extra new terms of interacting mesons by means of the derivative expansion of the effective action. It is shown that the self-consistent calculation yields zero eigenvalue for the spurious isoscalar-dipole state and also conserves the vector-current density.

  19. A self-consistent quasistatic equilibrium for non-neutral diamagnetic electron vortices

    SciTech Connect

    Church, B.W.; Sudan, R.N.

    1995-06-01

    A self-consistent quasistatic equilibrium for a non-neutral cylindrical electron vortex has been found using the two-dimensional relativistic electron fluid equations. While other work on electron vortices considered a regime where the vortex radius is much smaller than the collisionless skin depth {lambda}={ital c}/{omega}{sub {ital p}}, this equilibrium is valid for large-radius, diamagnetic vortices and predicts a maximum radius of 2{sup 3/2}{lambda} for a highly relativistic electron vortex. The vortex model shows good agreement with observations of diamagnetic electron vortices in two-dimensional electromagnetic particle-in-cell simulations of magnetically insulated transmission lines.

  20. Self-consistent Hartree-Fock approach for interacting bosons in optical lattices

    NASA Astrophysics Data System (ADS)

    Lü, Qin-Qin; Patton, Kelly R.; Sheehy, Daniel E.

    2014-12-01

    A theoretical study of interacting bosons in a periodic optical lattice is presented. Instead of the commonly used tight-binding approach (applicable near the Mott-insulating regime of the phase diagram), the present work starts from the exact single-particle states of bosons in a cubic optical lattice, satisfying the Mathieu equation, an approach that can be particularly useful at large boson fillings. The effects of short-range interactions are incorporated using a self-consistent Hartree-Fock approximation, and predictions for experimental observables such as the superfluid transition temperature, condensate fraction, and boson momentum distribution are presented.

  1. Interstellar turbulence model : A self-consistent coupling of plasma and neutral fluids

    SciTech Connect

    Shaikh, Dastgeer; Zank, Gary P.; Pogorelov, Nikolai

    2006-09-26

    We present results of a preliminary investigation of interstellar turbulence based on a self-consistent two-dimensional fluid simulation model. Our model describes a partially ionized magnetofluid interstellar medium (ISM) that couples a neutral hydrogen fluid to a plasma through charge exchange interactions and assumes that the ISM turbulent correlation scales are much bigger than the shock characteristic length-scales, but smaller than the charge exchange mean free path length-scales. The shocks have no influence on the ISM turbulent fluctuations. We find that nonlinear interactions in coupled plasma-neutral ISM turbulence are influenced substantially by charge exchange processes.

  2. Homogenization of Periodic Masonry Using Self-Consistent Scheme and Finite Element Method

    NASA Astrophysics Data System (ADS)

    Kumar, Nitin; Lambadi, Harish; Pandey, Manoj; Rajagopal, Amirtham

    2016-01-01

    Masonry is a heterogeneous anisotropic continuum, made up of the brick and mortar arranged in a periodic manner. Obtaining the effective elastic stiffness of the masonry structures has been a challenging task. In this study, the homogenization theory for periodic media is implemented in a very generic manner to derive the anisotropic global behavior of the masonry, through rigorous application of the homogenization theory in one step and through a full three-dimensional behavior. We have considered the periodic Eshelby self-consistent method and the finite element method. Two representative unit cells that represent the microstructure of the masonry wall exactly are considered for calibration and numerical application of the theory.

  3. Self-consistent, three-dimensional equilibrium effects on tokamak magnetic field ripple

    SciTech Connect

    Johnson, J.L.; Reiman, A.H.

    1987-10-01

    Self-consistent equilibrium effects on tokamak magnetic field ripple have been calculated using a three-dimensional equilibrium code. The effects are found to be large enough that they should be included in tokamak ignition experiment designs. Even the modification of the well depth associated with the flow of force-free plasma current along rippled field lines is substantial. An analysis of the results separates the contribution of the Shafranov shift to the ripple modification from the contributions of other finite-pressure effects. 5 refs., 10 figs., 1 tab.

  4. Study of transient self-consistent beam dynamics in RF linacs using a particle tracing code

    NASA Astrophysics Data System (ADS)

    Mytrochenko, V. V.; Opanasenko, A.

    2006-03-01

    The paper describes a simulation technique for study of unsteady self-consistent dynamics of charged particles in RF linacs that consist of cavities and travelling wave sections. The approach proposed is based on unsteady theories of the excitation of cavities and waveguides by a beam of charged particles and RF feeders. The theory of waveguide excitation is extended to the case of spatially inhomogeneous travelling-wave structures. The SUPERFISH code is used to evaluate the characteristics of the axisymmetric travelling-wave sections. The PARMELA code is applied for simulation of particle motion and to obtain data required for solving the equations for excitation of the RF structures by the beam.

  5. Self-consistent model of spin accumulation magnetoresistance in ferromagnet/insulator/semiconductor tunnel junctions

    NASA Astrophysics Data System (ADS)

    Appelbaum, Ian; Tinkey, Holly N.; Li, Pengke

    2014-12-01

    Spin accumulation in a paramagnetic semiconductor due to voltage-biased current tunneling from a polarized ferromagnet is experimentally manifest as a small additional spin-dependent resistance. We describe a rigorous model incorporating the necessary self-consistency between electrochemical potential splitting, spin-dependent injection current, and applied voltage that can be used to simulate this so-called "3 T " signal as a function of temperature, doping, ferromagnet bulk spin polarization, tunnel barrier features and conduction nonlinearity, and junction voltage bias.

  6. Self-Consistent Field Approach for Cross-Linked Copolymer Materials

    NASA Astrophysics Data System (ADS)

    Schmid, Friederike

    2013-07-01

    A generalized self-consistent field approach for polymer networks with a fixed topology is developed. It is shown that the theory reproduces the localization of cross-links, which is characteristic for gels. The theory is then used to study the order-disorder transition in regular networks of end-linked diblock copolymers. Compared to diblock copolymer melts, the transition is shifted towards lower values of the incompatibility parameter χ (the Flory- Huggins parameter). Moreover, the transition becomes strongly first order already at the mean-field level. If stress is applied, the transition is further shifted and finally vanishes in a critical point.

  7. Self-consistent description of coexistence phenomena in medium mass nuclei

    SciTech Connect

    Petrovici, A.; Schmid, K. W.; Faessler, Amand; Andrei, O.

    2010-11-24

    Shape coexistence and mixing, isospin mixing, the competition between neutron-proton and like-nucleon pairing correlations have been identified as the main characteristic features of nuclei near the N = Z line in the A{approx_equal}70 mass region. The self-consistent treatment of exotic phenomena dominated by their interplay represents a challenge for the nuclear many-body models. The realistic description of tiny effects in this mass region aiming to test the fundamental interactions and symmetries as well as the required theoretical predictions concerning the nuclear properties relevant for astrophysical scenarios are still open problems of the low-energy nuclear physics today.

  8. Solvent effects in time-dependent self-consistent field methods. I. Optical response calculations

    SciTech Connect

    Bjorgaard, J. A.; Kuzmenko, V.; Velizhanin, K. A.; Tretiak, S.

    2015-01-22

    In this study, we implement and examine three excited state solvent models in time-dependent self-consistent field methods using a consistent formalism which unambiguously shows their relationship. These are the linear response, state specific, and vertical excitation solvent models. Their effects on energies calculated with the equivalent of COSMO/CIS/AM1 are given for a set of test molecules with varying excited state charge transfer character. The resulting solvent effects are explained qualitatively using a dipole approximation. It is shown that the fundamental differences between these solvent models are reflected by the character of the calculated excitations.

  9. Self-consistent continuum random-phase approximation calculations of {sup 4}He electromagnetic responses

    SciTech Connect

    De Donno, V.; Co', G.; Anguiano, M.; Lallena, A. M.

    2011-09-15

    We study the electromagnetic responses of {sup 4}He within the framework of the self-consistent continuum random-phase approximation theory. In this approach, the ground-state properties are described by a Hartree-Fock calculation. The single-particle basis constructed in this manner is used in the calculations of the continuum responses of the system. Finite-range interactions are considered in the calculations. We compare our results with photon-absorption cross sections and electron-scattering quasielastic data. From this comparison, and also from the comparison with the results of microscopic calculations, we deduce that our approach describes well the continuum excitation.

  10. A self-consistent approach to the reflection component in 4U 1705-44

    NASA Astrophysics Data System (ADS)

    D'Aì, A.; di Salvo, T.; Iaria, R.; Papitto, A.; Ballantyne, D.

    2010-07-01

    High-resolution spectroscopy has recently revealed in many neutron-star Low-Mass X-ray binaries that the shape of the broad iron line observed in the 6.4-6.97 keV range is consistently well fitted by a relativistically smeared line profile. We show here spectral fitting results using a newly developed self-consistent reflection model on XMM-Newton data of the LMXB 4U 1705-44 during a period when the source was in a bright soft state. This reflection model adopts a blackbody prescription for the shape of the impinging radiation field, that we physically associate with the boundary layer emission.

  11. Accurate vibrational frequencies using the self-consistent-charge density-functional tight-binding method

    NASA Astrophysics Data System (ADS)

    Małolepsza, Edyta; Witek, Henryk A.; Morokuma, Keiji

    2005-09-01

    An optimization technique for enhancing the quality of repulsive two-body potentials of the self-consistent-charge density-functional tight-binding (SCC-DFTB) method is presented and tested. The new, optimized potentials allow for significant improvement of calculated harmonic vibrational frequencies. Mean absolute deviation from experiment computed for a group of 14 hydrocarbons is reduced from 59.0 to 33.2 cm -1 and maximal absolute deviation, from 436.2 to 140.4 cm -1. A drawback of the new family of potentials is a lower quality of reproduced geometrical and energetic parameters.

  12. Properties of hadronic systems according to the non-extensive self-consistent thermodynamics

    SciTech Connect

    Deppman, A.

    2014-11-11

    The non-extensive self-consistent theory describing the thermodynamics of hadronic systems at high temperatures is used to derive some thermodynamical quantities, as pressure, entropy, speed of sound and trace-anomaly. The calculations are free of fitting parameters, and the results are compared to lattice QCD calculations, showing a good agreement between theory and data up to temperatures around 175 MeV. Above this temperature the effects of a singularity in the partition function at T{sub o} = 192 MeV results in a divergent behaviour in respect with the lattice calculation.

  13. Symmetry breaking and self-consistent rotation of magnetic islands in neoclassical viscous regimes

    SciTech Connect

    Lazzaro, E.

    2009-09-15

    Classical or neoclassical tearing modes (NTMs) perturb the ideal axisymmetry of tokamaks. As a consequence of symmetry breaking a neoclassical toroidal viscosity (NTV) appears, that depends on the island amplitude. This work shows that in the low collisionality regimes NTV has a key role in determining self-consistently the magnetic island velocity and at the same time modifies significantly the ion polarization current effects on NTM instability. This finding can provide a better understanding of the mechanism of onset of NTMs, observed experimentally, and improve the concepts for their control or avoidance.

  14. Pressure tensor in the presence of velocity shear: Stationary solutions and self-consistent equilibria

    SciTech Connect

    Cerri, S. S.; Pegoraro, F.; Califano, F.; Jenko, F.

    2014-11-15

    Observations and numerical simulations of laboratory and space plasmas in almost collisionless regimes reveal anisotropic and non-gyrotropic particle distribution functions. We investigate how such states can persist in the presence of a sheared flow. We focus our attention on the pressure tensor equation in a magnetized plasma and derive analytical self-consistent plasma equilibria which exhibit a novel asymmetry with respect to the magnetic field direction. These results are relevant for investigating, within fluid models that retain the full pressure tensor dynamics, plasma configurations where a background shear flow is present.

  15. GW calculations of band offsets at AlN/GaN interfaces.

    NASA Astrophysics Data System (ADS)

    Cociorva, Daniel; Wilkins, John W.

    1998-03-01

    We study the quasiparticle bands at the interface of AlN and GaN using the GW approximation for the electron self-energy. Two structures are studied: the zinc-blende (cubic) interface in the (111) direction and the wurtzite (hexagonal) interface in the (001) direction. Our short superlattices are consistent with 4×4 LDA work(F. Bernardini et al.), Materials Research Society symposia proceedings 449, 923 (1997). and ``1×2'' GW work.(A. Rubio et al.), Phys. Rev. B 49, 1952 (1994). For comparison with recent experiments we are implementing a scalable parallel algorithm based on the reciprocal space formulation(M. S. Hybertsen and S. G. Louie, Phys. Rev. Lett. 55), 1418 (1985). to determine the self energy for thicker superlattices.

  16. Radio-frequency sheaths physics: Experimental characterization on Tore Supra and related self-consistent modeling

    NASA Astrophysics Data System (ADS)

    Jacquot, Jonathan; Milanesio, Daniele; Colas, Laurent; Corre, Yann; Goniche, Marc; Gunn, Jamie; Heuraux, Stéphane; Kubič, Martin

    2014-06-01

    During the 2011 experimental campaign, one of the three ion cyclotron resonance heating (ICRH) antennas in the Tore Supra tokamak was equipped with a new type of Faraday screen (FS). The new design aimed at minimizing the integrated parallel electric field over long field lines as well as increasing the heat exhaust capability of the actively cooled screen. It proved to be inefficient for attenuating the radio-frequency (RF)-sheaths on the screen itself on the contrary to the heat exhaust concept that allowed operation despite higher heat fluxes on the antenna. In parallel, a new approach has been proposed to model self-consistently RF sheaths: the SSWICH (Self-consistent Sheaths and Waves for IC Heating) code. Simulations results from SSWICH coupled with the TOPICA antenna code were able to reproduce the difference between the two FS designs and part of the spatial pattern of heat loads and Langmuir probe floating potential. The poloidal pattern is a reliable result that mainly depends on the electrical design of the antenna while the radial pattern is on the contrary highly sensitive to loosely constrained parameters such as perpendicular conductivity that generates a DC current circulation from the private region inside the antenna limiters to the free scrape off layer outside these limiters. Moreover, the cantilevered bars seem to be the element in the screen design that enhanced the plasma potential.

  17. A modified self-consistent computational model for an electrostrictive graft elastomer

    NASA Astrophysics Data System (ADS)

    Sun, Changjie; Wang, Youqi; Zhang, Congjian; Zhou, Eric; Su, Ji

    2008-04-01

    An electrostrictive graft elastomer, as recently developed by NASA, is a type of electro-active polymer. In this paper, a 2D computational model with a self-consistent boundary is developed. Firstly, three-dimensional deformations, induced by both bending angle and dihedral torsional angle changes, are projected onto a two-dimensional plane. Using both theoretical and numerical analyses, the projected 2D equilibrium bending angle is shown to have the same value as the 3D equilibrium bending angle. The 2D equivalent bending stiffness is derived using a series model based upon the fact that both bending and dihedral torsion produce a configurational change. Equivalent stiffness is justified by polymer chain end-to-end distance characteristics. Secondly, a self-consistent scheme is developed to eliminate the boundary effect. Eight images of the unit cell are created peripherally, with the original unit cell in the center. Thus the boundary can only affect the rotation of the eight images, not the central unit cell. A computational model is employed to determine the electromechanical properties of the electrostrictive graft elastomer. Relations between electric field induced strain and electric field strength are calculated. The effect of molecular scale factors, such as free volume fraction, graft weight percentage and graft orientation, are also discussed.

  18. Self-Consistent Simulations of Accretion-Induced Collapse of White Dwarfs

    NASA Astrophysics Data System (ADS)

    Kleiser, Io; Ott, Christian; Abdikamalov, Ernazar; O'Connor, Evan

    2013-04-01

    Accreting white dwarfs and white dwarf mergers are commonly thought to end in thermonuclear explosions that produce Type Ia supernovae (SNe Ia). However, there is an alternative outcome for these systems that has not been theoretically explored as thoroughly, nor has it been securely identified observationally. Some white dwarfs, rather than exploding, should undergo electron capture and collapse to neutron stars. This accretion-induced collapse (AIC) scenario is expected to be intrinsically rare compared to SNe Ia, and past studies indicate that the associated optical transient would be faint and short-lived, near the detection limits of current surveys. However, until now there have not been self-consistent numerical studies of AIC that examine the explosion dynamics, subsequent evolution, and all resulting observables. We use GR1D, a one-dimensional general-relativistic hydrodynamics code, to follow AIC through collapse, core bounce, explosion, and shock breakout and to present new results on its neutrino signature and nucleosynthetic yields. This study is preliminary to the goal of developing fully self-consistent three-dimensional models that will yield predictions for electromagnetic, neutrino, and gravitational-wave signals form AIC events.

  19. Self-Consistent Simulations of Accretion-Induced Collapse of White Dwarfs

    NASA Astrophysics Data System (ADS)

    Kleiser, Io; Ott, C. D.

    2013-01-01

    Accreting white dwarfs and white dwarf mergers are commonly thought to end in thermonuclear explosions that produce Type Ia supernovae (SNe Ia). However, there is an alternative outcome for these systems that has not been theoretically explored as thoroughly, nor has it been securely identified observationally. Some white dwarfs, rather than exploding, should undergo electron capture and collapse to neutron stars. This accretion-induced collapse (AIC) scenario is expected to be intrinsically rare compared to SNe Ia, and past studies indicate that the associated optical transient would be faint and short-lived, near the detection limits of current surveys. However, until now there have not been self-consistent numerical studies of AIC that examine the explosion dynamics, subsequent evolution, and all resulting observables. We use GR1D, a one-dimensional general-relativistic hydrodynamics code, to follow AIC through collapse, core bounce, explosion, and shock breakout and to present new results on its neutrino signature. This study is preliminary to the goal of developing fully self-consistent three-dimensional models that will yield predictions for electromagnetic, neutrino, and gravitational-wave signals form AIC events.

  20. False vacuum decay by self-consistent bounces in four dimensions

    NASA Astrophysics Data System (ADS)

    Baacke, Jürgen; Kevlishvili, Nina

    2007-02-01

    We compute bounce solutions describing false vacuum decay in a Φ4 model in four dimensions with quantum backreaction. The backreaction of the quantum fluctuations on the bounce profiles is computed in the one-loop and Hartree approximations. This is to be compared with the usual semiclassical approach where one computes the profile from the classical action and determines the one-loop correction from this profile. The computation of the fluctuation determinant is performed using a theorem on functional determinants, in addition we here need the Green’s function of the fluctuation operator in oder to compute the quantum backreaction. As we are able to separate from the determinant and from the Gree n’s function the leading perturbative orders, we can regularize and renormalize analytically, in analogy of standard perturbation theory. The iteration towards self-consistent solutions is found to converge for some range of the parameters. Within this range the corrections to the semiclassical action are at most a few percent, the corrections to the transition rate can amount to several orders of magnitude. The strongest deviations happen for large couplings, as to be expected. The transition rates are reduced for the one-loop backreaction, for the Hartree backreaction they are reduced for α≲0.5 and enhanced for larger values of α. Beyond some limit, there are no self-consistent bounce solutions.

  1. Self-Consistent Conversion of a Viscous Fluid to Particles and Heavy-Ion Physics Applications

    NASA Astrophysics Data System (ADS)

    Wolff, Zack J.

    The most widely used theoretical framework to model the early stages of a heavy-ion collision is viscous hydrodynamics. Comparing hydrodynamic simulations to heavy-ion data inevitably requires the conversion of the fluid to particles. This conversion, typically done in the Cooper-Frye formalism, is ambiguous for viscous fluids. In this thesis work, self-consistent phase space corrections are calculated by solving the linearized Boltzmann equation. These species-dependent solutions are contrasted with those obtained using the ad-hoc ''democratic Grad'' ansatz typically employed in the literature in which coefficients are independent of particle dynamics. Solutions are calculated analytically for a massless gas and numerically for the general case of a hadron resonance gas. For example, it is found that for a gas of massless particles interacting via isotropic, energy-independent 2 → 2 scatterings, the shear viscous corrections variationally prefer a momentum dependence close to p3/2 rather than the quadratic dependence assumed in the Grad ansatz. The self-consistent phase space distributions are then used to calculate transverse momentum spectra and differential flow coefficients, v n(pT), to study the effects on heavy-ion identified particle observables. Using additive quark model cross sections, it is found that proton flow coefficients are higher than those for pions at moderately high pT in Pb + Pb collisions at LHC, especially for the coefficients v 4 and v6.

  2. "Hairy" Nanoparticles in Block Copolymers and Homopolymers: Modeling using Hybrid Self-Consistent Field Theory

    NASA Astrophysics Data System (ADS)

    Ginzburg, Valeriy

    2011-03-01

    Today, dispersed nanoparticles play important role in various applications (toughened plastics, healthcare, personal care, etc.) Mesoscale simulations and theory are important in understanding what governs the morphology of nanoparticles under various conditions. In particular, for nanoparticle/block copolymer mixtures, two popular simulation methods are Self-Consistent Field/Density Functional Theory (SCF-DFT) (Thompson, Ginzburg, Matsen, and Balazs, Science 292, 2469 [2001]), and Hybrid Self-Consistent Field Theory (HSCFT) (Sides et al., Phys Rev Lett 96, 250601 [2006]). The two methods are shown to be very similar in their assumptions and end-results; the choice of the method to be used can depend on the specific problem. Here, we use modified HSCFT to explicitly account for the complicated role of short-chain ligands grafted onto nanoparticles to promote dispersion. In particular, we discuss the phase diagrams of such ``hairy'' nanoparticles in diblock copolymers as function of diblock composition, nanoparticle volume fraction, and ligand length. Depending on the particle size and ligand coverage, particles could segregate into favorable domain, stay close to the interface, or phase-separate from the block copolymer altogether. We also consider the dispersion of ``hairy'' nanoparticles in a homopolymer and analyze the morphologies of particle clusters as function of ligand length. The results could have interesting implications for the design of new nanocomposite materials.

  3. Radio-frequency sheaths physics: Experimental characterization on Tore Supra and related self-consistent modeling

    SciTech Connect

    Jacquot, Jonathan; Colas, Laurent Corre, Yann; Goniche, Marc; Gunn, Jamie; Kubič, Martin; Milanesio, Daniele; Heuraux, Stéphane

    2014-06-15

    During the 2011 experimental campaign, one of the three ion cyclotron resonance heating (ICRH) antennas in the Tore Supra tokamak was equipped with a new type of Faraday screen (FS). The new design aimed at minimizing the integrated parallel electric field over long field lines as well as increasing the heat exhaust capability of the actively cooled screen. It proved to be inefficient for attenuating the radio-frequency (RF)-sheaths on the screen itself on the contrary to the heat exhaust concept that allowed operation despite higher heat fluxes on the antenna. In parallel, a new approach has been proposed to model self-consistently RF sheaths: the SSWICH (Self-consistent Sheaths and Waves for IC Heating) code. Simulations results from SSWICH coupled with the TOPICA antenna code were able to reproduce the difference between the two FS designs and part of the spatial pattern of heat loads and Langmuir probe floating potential. The poloidal pattern is a reliable result that mainly depends on the electrical design of the antenna while the radial pattern is on the contrary highly sensitive to loosely constrained parameters such as perpendicular conductivity that generates a DC current circulation from the private region inside the antenna limiters to the free scrape off layer outside these limiters. Moreover, the cantilevered bars seem to be the element in the screen design that enhanced the plasma potential.

  4. Fractional charge and spin errors in self-consistent Green's function theory

    NASA Astrophysics Data System (ADS)

    Phillips, Jordan J.; Kananenka, Alexei A.; Zgid, Dominika

    2015-05-01

    We examine fractional charge and spin errors in self-consistent Green's function theory within a second-order approximation (GF2). For GF2, it is known that the summation of diagrams resulting from the self-consistent solution of the Dyson equation removes the divergences pathological to second-order Møller-Plesset (MP2) theory for strong correlations. In the language often used in density functional theory contexts, this means GF2 has a greatly reduced fractional spin error relative to MP2. The natural question then is what effect, if any, does the Dyson summation have on the fractional charge error in GF2? To this end, we generalize our previous implementation of GF2 to open-shell systems and analyze its fractional spin and charge errors. We find that like MP2, GF2 possesses only a very small fractional charge error, and consequently minimal many electron self-interaction error. This shows that GF2 improves on the critical failings of MP2, but without altering the positive features that make it desirable. Furthermore, we find that GF2 has both less fractional charge and fractional spin errors than typical hybrid density functionals as well as random phase approximation with exchange.

  5. Quantum rearrangement and self-consistent BCS-BEC crossover thermodynamics

    SciTech Connect

    Chen Jisheng; Li Jiarong; Wang Yanping; Cheng Chuanming

    2007-09-15

    Based on previous works, analytical calculational procedures for dealing with the strongly interacting fermion ground state are further developed through a medium-dependent potential in terms of the Bethe-Peierls contact interaction model. The methods are exact in the unitary limit regime and they lead to self-consistent equations analogous to the Hartree ones. The single particle energy spectrum rearrangement effects on the thermodynamics due to the Hugenholtz-van Hove theorem constraint are addressed. These effects lead to an additional instantaneous correlation potential contribution to the system physical chemical potential and pressure, i.e., equation of state, in order to enforce the classical thermodynamic consistency. The Dyson-Schwinger equations represent implicitly the various Bethe-Goldstone expansion ones. In a thermodynamically self-consistent way, the universal dimensionless factor is analytically calculated to be {xi}=(4/9), which defines the ratio of the unitary fermions energy density to that of the ideal noninteracting ones at T=0.

  6. Calculating gravitationally self-consistent sea level changes driven by dynamic topography

    NASA Astrophysics Data System (ADS)

    Austermann, J.; Mitrovica, J. X.

    2015-12-01

    We present a generalized formalism for computing gravitationally self-consistent sea level changes driven by the combined effects of dynamic topography, geoid perturbations due to mantle convection, ice mass fluctuations and sediment redistribution on a deforming Earth. Our mathematical treatment conserves mass of the surface (ice plus ocean) load and the solid Earth. Moreover, it takes precise account of shoreline migration and the associated ocean loading. The new formalism avoids a variety of approximations adopted in previous models of sea level change driven by dynamic topography, including the assumption that a spatially fixed isostatic amplification of `air-loaded' dynamic topography accurately accounts for ocean loading effects. While our approach is valid for Earth models of arbitrary complexity, we present numerical results for a set of simple cases in which a pattern of dynamic topography is imposed, the response to surface mass loading assumes that Earth structure varies only with depth and that isostatic equilibrium is maintained at all times. These calculations, involving fluid Love number theory, indicate that the largest errors in previous predictions of sea level change driven by dynamic topography occur in regions of shoreline migration, and thus in the vicinity of most geological markers of ancient sea level. We conclude that a gravitationally self-consistent treatment of long-term sea level change is necessary in any effort to use such geological markers to estimate ancient ice volumes.

  7. Self-consistent theory of nanodomain formation on non-polar surfaces of ferroelectrics

    DOE PAGESBeta

    Morozovska, Anna N.; Obukhovskii, Vyacheslav; Fomichov, Evhen; Varenyk, O. V.; Shur, Vladimir Ya.; Kalinin, Sergei V.; Eliseev, E. A.

    2016-04-28

    We propose a self-consistent theoretical approach capable of describing the features of the anisotropic nanodomain formation induced by a strongly inhomogeneous electric field of a charged scanning probe microscopy tip on nonpolar cuts of ferroelectrics. We obtained that a threshold field, previously regarded as an isotropic parameter, is an anisotropic function that is specified from the polar properties and lattice pinning anisotropy of a given ferroelectric in a self-consistent way. The proposed method for the calculation of the anisotropic threshold field is not material specific, thus the field should be anisotropic in all ferroelectrics with the spontaneous polarization anisotropy alongmore » the main crystallographic directions. The most evident examples are uniaxial ferroelectrics, layered ferroelectric perovskites, and low-symmetry incommensurate ferroelectrics. Obtained results quantitatively describe the differences at several times in the nanodomain length experimentally observed on X and Y cuts of LiNbO3 and can give insight into the anisotropic dynamics of nanoscale polarization reversal in strongly inhomogeneous electric fields.« less

  8. Microwave air plasmas in capillaries at low pressure I. Self-consistent modeling

    NASA Astrophysics Data System (ADS)

    Coche, P.; Guerra, V.; Alves, L. L.

    2016-06-01

    This work presents the self-consistent modeling of micro-plasmas generated in dry air using microwaves (2.45 GHz excitation frequency), within capillaries (<1 mm inner radius) at low pressure (300 Pa). The model couples the system of rate balance equations for the most relevant neutral and charged species of the plasma to the homogeneous electron Boltzmann equation. The maintenance electric field is self-consistently calculated adopting a transport theory for low to intermediate pressures, taking into account the presence of O‑ ions in addition to several positive ions, the dominant species being O{}2+ , NO+ and O+ . The low-pressure small-radius conditions considered yield very-intense reduced electric fields (∼600–1500 Td), coherent with species losses controlled by transport and wall recombination, and kinetic mechanisms strongly dependent on electron-impact collisions. The charged-particle transport losses are strongly influenced by the presence of the negative ion, despite its low-density (∼10% of the electron density). For electron densities in the range (1–≤ft. 4\\right)× {{10}12} cm‑3, the system exhibits high dissociation degrees for O2 (∼20–70%, depending on the working conditions, in contrast with the  ∼0.1% dissociation obtained for N2), a high concentration of O2(a) (∼1014 cm‑3) and NO(X) (5× {{10}14} cm‑3) and low ozone production (<{{10}-3}% ).

  9. Self-consistent calculation of hyperfine fields and adiabatic potential of impurities in iron

    NASA Astrophysics Data System (ADS)

    Kanamori, Junjiro; Akai, Hisazumi; Akai, Masako

    1984-01-01

    Hyperfine fields of impurities of the atomic number Z=1 56 at the substitutional site and those of light impurities of Z=1 9 at the interstitial sites in ferromagnetic iron are calculated by the KKR method adapted to the system containing a single impurity atom. The potential of the impurity atom is determined self-consistently by use of the local spin density functional formalism. The results for nonmagnetic sp valence impurities agree with those of the previous nonself-consistent calculation by Katayama-Yoshida, Terakura and Kanamori except for a few cases, confirming their theory of the systematic variation of hyperfine fields. The calculation for magnetic impurities of transition elements is presented for the first time in this paper. The calculations mentioned so far assume that impurities are situated at the center of each site. For the purpose of discussing the stability of the impurity positions, the change of the adiabatic potential due to displacements from the center is calculated by carrying out similar self-consistent calculations for off-center impurity positions. It is concluded that positive muon and some light impurities including boron will be displaced from the center when trapped in a vacancy.

  10. Efficient 3D/1D self-consistent integral-equation analysis of ICRH antennae

    NASA Astrophysics Data System (ADS)

    Maggiora, R.; Vecchi, G.; Lancellotti, V.; Kyrytsya, V.

    2004-08-01

    This work presents a comprehensive account of the theory and implementation of a method for the self-consistent numerical analysis of plasma-facing ion-cyclotron resonance heating (ICRH) antenna arrays. The method is based on the integral-equation formulation of the boundary-value problem, solved via a weighted-residual scheme. The antenna geometry (including Faraday shield bars and a recess box) is fairly general and three-dimensional (3D), and the plasma is in the one-dimensional (1D) 'slab' approximation; finite-Larmor radius effects, as well as plasma density and temperature gradients, are considered. Feeding via the voltages in the access coaxial lines is self-consistently accounted throughout and the impedance or scattering matrix of the antenna array obtained therefrom. The problem is formulated in both the dual space (physical) and spectral (wavenumber) domains, which allows the extraction and simple handling of the terms that slow the convergence in the spectral domain usually employed. This paper includes validation tests of the developed code against measured data, both in vacuo and in the presence of plasma. An example of application to a complex geometry is also given.