Quasilocal density functional theory and its application within the extended Thomas-Fermi approximation

In this paper we propose a generalization of the density functional theory. The theory leads to single-particle equations of motion with a quasilocal mean-field operator, which contains a quasiparticle position-dependent effective mass and a spin-orbit potential. The energy density functional is constructed using the extended Thomas-Fermi approximation and the ground-state properties of doubly magic nuclei are considered within the framework of this approach. Calculations were performed using the finite-range Gogny D1S forces and the results are compared with the exact Hartree-Fock calculations

Nuclear incompressibility in the quasilocal density functional theory

We explore the ability of the recently established quasilocal density functional theory for describing the isoscalar giant monopole resonance. Within this theory we use the scaling approach and perform constrained calculations for obtaining the cubic and inverse energy weighted moments (sum rules) of the RPA strength. The meaning of the sum rule approach in this case is discussed. Numerical calculations are carried out using Gogny forces and an excellent agreement is found with HF+RPA results previously reported in literature. The nuclear matter compression modulus predicted in our model lies in the range 210230 MeV which agrees with earlier findings. The information provided by the sum rule approach in the case of nuclei near the neutron drip line is also discussed.

An ab-initio calculation of the Coulomb explosion of N_2 after heavy-ion bombardement

Self-consistent-field calculations for the total potential energy of highly ionized N_2 molecules are presented. We compare these calculations to the experimentally observed energy released in the Coulomb explosion of ionized N_2 molecules created after collision with fast heavy ions. The most important electronic states of the fragment ions are determined.

Nuclear magnetic resonance chemical shifts with the statistical average of orbital-dependent model potentials in Kohn-Sham density functional theory

The performance of the SAOP potential for the calculation of NMR chemical shifts was evaluated. SAOP results show considerable improvement with respect to previous potentials, like VWN or BP86, at least for the carbon, nitrogen, oxygen, and fluorine chemical shifts. Furthermore, a few NMR calculations carried out on third period atoms (S, P, and Cl) improved when using the SAOP potential

Density functional energy decomposition into one- and two-atom contributions

The present work provides a generalization of Mayer's energy decomposition for the density-functional theory (DFT) case. It is shown that one- and two-atom Hartree-Fock energy components in Mayer's approach can be represented as an action of a one-atom potential VA on a one-atom density ρ A or ρ B. To treat the exchange-correlation term in the DFT energy expression in a similar way, the exchange-correlation energy density per electron is expanded into a linear combination of basis functions. Calculations carried out for a number of density functionals demonstrate that the DFT and Hartree-Fock two-atom energies agree to a reasonable extent with each other. The two-atom energies for strong covalent bonds are within the range of typical bond dissociation energies and are therefore a convenient computational tool for assessment of individual bond strength in polyatomic molecules. For nonspecific nonbonding interactions, the two-atom energies are low. They can be either repulsive or slightly attractive, but the DFT results more frequently yield small attractive values compared to the Hartree-Fock case. The hydrogen bond in the water dimer is calculated to be between the strong covalent and nonbonding interactions on the energy scale

Low-lying excited states and primary photoproducts of [Os-3(CO)(10)-(s-cis-L)] (L = cyclohexa-1,3-diene, buta-1,3-diene)] clusters studied by picosecond time-resolved UV/Vis and IR spectroscopy and by density functional theory

Combined picosecond transient absorption and time-resolved infrared studies were performed, aimed at characterising low-lying excited states of the cluster [Os-3(CO)(10)(s-cis-L)] (L= cyclohexa-1,3-diene, 1) and monitoring the formation of its photoproducts. Theoretical (DFT and TD-DFT) calculations on the closely related cluster with L=buta-1,3-diene (2') have revealed that the low-lying electronic transitions of these [Os-3(CO)(10)(s-cis-1,3-diene)] clusters have a predominant sigma(core)pi*(CO) character. From the lowest sigmapi* excited state, cluster 1 undergoes fast Os-Os(1,3-diene) bond cleavage (tau=3.3 ps) resulting in the formation of a coordinatively unsaturated primary photoproduct (1a) with a single CO bridge. A new insight into the structure of the transient has been obtained by DFT calculations. The cleaved Os-Os(1,3-diene) bond is bridged by the donor 1,3-diene ligand, compensating for the electron deficiency at the neighbouring Os centre. Because of the unequal distribution of the electron density in transient la, a second CO bridge is formed in 20 ps in the photoproduct [Os-3(CO)(8)(mu-CO)(2)- (cyclohexa-1,3-diene)] (1b). The latter compound, absorbing strongly around 630 nm, mainly regenerates the parent cluster with a lifetime of about 100 ns in hexane. Its structure, as suggested by the DFT calculations, again contains the 1,3-diene ligand coordinated in a bridging fashion. Photoproduct 1b can therefore be assigned as a high-energy coordination isomer of the parent cluster with all Os-Os bonds bridged.

Strongly nucleophilic RhI centre in square-planar complexes with terdentate (κ3) 2,2′:6′,2′′-terpyridine ligands: crystallographic, electrochemical and density functional theoretical studies

Rh-I-terpyridine complexes have been unambiguously formed for the first time. The 2,21:6',2"-terpyridine (tpy), 4'-chloro-2,2':6',2"-terpyridine (4'-Cl-tpy) and 4'-(tert-butyldimethylsilyl-ortho-carboranyl)-2,2':6',2"-terpyridine (carboranyl-tpy) ligands were used for successful syntheses and characterisation of the corresponding Rh-I complexes with halide coligands, [Rh(X)(4'-Y-terpyridine)] (X = Cl, Y = H, Cl, carboranyl; X = Br, Y = H). All four neutral Rh-tpy complexes are square planar, with Rh-X bonds in the plane of the 4'-Y-terpyridine ligands. Full characterisation of these dark blue, highly air-sensitive compounds was hampered by their poor solubility in various organic solvents. This is mainly due to the formation of pi-stacked aggregates, as evidenced by the crystal structure of [Rh(Cl)(tpy)]; in addition, [Rh(Cl)(carboranyl-tpy)] merely forms discrete dimers. The (bonding) properties of the novel Rh-I-terpyridine complexes have been studied with single-crystal X-ray diffraction, (time-dependent) density functional theoretical (DFT) calculations, far-infrared spectroscopy, electronic absorption spectroscopy and cyclic voltammetry. From DFT calculations, the HOMO of the studied Rh-I-terpyridine complexes involves predominantly the metal centre, while the LUMO resides on the terpyridine ligand. Absorption bands of the studied complexes in the visible region (400-900 nm) can be assigned to MLCT and MLCT/XLCT transitions. The relatively low oxidation potentials of [Rh(X)(tpy)] (X = Cl, Br) point to a high electron density on the metal centre. This makes the Rh-I-terpyridine complexes strongly nucleophilic and (potentially) highly reactive towards various (small) substrate molecules containing carbon-halide bonds.

Synthesis, chemical, electrochemical characterization of oxomolybdenum(V) complexes of 2-(3,5-dimethyl pyrazol-1-yl) benzothiazole ligand: crystal structure of ligand and oxomolybdenum(VI) complex and density functional theory (DFT) calculation

This work reports the ligational behavior of the neutral bidentate chelating molecule 2-(3,5-dimethyl pyrazol-1-yl) benzothiazole towards the oxomolybdenum(V) center. Both mononuclear complexes of the type (MoOX3L)-O-V and binuclear complexes of the formula (Mo2O4X2L2)-O-V (where X = Cl, Br) are isolated in the solid state. The complexes are characterized by elemental analyses, various spectroscopic techniques (UV-Vis IR), magnetic susceptibility measurement at room temperature, and cyclic voltammetry for their redox behavior at a platinum electrode in CH3CN. The mononuclear complexes (MoOX3L)-O-V are found to be paramagnetic while the binuclear complexes Mo2O4X2L2 are diamagnetic. Crystal and molecular structure of the ligand and the dioxomolybdenum complex (MoO2Br2L)-O-VI (obtained from the complex MoOBr3L during crystallization) have been solved by single crystal X-ray diffraction technique. Relevant DFT calculations of the ligand and the complex (MoO2Br2L)-O-VI are also carried out.

The interplay of secondary Hg⋯S, Hg⋯N and Hg⋯π bonding interactions in supramolecular structures of phenylmercury(ii) dithiocarbamates

Three new phenylmercury(II) and one mercury(II) dithiocarbamate complexes viz. PhHg S2CN(PyCH2) Bz (1), PhHg S2CN(PyCH2)CH3 (2), PhHg S2CN(Bz)CH3 (3), and [Hg (NCS2(PyCH2)Bz)(2)] (4) (Py = pyridine; Bz = benzyl) have been synthesized and characterized by elemental analyses, IR, electronic absorption, H-1 and C-13 NMR spectroscopy. The crystal structures of 1, 2 and 3 showed a linear S-Hg-C core at the centre of the molecule, in which the metal atom is bound to the sulfur atom of the dithiocarbamate ligand and a carbon atom of the aromatic ring. In contrast the crystal structure of 4 showed a linear S-Hg-S core at the Hg(II) centre of the molecule. Weak intermolecular Hg center dot center dot center dot N (Py) interactions link molecules into a linear chain in the case of 1, whereas chains of dimers are formed in 2 through intermolecular Hg center dot center dot center dot N (Py) and Hg center dot center dot center dot S interactions. 3 forms a conventional face-to-edge dimeric structure through intermolecular Hg center dot center dot center dot S secondary bonding and 4 forms a linear chain of dimers through face-to-face Hg center dot center dot center dot S secondary bonding. In order to elucidate the nature of these secondary bonding interactions and the electronic absorption spectra of the complexes, ab initio quantum chemical calculations at the MP2 level and density functional theory calculations were carried out for 1-3. Complexes 1 and 2 exhibited photoluminescent properties in the solid state as well as in the solution phase. Studies indicate that Hg center dot center dot center dot S interactions decrease and Hg center dot center dot center dot N interactions increase the chances of photoluminescence in the solid phase

Cobalt incorporation in calcite: Thermochemistry of (Ca,Co)CO3 solid solutions from density functional theory simulations

The incorporation of cobalt in mixed metal carbonates is a possible route to the immobilization of this toxic element in the environment. However, the thermodynamics of (Ca,Co)CO3 solid solutions are still unclear due to conflicting data from experiment and from the observation of natural ocurrences. We report here the results of a computer simulation study of the mixing of calcite (CaCO3) and spherocobaltite (CoCO3), using density functional theory calculations. Our simulations suggest that previously proposed thermodynamic models, based only on the range of observed compositions, significantly overestimate the solubility between the two solids and therefore underestimate the extension of the miscibility gap under ambient conditions. The enthalpy of mixing of the disordered solid solution is strongly positive and moderately asymmetric: calcium incorporation in spherocobaltite is more endothermic than cobalt incorporation in calcite. Ordering of the impurities in (0001) layers is energetically favourable with respect to the disordered solid solution at low temperatures and intermediate compositions, but the ordered phase is still unstable to demixing. We calculate the solvus and spinodal lines in the phase diagram using a sub-regular solution model, and conclude that many Ca1-xCoxCO3 mineral solid solutions (with observed compositions of up to x=0.027, and above x=0.93) are metastable with respect to phase separation. We also calculate solid/aqueous distribution coefficients to evaluate the effect of the strong non-ideality of mixing on the equilibrium with aqueous solution, showing that the thermodynamically-driven incorporation of cobalt in calcite (and of calcium in spherocobaltite) is always very low, regardless of the Co/Ca ratio of the aqueous environment.

Ab initio study of SO2 molecules interacting with pristine and transition metal covered fullerenes as a possible route for nanofilters

In this work, a systematic study of SO2 molecules interacting with pristine and transition metal (TM) covered C-60 is presented by means of first principles calculations. It is observed that the SO2 molecule interacts weakly with the pristine C-60 fullerene, although the resulting interaction is largely increased when the C-60 structure is covered with Fe, Mn, or Ti atoms and the SO2 Molecules are bounded through the TM atoms. The number of bounded SO2 molecules per TM atoms, in addition to the elevated binding energies per molecules, allows us to conclude that such composites can be used as a template for efficient devices to remove SO2 molecules or, alternatively, as SO2 gas sensor.

DFT study of the electronic, vibrational, and optical properties of SnO(2)

We report results on the electronic, vibrational, and optical properties of SnO(2) obtained using first-principles calculations performed within the density functional theory. All the calculated phonon frequencies, real and imaginary parts of complex dielectric function, the energy-loss spectrum, the refractive index, the extinction, and the absorption coefficients show good agreement with experimental results. Based on our calculations, the SnO(2) electron and hole effective masses were found to be strongly anisotropic. The lattice contribution to the low-frequency region of the SnO(2) dielectric function arising from optical phonons was also determined resulting the values of E > (1aSyen) (latt) (0) = 14.6 and E > (1ayen) (latt) (0) = 10.7 for directions perpendicular and parallel to the tetragonal c-axis, respectively. This is in excellent agreement with the available experimental data. After adding the electronic contribution to the lattice contribution, a total average value of E >(1)(0) = 18.2 is predicted for the static permittivity constant of SnO(2).

Ab initio calculation of the BCC Fe-Al-Mo (Iron-Aluminum-Molybdenum) phase diagram: Implications for the nature of the tau(2) phase

The metastable phase diagram of the BCC-based ordering equilibria in the Fe-Al-Mo system has been calculated via a truncated cluster expansion, through the combination of Full-Potential-Linear augmented Plane Wave (FP-LAPW) electronic structure calculations and of Cluster Variation Method (CVM) thermodynamic calculations in the irregular tetrahedron approximation. Four isothermal sections at 1750 K, 2000 K, 2250 K and 2500 K are calculated and correlated with recently published experimental data on the system. The results confirm that the critical temperature for the order-disorder equilibrium between Fe(3)Al-D0(3) and FeAl-B2 is increased by Mo additions, while the critical temperature for the FeAl-B2/A2 equilibrium is kept approximately invariant with increasing Mo contents. The stabilization of the Al-rich A2 phase in equilibrium with overstoichiometric B2-(Fe,Mo)Al is also consistent with the attribution of the A2 structure to the tau(2) phase, stable at high temperatures in overstoichiometric B2-FeAl. (C) 2009 Elsevier Ltd. All rights reserved.

Theoretical investigation of electronic and optical properties of andalusite within density functional theory

The electronic and optical properties of andalusite were studied by using quantum-mechanical calculations based on the density functional theory (DFT). The electronic structure shows that andalusite has a direct band gap of 5.01 eV. The complex dielectric function and optical constants, such as extinction coefficient, refractive index, reflectivity and energy-loss spectrum, are calculated. The optical properties of andalusite are discussed based on the band structure calculations. It is shown that the O-2p states and Al-3s states play a major role in optical transitions as initial and final states, respectively. (C) 2010 Elsevier Ltd. All rights reserved.

Large Changes of Static Electric Properties Induced by Hydrogen Bonding: An ab Initio Study of Linear HCN Oligomers

We report the partitioning of the interaction-induced static electronic dipole (hyper)polarizabilities for linear hydrogen cyanide complexes into contributions arising from various interaction energy terms. We analyzed the nonadditivities of the studied properties and used these data to predict the electric properties of an infinite chain. The interaction-induced static electric dipole properties and their nonadditivities were analyzed using an approach based on numerical differentiation of the interaction energy components estimated in an external electric field. These were obtained using the hybrid variational-perturbational interaction energy decomposition scheme, augmented with coupled-cluster calculations, with singles, doubles, and noniterative triples. Our results indicate that the interaction-induced dipole moments and polarizabilities are primarily electrostatic in nature; however, the composition of the interaction hyperpolarizabilities is much more complex. The overlap effects substantially quench the contributions due to electrostatic interactions, and therefore, the major components are due to the induction and exchange induction terms, as well as the intramolecular electron-correlation corrections. A particularly intriguing observation is that the interaction first hyperpolarizability in the studied systems not only is much larger than the corresponding sum of monomer properties, but also has the opposite sign. We show that this effect can be viewed as a direct consequence of hydrogen-bonding interactions that lead to a decrease of the hyperpolarizability of the proton acceptor and an increase of the hyperpolarizability of the proton donor. In the case of the first hyperpolarizability, we also observed the largest nonadditivity of interaction properties (nearly 17%) which further enhances the effects of pairwise interactions.