Acetone and Toluene Gas Sensing by WO3: Focusing on the Selectivity from First Principle Calculations.

Sensitivity and selectivity are the two major parameters that should be optimized in chemiresistive devices with boosted performances towards Volatile Organic Compounds (VOCs). Notwithstanding a plethora of metal oxides/VOCs combinations that have been investigated so far, a close inspection based on theoretical models to provide guidelines to enhance sensors features has been scarcely explored. In this work, we measured experimentally the sensor response of a WO3 chemiresistor towards gaseous acetone and toluene, observing a two orders of magnitude higher signal for the former. In order to gain insight on the observed selectivity, Density Functional Theory was then adopted to elucidate how acetone and toluene molecules adsorption may perturb the electronic structure of WO3 due to electrostatic interactions with the surface and hybridization with its electronic structure. The results of acetone adsorption suggest the activation of the carbonyl group for reactions, while an overall lower charge redistribution on the surface and the molecule was observed for toluene. This, combined with acetone's higher binding energy, justifies the difference in the final responses. Notably, the presence of surface oxygen vacancies, characterizing the nanostructure of the oxide, influences the sensing performances.

Ab Initio Many-Body Perturbation Theory Calculations of the Electronic and Optical Properties of Cyclometalated Ir(III) Complexes.

Cyclometalated Ir(III) compounds are the preferred choice as organic emitters in organic light-emitting diodes. In practice, the presence of the transition metal surrounded by carefully designed ligands allows fine-tuning of the emission frequency as well as good efficiency of the device. To support the development of new compounds, experimental measurements are generally compared with absorption and emission spectra obtained from ab initio calculations. The standard approach for these calculations is time-dependent density functional theory (TDDFT) with a hybrid exchange-correlation functional like B3LYP. Because of the size of these compounds, the application of more complex quantum chemistry approaches can be challenging. In this work, we used many-body perturbation theory approaches, in particular the GW approximation with the Bethe-Salpeter equation (BSE) implemented in Gaussian basis sets, to calculate the quasiparticle properties and the absorption spectra of six cyclometalated Ir(III) complexes, going beyond TDDFT. In the presented results, we compared standard TDDFT simulations with BSE calculations performed on top of perturbative G0W0 and accounting for eigenvalue self-consistency. Moreover, in order to investigate in detail the effect of the DFT starting point, we concentrated on Ir(ppy)3 and performed GW-BSE simulations starting from different DFT exchange-correlation potentials.

Magnetic Moments and Electron Transport through Chromium-Based Antiferromagnetic Nanojunctions.

We report the electronic, magnetic and transport properties of a prototypical antiferromagnetic (AFM) spintronic device. We chose Cr as the active layer because it is the only room-temperature AFM elemental metal. We sandwiched Cr between two non-magnetic metals (Pt or Au) with large spin-orbit coupling. We also inserted a buffer layer of insulating MgO to mimic the structure and finite resistivity of a real device. We found that, while spin-orbit has a negligible effect on the current flowing through the device, the MgO layer plays a crucial role. Its effect is to decouple the Cr magnetic moment from Pt (or Au) and to develop an overall spin magnetization. We have also calculated the spin-polarized ballistic conductance of the device within the Büttiker⁻Landauer framework, and we have found that for small applied bias our Pt/Cr/MgO/Pt device presents a spin polarization of the current amounting to ≃25%.

Helium-induced electronic transitions in photo-excited Ba+-Hen exciplexes.

The possibility for helium-induced electronic transitions in a photo-excited atom is investigated using Ba+ excited to the 6p 2P state as a prototypical example. A diabatization scheme has been designed to obtain the necessary potential energy surfaces and couplings for complexes of Ba+ with an arbitrary number of helium atoms. It involves computing new He-Ba+ electronic wave functions and expanding them in determinants of the non-interacting complex. The 6p 2P ← 6s 2S photodissociation spectrum of He⋯Ba+ calculated with this model shows very weak coupling for a single He atom. However, several electronic relaxation mechanisms are identified, which could potentially explain the expulsion of barium ions from helium nanodroplets observed experimentally upon Ba+ photoexcitation. For instance, an avoided crossing in the ring-shaped He7Ba+ structure is shown to provide an efficient pathway for fine structure relaxation. Symmetry breaking by either helium density fluctuations or vibrations can also induce efficient relaxation in these systems, e.g., bending vibrations in the linear He2Ba+ excimer. The identified relaxation mechanisms can provide insight into helium-induced non-adiabatic transitions observed in other systems.

Dynamics of photoexcited Ba(+) cations in (4)He nanodroplets.

We present a joint experimental and theoretical study on the desolvation of Ba(+) cations in (4)He nanodroplets excited via the 6p ← 6s transition. The experiments reveal an efficient desolvation process yielding mainly bare Ba(+) cations and Ba(+)Hen exciplexes with n = 1 and 2. The speed distributions of the ions are well described by Maxwell-Boltzmann distributions with temperatures ranging from 60 to 178 K depending on the excitation frequency and Ba(+) Hen exciplex size. These results have been analyzed by calculations based on a time-dependent density functional description for the helium droplet combined with classical dynamics for the Ba(+). In agreement with experiment, the calculations reveal the dynamical formation of exciplexes following excitation of the Ba(+) cation. In contrast to experimental observation, the calculations do not reveal desolvation of excited Ba(+) cations or exciplexes, even when relaxation pathways to lower lying states are included.

Exciplexes with ionic dopants: stability, structure, and experimental relevance of M(+)((2)P)(4)He(n) (M = Sr, Ba).

M(+)((2)P)(4)Hen species, possibly involved in the post (2)P ← (2)S excitation dynamics of Sr(+) and Ba(+) in cold (4)He gas or droplets, are studied employing both high level ab initio calculations to determine the potential energy curves (PEC) and diffusion Monte Carlo (DMC) to obtain information on their ground state structure and energetics. PEC for the excited M(+)((2)P)He dimers were obtained using MRCI calculations with extended basis sets. Potential energy surfaces (PES) for M(+)((2)P)Hen were built with the DIM model including spin-orbit coupling via a perturbative procedure. DMC simulations indicated several exciplexes (n > 2) to be stable against He dissociation whatever the ion state, a finding that is at variance with what was previously suggested for Ba(+)((2)P1/2) due to the repulsive nature of the interaction potential obtained in [ Phys. Rev. A 2004 , 69 , 042505 ]. Our results, instead, support the suggestion made in [J. Chem. Phys. 2012, 137, 051102] for the existence of Ba(+)((2)P1/2)Hen exciplexes emitted following the excitation of the barium cation solvated into He droplets. In the (2)P1/2 state, the Ba ion also shows a peculiar behavior as a function of n with respect to the location and binding strength of the attached He atoms compared to Sr(+). Although the latter forms the usual equatorial He ring, Ba(+) deviates from this geometry for 1 ≤ n ≤ 4, with the radial distribution functions strongly depending on the number of solvent atoms. Finally, a putative species is proposed to explain the emission band at 523 nm that follows D1 or D2 excitation of Ba(+) in superfluid bulk helium.

Communication: nucleation of quantized vortex rings in 4He nanodroplets.

Whereas most of the phenomena associated with superfluidity have been observed in finite-size helium systems, the nucleation of quantized vortices has proven elusive. Here we show using time-dependent density functional simulations that the solvation of a Ba(+) ion created by photoionization of neutral Ba at the surface of a (4)He nanodroplet leads to the nucleation of a quantized ring vortex. The vortex is nucleated on a 10 ps timescale at the equator of a solid-like solvation structure that forms around the Ba(+) ion. The process is expected to be quite general and very efficient under standard experimental conditions.

Coinage metal exciplexes with helium atoms: a theoretical study of M*(2L)He(n) (M = Cu, Ag, Au; L = P,D).

The structure and energetics of exciplexes M*((2)L)He(n) (M = Cu, Ag and Au; L = P and D) in their vibrational ground state are studied by employing diffusion Monte Carlo (DMC). Interaction potentials between the excited coinage metals and He atoms are built using the Diatomics-in-Molecule (DIM) approach and ab initio potential curves for the M((2)L)-He dimers. Extending our previous work [Cargnoni et al., J. Phys. Chem. A, 2011, 115, 7141], we computed the dimer potential for Au in the (2)P and (2)D states, as well for Cu and Ag in the (2)D state, employing basis set superposition error-corrected Configuration Interaction calculations. We found that the (2)Π potential correlating with the (2)P state of Au is substantially less binding than for Ag and Cu, a trend well supported by the M(+) ionic radiuses. Conversely, the interaction potentials between a (n - 1)d(9)ns(2 2)D metal and He present a very weak dependency on M itself or the projection of the angular momentum along the dimer axis. This is due to the screening exerted by the ns(2) electrons on the hole in the (n - 1)d shell. Including the spin-orbit coupling perturbatively in the DIM energy matrix has a major effect on the lowest potential energy surface of the (2)P manifold, the one for Cu allowing the formation of a "belt" of five He atoms while the one for Au being completely repulsive. Conversely, spin-orbit coupling has only a weak effect on the (2)D manifold due to the nearly degenerate nature of the diatomic potentials. Structural and energetic results from DMC have been used to support experimental indications for the formation of metastable exciplexes or the opening of non-radiative depopulation channels in bulk and cold gaseous He.

Solubility of metal atoms in helium droplets: exploring the effect of the well depth using the coinage metals Cu and Ag.

We report a theoretical investigation of the solution properties of Cu and Ag atoms dissolved in He clusters. Employing our recent ab initio ground state pair potential for Me-He (Me = Ag, Cu), we simulated the species Me@He (n) (n = 2-100) by means of diffusion Monte Carlo (DMC) obtaining exact information on their energetics and the structural properties. In particular, we investigated the sensitivity of structural details on the well depth of the two interaction potentials. Whereas Ag structures the first He solvation layer similarly, to some extent, to a positive ion such as Na(+), Cu appears to require the onset of a second solvation shell for a similar dense structure to be formed despite an interaction well of 28.4 µhartree. An additional signature of the different solution behavior between Ag and Cu appears also in the dependence of the energy required to evaporate a single He atom on the size of the MeHe(n) clusters. The absorption spectrum for the (2)P ← (2)S excitation of the metals was also simulated employing the semi-classical Lax approximation to further characterize Me@He(n) (n = 2-100) using novel accurate interaction potentials between He and the lowest (2)P state of Ag and Cu in conjunction with the Diatomic-in-Molecules approach. The results indicated that Ag exciplexes should not form via a direct vertical excitation into an attractive region of the excited manifolds and that there is an interesting dependence of the shape of the Cu excitation bands on the local structure of the first solvation shell.

Results and perspectives of the MO-VB method. Application examples on the He2 and the LiH-He complexes.

We present a short overview and two benchmark applications of the Molecular Orbital-Valence Bond (MO-VB) method, a quantum mechanical scheme developed within the framework of modern Valence Bond theory. The MO-VB has been especially designed to deal with weak intermolecular interactions, and in recent years it has been successfully applied to compute the potential energy surface of several complexes, namely He(2), He-H(2)O, He-CH(4) and Ne-CH(4). In this investigation we test extensively the performance of the MO-VB on two limit systems, the He dimer and the LiH-He complex, which span three extremes in the field of van der Waals interactions: a purely dispersive system (He...He) and, thanks to the highly polar Li-H bond, a noble gas approaching a cation (He...Li(+)), and an anion (He...H(-)). Very accurate computations are available in the literature for He(2) and LiH-He, performed either with conventional Molecular Orbital or Monte Carlo approaches, and hence these systems are the ideal candidates to establish the capabilities of the MO-VB. Based on the results of our study, we conclude that in He(2) and LiH-He the MO-VB overestimates the correlation energy contribution to the interaction energy of about 5%, and hence it is a valid option for computing accurate lower bounds to the potential energy surface of these complexes.

Ground state potential energy surfaces and bound states of M-He dimers (M=Cu,Ag,Au): a theoretical investigation.

We present an ab initio investigation on the ground state interaction potentials [potential energy surface (PES)] between helium and the group 11 metal atoms: copper, silver, and gold. To the best of our knowledge, there are no previous theoretical PESs proposed for Cu-He and Au-He, and a single one for Ag-He [Z. J. Jakubek and M. Takami, Chem. Phys. Lett. 265, 653 (1997)], computed about 10 years ago at MP2 level and significantly improved by our study. To reach a high degree of accuracy in the determination of the three M-He potentials (M=Cu,Ag,Au), we performed extensive series of test computations to establish the appropriate basis set, the theoretical method, and the computational scheme for these systems. For each M-He dimer we computed the PES at the CCSD(T) level of theory, starting from the reference unrestricted Hartree-Fock wave function. We described the inner shells with relativistic small core pseudopotentials, and we adopted high quality basis sets for the valence electrons. We also performed CCSDT computations in a limited set of M-He internuclear distances, adopting a medium-sized basis set, such as to define for each dimer a CCSD(T) to CCSDT correction term and to improve further the quality of the CCSD(T) interaction potentials. The Cu-He complex has minimum interaction energy (E(min)) of -28.4 microhartree at the internuclear distance of 4.59 A (R(min)), and the short-range repulsive wall starts at 4.04 A (R(E=0)). Quite interestingly, the PES of Ag-He is more attractive (E(min)=-33.8 microhartree) but presents nearly the same R(min) and R(E=0) values, 4.60 and 4.04 A, respectively. The interaction potential for Au-He is markedly deeper and shifted at shorter distances as compared to the lighter complexes, with E(min)=-69.6 microhartree, R(min)=4.09 A and R(E=0)=3.60 A. As a first insight in the structure of M-He(n) aggregates, we determined the rovibrational structure of the three M-He dimers. The Cu-He and Ag-He potentials support just few rotational excitations, while the Au-He PES admits also a bound vibrational excitation.

The electronic structure of nitrilimine: absence of the carbenic form.

The electronic structure of nitrilimine HCNNH is shown to essentially be propargylic by CASSCF and Spin-Coupled (modern VB) calculations; in contrast to a recent claim, the carbenic resonance form is absent.

Predicting atomic dopant solvation in helium clusters: the MgHe(n) case.

We present a quantum Monte Carlo study of the solvation and spectroscopic properties of the Mg-doped helium clusters MgHe(n) with n=2-50. Three high-level [MP4, CCSD(T), and CCSDT] MgHe interaction potentials have been used to study the sensitivity of the dopant location on the shape of the pair interaction. Despite the similar MgHe well depth, the pair distribution functions obtained in the diffusion Monte Carlo simulations markedly differ for the three pair potentials, therefore indicating different solubility properties for Mg in He(n). Moreover, we found interesting size effects for the behavior of the Mg impurity. As a sensitive probe of the solvation properties, the Mg excitation spectra have been simulated for various cluster sizes and compared with the available experimental results. The interaction between the excited 1P Mg atom and the He moiety has been approximated using the diatomics-in-molecules method and the two excited 1pi and 1sigma MgHe potentials. The shape of the simulated MgHe50 spectra shows a substantial dependency on the location of the Mg impurity, and hence on the MgHe pair interaction employed. To unravel the dependency of the solvation behavior on the shape of the computed potentials, exact density-functional theory has been adapted to the case of doped He(n) and various energy distributions have been computed. The results indicate the shape of the repulsive part of the MgHe potential as an important cause of the different behaviors.

Interstitial Zn atoms do the trick in thermoelectric zinc antimonide, Zn4Sb3: a combined maximum entropy method X-ray electron density and ab initio electronic structure study.

The experimental electron density of the high-performance thermoelectric material Zn4Sb3 has been determined by maximum entropy (MEM) analysis of short-wavelength synchrotron powder diffraction data. These data are found to be more accurate than conventional single-crystal data due to the reduction of common systematic errors, such as absorption, extinction and anomalous scattering. Analysis of the MEM electron density directly reveals interstitial Zn atoms and a partially occupied main Zn site. Two types of Sb atoms are observed: a free spherical ion (Sb3-) and Sb2(4-) dimers. Analysis of the MEM electron density also reveals possible Sb disorder along the c axis. The disorder, defects and vacancies are all features that contribute to the drastic reduction of the thermal conductivity of the material. Topological analysis of the thermally smeared MEM density has been carried out. Starting with the X-ray structure ab initio computational methods have been used to deconvolute structural information from the space-time data averaging inherent to the XRD experiment. The analysis reveals how interstitial Zn atoms and vacancies affect the electronic structure and transport properties of beta-Zn4Sb3. The structure consists of an ideal A12Sb10 framework in which point defects are distributed. We propose that the material is a 0.184:0.420:0.396 mixture of A12Sb10, A11BCSb10 and A10BCDSb10 cells, in which A, B, C and D are the four Zn sites in the X-ray structure. Given the similar density of states (DOS) of the A12Sb10, A11BCSb10 and A10BCDSb10 cells, one may electronically model the defective stoichiometry of the real system either by n-doping the 12-Zn atom cell or by p-doping the two 13-Zn atom cells. This leads to similar calculated Seebeck coefficients for the A12Sb10, A11BCSb10 and A10BCDSb10 cells (115.0, 123.0 and 110.3 microV K(-1) at T=670 K). The model system is therefore a p-doped semiconductor as found experimentally. The effect is dramatic if these cells are doped differently with respect to the experimental electron count. Thus, 0.33 extra electrons supplied to either kind of cell would increase the Seebeck coefficient to about 260 microV K(-1). Additional electrons would also lower sigma, so the resulting effect on the thermoelectric figure of merit of Zn4Sb3 challenges further experimental work.

Chemical information from the source function.

The source function, which enables one to equate the value of the electron density at any point within a molecule to a sum of atomic contributions, has been applied to a number of cases. The source function is a model-independent, quantitative measure of the relative importance of an atom's or group's contribution to the density at any point in a system, and it represents a potentially interesting tool to provide chemical information. It is shown that the source contribution from H to the electron density rho(b) at the bond critical point in HX diatomics decreases with increasing X's electronegativity, and that this decrease is a result of significant changes in the Laplacian distribution within the H-basin. It is also demonstrated that the source function from Li to rho(b) in LiX diatomics is a more sensitive index of atomic transferability than it is the lithium atomic energy or population. The observed changes are such as to ensure a constant percentage source contribution from Li to rho(b) throughout the LiX series, rather than a constant source as one would expect in the limit of perfect atomic transferability. Application of the source function to planar lithium clusters has revealed that the source function clearly discriminates between a nonnuclear electron density maximum and a maximum associated to a nucleus, on the basis of the relative weight of the source contributions from the basin associated to the maximum and from the remaining basins in the cluster. The source function has also allowed for a classification of hydrogen bonds in terms of characteristic source contributions to the density at the H-bond critical point from the H involved in the H-bond, the H-donor D, and the H-acceptor A. The source contribution from the H appears as the most distinctive marker of the H-bond strength, being highly negative for isolated H-bonds, slightly negative for polarized assisted H-bonds, close to zero for resonance-assisted H-bonds, and largely positive for charge-assisted H-bonds. The contributions from atoms other than H, D, and A strongly increase with decreasing H-bond strength, consistently with the parallel increased electrostatic character of the interaction. The correspondence between the classification provided by the Electron Localization Function topologic approach and by the source function has been highlighted. It is concluded that the source function represents a practical tool to disclose the local and nonlocal character of the electron density distributions and to quantify such a locality and nonlocality in terms of a physically sound and appealing chemical partitioning.