Relativistic LCAO with Minimax Principle and New Balanced Basis Sets

Relativistic density functional theory is widely applied in molecular calculations with heavy atoms, where relativistic and correlation effects are on the same footing. Variational stability of the Dirac Hamiltonian is a very important field of research from the beginning of relativistic molecular calculations on, among efforts for accuracy, efficiency, and density functional formulation, etc. Approximations of one- or two-component methods and searching for suitable basis sets are two major means for good projection power against the negative continuum. The minimax two-component spinor linear combination of atomic orbitals (LCAO) is applied in the present work for both light and super-heavy one-electron systems, providing good approximations in the whole energy spectrum, being close to the benchmark minimax finite element method (FEM) values and without spurious and contaminated states, in contrast to the presence of these artifacts in the traditional four-component spinor LCAO. The variational stability assures that minimax LCAO is bounded from below. New balanced basis sets, kinetic and potential defect balanced (TVDB), following the minimax idea, are applied with the Dirac Hamiltonian. Its performance in the same super-heavy one-electron quasi-molecules shows also very good projection capability against variational collapse, as the minimax LCAO is taken as the best projection to compare with. The TVDB method has twice as many basis coefficients as four-component spinor LCAO, which becomes now linear and overcomes the disadvantage of great time-consumption in the minimax method. The calculation with both the TVDB method and the traditional LCAO method for the dimers with elements in group 11 of the periodic table investigates their difference. New bigger basis sets are constructed than in previous research, achieving high accuracy within the functionals involved. Their difference in total energy is much smaller than the basis incompleteness error, showing that the traditional four-spinor LCAO keeps enough projection power from the numerical atomic orbitals and is suitable in research on relativistic quantum chemistry. In scattering investigations for the same comparison purpose, the failure of the traditional LCAO method of providing a stable spectrum with increasing size of basis sets is contrasted to the TVDB method, which contains no spurious states already without pre-orthogonalization of basis sets. Keeping the same conditions including the accuracy of matrix elements shows that the variational instability prevails over the linear dependence of the basis sets. The success of the TVDB method manifests its capability not only in relativistic quantum chemistry but also for scattering and under the influence of strong external electronic and magnetic fields. The good accuracy in total energy with large basis sets and the good projection property encourage wider research on different molecules, with better functionals, and on small effects.

Many-body theory of laser-induced ultrafast demagnetization and angular momentum transfer in ferromagnetic transition metals

An electronic theory is developed, which describes the ultrafast demagnetization in itinerant ferromagnets following the absorption of a femtosecond laser pulse. The present work intends to elucidate the microscopic physics of this ultrafast phenomenon by identifying its fundamental mechanisms. In particular, it aims to reveal the nature of the involved spin excitations and angular-momentum transfer between spin and lattice, which are still subjects of intensive debate. In the first preliminary part of the thesis the initial stage of the laser-induced demagnetization process is considered. In this stage the electronic system is highly excited by spin-conserving elementary excitations involved in the laser-pulse absorption, while the spin or magnon degrees of freedom remain very weakly excited. The role of electron-hole excitations on the stability of the magnetic order of one- and two-dimensional 3d transition metals (TMs) is investigated by using ab initio density-functional theory. The results show that the local magnetic moments are remarkably stable even at very high levels of local energy density and, therefore, indicate that these moments preserve their identity throughout the entire demagnetization process. In the second main part of the thesis a many-body theory is proposed, which takes into account these local magnetic moments and the local character of the involved spin excitations such as spin fluctuations from the very beginning. In this approach the relevant valence 3d and 4p electrons are described in terms of a multiband model Hamiltonian which includes Coulomb interactions, interatomic hybridizations, spin-orbit interactions, as well as the coupling to the time-dependent laser field on the same footing. An exact numerical time evolution is performed for small ferromagnetic TM clusters. The dynamical simulations show that after ultra-short laser pulse absorption the magnetization of these clusters decreases on a time scale of hundred femtoseconds. In particular, the results reproduce the experimentally observed laser-induced demagnetization in ferromagnets and demonstrate that this effect can be explained in terms of the following purely electronic non-adiabatic mechanism: First, on a time scale of 10–100 fs after laser excitation the spin-orbit coupling yields local angular-momentum transfer between the spins and the electron orbits, while subsequently the orbital angular momentum is very rapidly quenched in the lattice on the time scale of one femtosecond due to interatomic electron hoppings. In combination, these two processes result in a demagnetization within hundred or a few hundred femtoseconds after laser-pulse absorption.

A quantum molecular similarity analysis of changes in molecular electron density caused by basis set flotation and electric field application

Quantum molecular similarity (QMS) techniques are used to assess the response of the electron density of various small molecules to application of a static, uniform electric field. Likewise, QMS is used to analyze the changes in electron density generated by the process of floating a basis set. The results obtained show an interrelation between the floating process, the optimum geometry, and the presence of an external field. Cases involving the Le Chatelier principle are discussed, and an insight on the changes of bond critical point properties, self-similarity values and density differences is performed

Exploring chromium (VI) dioxodihalides chemistry: is density functional theory the most suitable tool?

A comparative systematic study of the CrO2F2 compound has been performed using different conventional ab initio methodologies and density functional procedures. Two points have been analyzed: first, the accuracy of results yielded by each method under study, and second, the computational cost required to reach such results. Weighing up both aspects, density functional theory has been found to be more appropriate than the Hartree-Fock (HF) and the analyzed post-HF methods. Hence, the structural characterization and spectroscopic elucidation of the full CrO2X2 series (X=F,Cl,Br,I) has been done at this level of theory. Emphasis has been given to the unknown CrO2I2 species, and specially to the UV/visible spectra of all four compounds. Furthermore, a topological analysis in terms of charge density distributions has revealed why the valence shell electron pair repulsion model fails in predicting the molecular shape of such CrO2X2 complexes

Estructura computacional i aplicacions de la semblança molecular quàntica

La tesis tracta diferents aspectes relacionats amb el càlcul de la semblança quàntica, així com la seva aplicació en la racionalització i predicció de l'activitat de fàrmacs. Es poden destacar dos progressos importants en el desenvolupament de noves metodologies que faciliten el càlcul de les mesures de semblança quàntica. En primer lloc, la descripció de les molècules mitjançant les funciones densitat aproximades PASA (Promolecular Atomic Shell Approximation) ha permès descriure amb suficient precisió la densitat electrònica dels sistemes moleculars analitzats, reduint substancialment el temps de càlcul de les mesures de semblança. En segon lloc, el desenvolupament de tècniques de superposició molecular específiques de les mesures de semblança quàntica ha permès resoldre el problema de l'alineament en l'espai dels compostos comparats. El perfeccionament d'aquests nous procediments i algoritmes matemàtics associats a les mesures de semblança molecular quàntica, ha estat essencial per poder progressar en diferents disciplines de la química computacional, sobretot les relacionades amb les anàlisis quantitatives entre les estructures moleculars i les seves activitats biològiques, conegudes amb les sigles angleses QSAR (Quantitative Structure-Activity Relationships). Precisament en l'àrea de les relacions estructura-activitat s'han presentat dues aproximacions fonamentades en la semblança molecular quàntica que s'originen a partir de dues representacions diferents de les molècules. La primera descripció considera la densitat electrònica global de les molècules i és important, entre altres, la disposició dels objectes comparats en l'espai i la seva conformació tridimensional. El resultat és una matriu de semblança amb les mesures de semblança de tots els parells de compostos que formen el conjunt estudiat. La segona descripció es fonamenta en la partició de la densitat global de les molècules en fragments. S'utilitzen mesures d'autosemblança per analitzar els requeriments bàsics d'una determinada activitat des del punt de vista de la semblança quàntica. El procés permet la detecció de les regions moleculars que són responsables d'una alta resposta biològica. Això permet obtenir un patró amb les regions actives que és d'evident interès per als propòsits del disseny de fàrmacs. En definitiva, s'ha comprovat que mitjançant la simulació i manipulació informàtica de les molècules en tres dimensions es pot obtenir una informació essencial en l'estudi de la interacció entre els fàrmacs i els seus receptors macromoleculars.

Electronic transitions and bonding properties in a series of five-coordinate “16-electron” complexes [Mn(CO)3(L2)]− (L2=chelating redox-active π-donor ligand)

In this article we present for the first time accurate density functional theory (DFT) and time-dependent (TD) DFT data for a series of electronically unsaturated five-coordinate complexes [Mn(CO)(3)(L-2)](-), where L-2 stands for a chelating strong pi-donor ligand represented by catecholate, dithiolate, amidothiolate, reduced alpha-diimine (1,4-dialkyl-1,4-diazabutadiene (R-DAB), 2,2'-bipyridine) and reduced 2,2'-biphosphinine types. The single-crystal X-ray structure of the unusual compound [Na(BPY)][Mn(CO)(3)(BPY)]center dot Et2O and the electronic absorption spectrum of the anion [Mn(CO)(3)(BPY)](-) are new in the literature. The nature of the bidentate ligand determines the bonding in the complexes, which varies between two limiting forms: from completely pi-delocalized diamagnetic {(CO)(3)Mn-L-2}(-) for L-2 = alpha-diimine or biphosphinine, to largely valence-trapped {(CO)(3)Mn-1-L-2(2-)}(-) for L-2(2-) = catecholate, where the formal oxidation states of Mn and L-2 can be assigned. The variable degree of the pi-delocalization in the Mn(L-2) chelate ring is indicated by experimental resonance Raman spectra of [Mn(CO)(3)(L-2)](-) (L-2=3,5-di-tBu-catecholate and iPr-DAB), where accurate assignments of the diagnostically important Raman bands have been aided by vibrational analysis. The L-2 = catecholate type of complexes is known to react with Lewis bases (CO substitution, formation of six-coordinate adducts) while the strongly pi-delocalized complexes are inert. The five-coordinate complexes adopt usually a distorted square pyramidal geometry in the solid state, even though transitions to a trigonal bipyramid are also not rare. The experimental structural data and the corresponding DFT-computed values of bond lengths and angles are in a very good agreement. TD-DFT calculations of electronic absorption spectra of the studied Mn complexes and the strongly pi-delocalized reference compound [Fe(CO)(3)(Me-DAB)] have reproduced qualitatively well the experimental spectra. Analyses of the computed electronic transitions in the visible spectroscopic region show that the lowest-energy absorption band always contains a dominant (in some cases almost exclusive) contribution from a pi(HOMO) -> pi*(LUMO) transition within the MnL2 metallacycle. The character of this optical excitation depends strongly on the composition of the frontier orbitals, varying from a partial L-2 -> Mn charge transfer (LMCT) through a fully delocalized pi(MnL2) -> pi*(MnL2) situation to a mixed (CO)Mn -> L-2 charge transfer (LLCT/MLCT). The latter character is most apparent in the case of the reference complex [Fe(CO)(3)(Me-DAB)]. The higher-lying, usually strongly mixed electronic transitions in the visible absorption region originate in the three lower-lying occupied orbitals, HOMO - 1 to HOMO - 3, with significant metal-d contributions. Assignment of these optical excitations to electronic transitions of a specific type is difficult. A partial LLCT/MLCT character is encountered most frequently. The electronic absorption spectra become more complex when the chelating ligand L-2, such as 2,2'-bipyridine, features two or more closely spaced low-lying empty pi* orbitals.

Wave-activity conservation laws and stability theorems for semi-geostrophic dynamics: part 1. pseudomomentum-based theory

There exists a well-developed body of theory based on quasi-geostrophic (QG) dynamics that is central to our present understanding of large-scale atmospheric and oceanic dynamics. An important question is the extent to which this body of theory may generalize to more accurate dynamical models. As a first step in this process, we here generalize a set of theoretical results, concerning the evolution of disturbances to prescribed basic states, to semi-geostrophic (SG) dynamics. SG dynamics, like QG dynamics, is a Hamiltonian balanced model whose evolution is described by the material conservation of potential vorticity, together with an invertibility principle relating the potential vorticity to the advecting fields. SG dynamics has features that make it a good prototype for balanced models that are more accurate than QG dynamics. In the first part of this two-part study, we derive a pseudomomentum invariant for the SG equations, and use it to obtain: (i) linear and nonlinear generalized Charney–Stern theorems for disturbances to parallel flows; (ii) a finite-amplitude local conservation law for the invariant, obeying the group-velocity property in the WKB limit; and (iii) a wave-mean-flow interaction theorem consisting of generalized Eliassen–Palm flux diagnostics, an elliptic equation for the stream-function tendency, and a non-acceleration theorem. All these results are analogous to their QG forms. The pseudomomentum invariant – a conserved second-order disturbance quantity that is associated with zonal symmetry – is constructed using a variational principle in a similar manner to the QG calculations. Such an approach is possible when the equations of motion under the geostrophic momentum approximation are transformed to isentropic and geostrophic coordinates, in which the ageostrophic advection terms are no longer explicit. Symmetry-related wave-activity invariants such as the pseudomomentum then arise naturally from the Hamiltonian structure of the SG equations. We avoid use of the so-called ‘massless layer’ approach to the modelling of isentropic gradients at the lower boundary, preferring instead to incorporate explicitly those boundary contributions into the wave-activity and stability results. This makes the analogy with QG dynamics most transparent. This paper treats the f-plane Boussinesq form of SG dynamics, and its recent extension to β-plane, compressible flow by Magnusdottir & Schubert. In the limit of small Rossby number, the results reduce to their respective QG forms. Novel features particular to SG dynamics include apparently unnoticed lateral boundary stability criteria in (i), and the necessity of including additional zonal-mean eddy correlation terms besides the zonal-mean potential vorticity fluxes in the wave-mean-flow balance in (iii). In the companion paper, wave-activity conservation laws and stability theorems based on the SG form of the pseudoenergy are presented.

Direct probing of photoinduced electron transfer in a self- assembled biomimetic [2Fe2S]-hydrogenase complex using ultrafast vibrational spectroscopy

A pyridyl-functionalized diiron dithiolate complex, [μ-(4-pyCH2−NMI-S2)Fe2(CO)6] (3, py = pyridine(ligand), NMI = naphthalene monoimide) was synthesized and fully characterized. In the presence of zinc tetraphenylporphyrin (ZnTPP), a self-assembled 3·ZnTPP complex was readily formed in CH2Cl2 by the coordination of the pyridyl nitrogen to the porphyrin zinc center. Ultrafast photoinduced electron transfer from excited ZnTPP to complex 3 in the supramolecular assembly was observed in real time by monitoring the ν(CO) and ν(CO)NMI spectral changes with femtosecond time-resolved infrared (TRIR) spectroscopy. We have confirmed that photoinduced charge separation produced the monoreduced species by comparing the time-resolved IR spectra with the conventional IR spectra of 3•− generated by reversible electrochemical reduction. The lifetimes for the charge separation and charge recombination processes were found to be τCS = 40 ± 3 ps and τCR = 205 ± 14 ps, respectively. The charge recombination is much slower than that in an analogous covalent complex, demonstrating the potential of a supramolecular approach to extend the lifetime of the chargeseparated state in photocatalytic complexes. The observed vibrational frequency shifts provide a very sensitive probe of the delocalization of the electron-spin density over the different parts of the Fe2S2 complex. The TR and spectro-electrochemical IR spectra, electron paramagnetic resonance spectra, and density functional theory calculations all show that the spin density in 3•− is delocalized over the diiron core and the NMI bridge. This delocalization explains why the complex exhibits low catalytic dihydrogen production even though it features a very efficient photoinduced electron transfer. The ultrafast porphyrin-to-NMIS2−Fe2(CO)6 photoinduced electron transfer is the first reported example of a supramolecular Fe2S2-hydrogenase model studied by femtosecond TRIR spectroscopy. Our results show that TRIR spectroscopy is a powerful tool to investigate photoinduced electron transfer in potential dihydrogen-producing catalytic complexes, and that way to optimize their performance by rational approaches.

Furan interaction with the Si(001)-(2 x 2) surface: structural, energetics, and vibrational spectra from first-principles

In this work we employ the state of the art pseudopotential method, within a generalized gradient approximation to the density functional theory, to investigate the adsorption process of furan on the silicon (001) surface. A direct comparison of different adsorption structures with x-ray photoelectron spectroscopy (XPS), ultra-violet photoelectron spectroscopy (UPS), high resolution electron energy loss spectroscopy (HREELS), near edge x-ray absorption fine structure (NEXAFS), and high resolution spectroscopy experimental data allows us to identify the [4 + 2] cycloaddition reaction as the most probable adsorbate. In addition, theoretical scanning tunnelling microscopy (STM) images are presented, with a view to contributing to further experimental investigations.

Oxygen-mediated electron transport through hybrid silicon-organic interfaces

We investigate from first principles the electronic and transport properties of hybrid organic/silicon interfaces of relevance to molecular electronics. We focus on conjugated molecules bonded to hydrogenated Si through hydroxyl or thiol groups. The electronic structure of the systems is addressed within density functional theory, and the electron transport across the interface is directly evaluated within the Landauer approach. The microscopic effects of molecule-substrate bonding on the transport efficiency are explicitly analyzed, and the oxygen-bonded interface is identified as a candidate system when preferential hole transfer is needed.

Structural parameters of polyacrylamide hydrogels obtained by the Equilibrium Swelling Theory

In this paper, the synthesis and structural characterization of a series of polyacrylamide hydrogels with different degrees of reticulation are reported. Although the Equilibrium Swelling Theory was recognized as a simple and reliable tool for the determination of structural hydrogels network parameters like equilibrium degree of swelling, cross-link ratio and mesh size, this is the first application of this methodology for polyacrylamide hydrogels. By changing the total monomer content in the synthesis solution (%T) from 5 to 30%, at a fixed value of cross-linker content in the total monomer amount (%C) of 5%, the final parameter obtained, the mesh size, can be tuned from 2 to 0.3 nm. It was also possible to change the mesh size (0.19-0.35) by varying %C from 5 to 12% (at %T = 20%). Scanning Electron Microscopy images for the most different formulations are shown and corroborate data obtained from the theory. (c) 2008 Elsevier Ltd. All rights reserved.

An existence theorem for an analytic first order PDE in the framework of Colombeau's theory

We study the existence of a holomorphic generalized solution u of the PDE[GRAPHICS]where f is a given holomorphic generalized function and (alpha (1),...alpha (m)) is an element of C-m\{0}.

Statistics and kinetics of single-molecule electron transfer dynamics in complex environments: A simulation model study

Fundação de Amparo à Pesquisa do Estado de São Paulo (FAPESP)

Effect of electron correlation on positronium formation in positron-helium scattering

A three-parameter correlated wave function for the helium ground state is used to study the scattering reaction e(+) + He --> He+ + Ps, where Ps stands for positronium atom. An exact analytical expression is constructed for the first Born scattering amplitude for Ps formation from helium. Based on this numerical results are presented for both differential and total cross-sections. It is demonstrated that the inner electronic correlation of the target atom plays a crucial role in explaining the discrepency between theory and experiment.