Coulomb-hole summations and energies for GW calculations with limited number of empty orbitals: a modified static remainder approach

Ab initio GW calculations are a standard method for computing the spectroscopic properties of many materials. The most computationally expensive part in conventional implementations of the method is the generation and summation over the large number of empty orbitals required to converge the electron self-energy. We propose a scheme to reduce the summation over empty states by the use of a modified static remainder approximation, which is simple to implement and yields accurate self-energies for both bulk and molecular systems requiring a small fraction of the typical number of empty orbitals.

Valence band electronic structure of Nd1-xYxMnO3 using X-ray absorption, photoemission and GGA plus U calculations

The electronic structures of Nd1-xYxMnO3 (x=0-0.5) were studied using X-ray absorption near-edge structure (XANES) at the Mn L-3,L-2- and O K-edge along with valence-band photoemission spectroscopy (VB-PES). The systematic increase in white-line intensity of the Mn L-3,L-2-edge with doping, suggests a decrease in the occupancy of Mn 3d orbitals. The O K-edge XANES shows a depletion of unoccupied states above the Fermi energy. The changes in the O K-edge spectra due to doping reflects an increase in the Jahn-Teller distortion. The VB-PES shows broadening of the features associated with Mn 3d and O 2p hybridized states and the shift of these features to a slightly higher binding energy in agreement with our GGA + U calculations. The system shows a net shift of the occupied and unoccupied states away from the Fermi energy with doping. The shift in theoretical site-projected density of states of x=0.5 composition with respect to x=0 suggest a subtle change from a charge transfer to Mott-Hubbard type insulator. (C) 2013 Elsevier B.V. All rights reserved.

Improved quasiparticle wave functions and mean field for G(0)W(0) calculations: Initialization with the COHSEX operator

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 G W self-energy operator, E, is taken to be diagonal in the density functional theory (DFT) Kohn-Sham basis within the G0 W0 scheme. However, there are known situations in which this diagonal Go Wo approximation starting from DFT is inadequate. We present two schemes to resolve such problems. The first, which we called sc-COHSEX-PG W, 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 W. In this scheme, frequency-dependent self energy E(N), 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 G W (od-COHSEX-PG W). 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 G W E 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-PGW method solves the qualitatively incorrect DFT mean-field starting point (having a band overlap) in bulk Ge under pressure.

Intramolecular hydrogen bonds involving organic fluorine in the derivatives of hydrazides: an NMR investigation substantiated by DFT based theoretical calculations

The rare examples of intramolecular hydrogen bonds (HB) of the type the N-H center dot center dot center dot F-C, detected in a low polarity solvent in the derivatives of hydrazides, by utilizing one and two-dimensional solution state multinuclear NMR techniques, are reported. The observation of through-space couplings, such as, (1h)J(FH), and (1h)J(FN), provides direct evidence for the existence of intra-molecular HB. Solvent induced perturbations and the variable temperature NMR experiments unambiguously establish the presence of intramolecular HB. The existence of multiple conformers in some of the investigated molecules is also revealed by two dimensional HOESY and N-15-H-1 HSQC experiments. The H-1 DOSY experimental results discard any possibility of self or cross dimerization of the molecules. The derived NMR experimental results are further substantiated by Density Function Theory (DFT) based Non Covalent Interaction (NCI), and Quantum Theory of Atom in Molecule (QTAIM) calculations. The NCI calculations served as a very sensitive tool for detection of non-covalent interactions and also confirm the presence of bifurcated HBs.

Organic fluorine involved intramolecular hydrogen bonds in the derivatives of imides: NMR evidence corroborated by DFT based theoretical calculations

The rare occurrence of intramolecular hydrogen bonds (HBs) of the type N-H center dot center dot center dot F-C is detected in the derivatives of imides in a low polarity solvent by using multi-dimensional and multinuclear NMR experiments. The observation of (1h)J(FH), (2h)J(FN), and (2h)J(FF), where the spin magnetization is transmitted through space among the interacting NMR active nuclei, provided strong and unambiguous evidence for the existence of intra-molecular HBs. The variation in the chemical shifts of labile protons depending on physical conditions, such as the solvent dilution and the systematic alteration of temperature confirmed the presence of weak interactions through intramolecular HBs in all the investigated fluorine substituted molecules. The self or cross dimerization of molecules is unequivocally discarded by the analysis of the rates of diffusion obtained using pseudo-two dimensional DOSY experiments. The Density Function Theory (DFT) calculations based on the Quantum Theory of Atoms In Molecules (QTAIM) and Non Covalent Interaction (NCI), are in close agreement with the NMR experimental findings.

H2O-CH4 and H2S-CH4 complexes: a direct comparison through molecular beam experiments and ab initio calculations

New molecular beam scattering experiments have been performed to measure the total ( elastic plus inelastic) cross sections as a function of the velocity in collisions between water and hydrogen sulfide projectile molecules and the methane target. Measured data have been exploited to characterize the range and strength of the intermolecular interaction in such systems, which are of relevance as they drive the gas phase molecular dynamics and the clathrate formation. Complementary information has been obtained by rotational spectra, recorded for the hydrogen sulfide-methane complex, with a pulsed nozzle Fourier transform microwave spectrometer. Extensive ab initio calculations have been performed to rationalize all the experimental findings. The combination of experimental and theoretical information has established the ground for the understanding of the nature of the interaction and allows for its basic components to be modelled, including charge transfer, in these weakly bound systems. The intermolecular potential for H2S-CH4 is significantly less anisotropic than for H2O-CH4, although both of them have potential minima that can be characterized as `hydrogen bonded'.

Ray tracing based path-length calculations for polarized light tomographic imaging

A ray tracing based path length calculation is investigated for polarized light transport in a pixel space. Tomographic imaging using polarized light transport is promising for applications in optical projection tomography of small animal imaging and turbid media with low scattering. Polarized light transport through a medium can have complex effects due to interactions such as optical rotation of linearly polarized light, birefringence, diattenuation and interior refraction. Here we investigate the effects of refraction of polarized light in a non-scattering medium. This step is used to obtain the initial absorption estimate. This estimate can be used as prior in Monte Carlo (MC) program that simulates the transport of polarized light through a scattering medium to assist in faster convergence of the final estimate. The reflectance for p-polarized (parallel) and s-polarized (perpendicular) are different and hence there is a difference in the intensities that reach the detector end. The algorithm computes the length of the ray in each pixel along the refracted path and this is used to build the weight matrix. This weight matrix with corrected ray path length and the resultant intensity reaching the detector for each ray is used in the algebraic reconstruction (ART) method. The proposed method is tested with numerical phantoms for various noise levels. The refraction errors due to regions of different refractive index are discussed, the difference in intensities with polarization is considered. The improvements in reconstruction using the correction so applied is presented. This is achieved by tracking the path of the ray as well as the intensity of the ray as it traverses through the medium.

Timing Calculations for a General N-Level Dodecagonal Space Vector Structure Using Only Reference Phase Voltages

Multilevel inverters with dodecagonal (12-sided polygon) voltage space vector (SV) structures have advantages like extension of linear modulation range, elimination of fifth and seventh harmonics in phase voltages and currents for the full modulation range including extreme 12-step operation, reduced device voltage ratings, lesser dv/dt stresses on devices and motor phase windings resulting in lower EMI/EMC problems, and lower switching frequency-making it more suitable for high-power drive applications. This paper proposes a simple method to obtain pulsewidth modulation (PWM) timings for a dodecagonal voltage SV structure using only sampled reference voltages. In addition to this, a carrier-based method for obtaining the PWM timings for a general N-level dodecagonal structure is proposed in this paper for the first time. The algorithm outputs the triangle information and the PWM timing values which can be set as the compare values for any carrier-based hardware PWM module to obtain SV PWM like switching sequences. The proposed method eliminates the need for angle estimation, computation of modulation indices, and iterative search algorithms that are typical in multilevel dodecagonal SV systems. The proposed PWM scheme was implemented on a five-level dodecagonal SV structure. Exhaustive simulation and experimental results for steady-state and transient conditions are presented to validate the proposed method.

Multiphase technique improves mud-pulse velocity calculations

Measurement while drilling (MWD) has become a popular survey technology to monitor directional data, drilling data, formation evaluation data and safety data in the world. And closed loop drilling shows promise in recent years. Obviously, the method of tr

Elasticity, stability, and ideal strength of beta-SiC in plane-wave-based ab initio calculations

On the basis of the pseudopotential plane-wave method and the local-density-functional theory, this paper studies energetics, stress-strain relation, stability, and ideal strength of beta-SiC under various loading modes, where uniform uniaxial extension and tension and biaxial proportional extension are considered along directions [001] and [111]. The lattice constant, elastic constants, and moduli of equilibrium state are calculated and the results agree well with the experimental data. As the four SI-C bonds along directions [111], [(1) over bar 11], [11(1) over bar] and [111] are not the same under the loading along [111], internal relaxation and the corresponding internal displacements must be considered. We find that, at the beginning of loading, the effect of internal displacement through the shuffle and glide plane diminishes the difference among the four Si-C bonds lengths, but will increase the difference at the subsequent loading, which will result in a crack nucleated on the {111} shuffle plane and a subsequently cleavage fracture. Thus the corresponding theoretical strength is 50.8 GPa, which agrees well with the recent experiment value, 53.4 GPa. However, with the loading along [001], internal relaxation is not important for tetragonal symmetry. Elastic constants during the uniaxial tension along [001] are calculated. Based on the stability analysis with stiffness coefficients, we find that the spinodal and Born instabilities are triggered almost at the same strain, which agrees with the previous molecular-dynamics simulation. During biaxial proportional extension, stress and strength vary proportionally with the biaxial loading ratio at the same longitudinal strain.

Calculations of longitudinal and transverse velocity structure functions using a vortex model of isotropic turbulence

The longitudinal structure function (LSF) and the transverse structure function (TSF) in isotropic turbulence are calculated using a vortex model. The vortex model is composed of the Rankine and Burgers vortices which have the exponential distributions in the vortex Reynolds number and vortex radii. This model exhibits a power law in the inertial range and satisfies the minimal condition of isotropy that the second-order exponent of the LSF in the inertial range is equal to that of the TSF. Also observed are differences between longitudinal and transverse structure functions caused by intermittency. These differences are related to their scaling differences which have been previously observed in experiments and numerical simulations.

Analytical calculations of Eulerian and Lagrangian time correlations in turbulent shear flows

The Taylor series expansion method is used to analytically calculate the Eulerian and Lagrangian time correlations in turbulent shear flows. The short-time behaviors of those correlation functions can be obtained from the series expansions. Especially, the propagation velocity and sweeping velocity in the elliptic model of space-time correlation are analytically calculated and further simplified using the sweeping hypothesis and straining hypothesis. These two characteristic velocities mainly determine the space-time correlations.