Variational calculation of vibrational linear and nonlinear optical properties

A variational approach for reliably calculating vibrational linear and nonlinear optical properties of molecules with large electrical and/or mechanical anharmonicity is introduced. This approach utilizes a self-consistent solution of the vibrational Schrödinger equation for the complete field-dependent potential-energy surface and, then, adds higher-level vibrational correlation corrections as desired. An initial application is made to static properties for three molecules of widely varying anharmonicity using the lowest-level vibrational correlation treatment (i.e., vibrational Møller-Plesset perturbation theory). Our results indicate when the conventional Bishop-Kirtman perturbation method can be expected to break down and when high-level vibrational correlation methods are likely to be required. Future improvements and extensions are discussed

A hydrogen-1 nuclear magnetic resonance, mass spectral and extended Hückel study of some olefin complexes of rhodium(I) and iridium(I) and the crystal and molecular structure of η4-2,4-dimethylpenta-1,4-diene(η5-formylcyclopentadienyl)rhodium(I)

A 1H NMR study of monosubstituted η-cyclopentadienyl-rhodium(I) complexes of type LLRh(C5H4X) and -iridium(I) complexes of type L2Ir(C5H4X) (L = ethene, LL = 1,3- or 1,5-diolefin; X = C(C6H5)3, CHO, or COOCH3) has been carried out. For complexes of both metals in which the neutral ligand is ethene or a non-conjugated diolefin the NMR spectra of the cyclopentadienyl protons are unusual in that H(2), H(5) resonate to high field either at room temperature or below. The corresponding NMR spectra for the cyclopentadienyl ring protons of complexes where the neutral ligand is a conjugated diene are, with one exception, normal. A single crystal X-ray structural analysis of (η4-2,4-dimethylpenta-1,4-diene)(η5-formylcyclopentadienyl)rhodium(I) (which exhibits an abnormal 1H NMR spectrum) reveals substantial localisation of electron density in the C(3)C(4) Cp ring bond (1.283(33) Å) which may be consistent with a contribution from an ‘allyl-ene’ rotamer to the ring—metal bonding scheme. An extended Hückel calculation with self consistent charge iteration was performed on this complex. The results predict a greater Mulliken overlap population for the C(3)C(4) bond in the cyclopentadienyl ring and show that the localisation is dependent on both the Cp ring substituent and the nature of the diolefin. The mass spectral fragmentation patterns of some representative diene complexes of iridium(I) and rhodium(I) are presented.

Electron trapping in shear Alfvén waves that power the aurora

Results from 1D Vlasov drift-kinetic plasma simulations reveal how and where auroral electrons are accelerated along Earth’s geomagnetic field. In the warm plasma sheet, electrons become trapped in shear Alfven waves, preventing immediate wave damping. As waves move to regions with larger vTe=vA, their parallel electric field decreases, and the trapped electrons escape their influence. The resulting electron distribution functions compare favorably with in situ observations, demonstrating for the first time a self-consistent link between Alfven waves and electrons that form aurora.

Coupled Kelvin-wave and mirage-wave instabilities in semigeostrophic dynamics

A weak instability mode, associated with phase-locked counterpropagating coastal Kelvin waves in horizontal anticyclonic shear, is found in the semigeostrophic (SG) equations for stratified flow in a channel. This SG instability mode approximates a similar mode found in the Euler equations in the limit in which particle-trajectory slopes are much smaller than f/N, where f is the Coriolis frequency and N > f the buoyancy frequency. Though weak under normal parameter conditions, this instability mode is of theoretical interest because its existence accounts for the failure of an Arnol’d-type stability theorem for the SG equations. In the opposite limit, in which the particle motion is purely vertical, the Euler equations allow only buoyancy oscillations with no horizontal coupling. The SG equations, on the other hand, allow a physically spurious coastal “mirage wave,” so called because its velocity field vanishes despite a nonvanishing disturbance pressure field. Counterpropagating pairs of these waves can phase-lock to form a spurious “mirage-wave instability.” Closer examination shows that the mirage wave arises from failure of the SG approximations to be self-consistent for trajectory slopes f/N.

Electron acceleration and parallel electric fields due to kinetic Alfvén waves in plasma with similar thermal and Alfvén speeds

We investigate electron acceleration due to shear Alfven waves in a collissionless plasma for plasma parameters typical of 4–5RE radial distance from the Earth along auroral field lines. Recent observational work has motivated this study, which explores the plasma regime where the thermal velocity of the electrons is similar to the Alfven speed of the plasma, encouraging Landau resonance for electrons in the wave fields. We use a self-consistent kinetic simulation model to follow the evolution of the electrons as they interact with a short-duration wave pulse, which allows us to determine the parallel electric field of the shear Alfven wave due to both electron inertia and electron pressure effects. The simulation demonstrates that electrons can be accelerated to keV energies in a modest amplitude sub-second period wave. We compare the parallel electric field obtained from the simulation with those provided by fluid approximations.

Diagnostic tool to analyse colour-magnitude diagrams of poorly populated stellar concentrations

The dynamical processes that lead to open cluster disruption cause its mass to decrease. To investigate such processes from the observational point of view, it is important to identify open cluster remnants (OCRs), which are intrinsically poorly populated. Due to their nature, distinguishing them from field star fluctuations is still an unresolved issue. In this work, we developed a statistical diagnostic tool to distinguish poorly populated star concentrations from background field fluctuations. We use 2MASS photometry to explore one of the conditions required for a stellar group to be a physical group: to produce distinct sequences in a colour-magnitude diagram (CMD). We use automated tools to (i) derive the limiting radius; (ii) decontaminate the field and assign membership probabilities; (iii) fit isochrones; and (iv) compare object and field CMDs, considering the isochrone solution, in order to verify the similarity. If the object cannot be statistically considered as a field fluctuation, we derive its probable age, distance modulus, reddening and uncertainties in a self-consistent way. As a test, we apply the tool to open clusters and comparison fields. Finally, we study the OCR candidates DoDz 6, NGC 272, ESO 435 SC48 and ESO 325 SC15. The tool is optimized to treat these low-statistic objects and to separate the best OCR candidates for studies on kinematics and chemical composition. The study of the possible OCRs will certainly provide a deep understanding of OCR properties and constraints for theoretical models, including insights into the evolution of open clusters and dissolution rates.

THREE-DIMENSIONAL NUMERICAL SIMULATIONS OF MAGNETIZED WINDS OF SOLAR-LIKE STARS

By means of self-consistent three-dimensional magnetohydrodynamics (MHD) numerical simulations, we analyze magnetized solar-like stellar winds and their dependence on the plasma-beta parameter (the ratio between thermal and magnetic energy densities). This is the first study to perform such analysis solving the fully ideal three-dimensional MHD equations. We adopt in our simulations a heating parameter described by gamma, which is responsible for the thermal acceleration of the wind. We analyze winds with polar magnetic field intensities ranging from 1 to 20 G. We show that the wind structure presents characteristics that are similar to the solar coronal wind. The steady-state magnetic field topology for all cases is similar, presenting a configuration of helmet streamer-type, with zones of closed field lines and open field lines coexisting. Higher magnetic field intensities lead to faster and hotter winds. For the maximum magnetic intensity simulated of 20 G and solar coronal base density, the wind velocity reaches values of similar to 1000 km s(-1) at r similar to 20r(0) and a maximum temperature of similar to 6 x 10(6) K at r similar to 6r(0). The increase of the field intensity generates a larger ""dead zone"" in the wind, i.e., the closed loops that inhibit matter to escape from latitudes lower than similar to 45 degrees extend farther away from the star. The Lorentz force leads naturally to a latitude-dependent wind. We show that by increasing the density and maintaining B(0) = 20 G the system recover back to slower and cooler winds. For a fixed gamma, we show that the key parameter in determining the wind velocity profile is the beta-parameter at the coronal base. Therefore, there is a group of magnetized flows that would present the same terminal velocity despite its thermal and magnetic energy densities, as long as the plasma-beta parameter is the same. This degeneracy, however, can be removed if we compare other physical parameters of the wind, such as the mass-loss rate. We analyze the influence of gamma in our results and we show that it is also important in determining the wind structure.

THE VORONOI TESSELLATION CLUSTER FINDER IN 2+1 DIMENSIONS

We present a detailed description of the Voronoi Tessellation (VT) cluster finder algorithm in 2+1 dimensions, which improves on past implementations of this technique. The need for cluster finder algorithms able to produce reliable cluster catalogs up to redshift 1 or beyond and down to 10(13.5) solar masses is paramount especially in light of upcoming surveys aiming at cosmological constraints from galaxy cluster number counts. We build the VT in photometric redshift shells and use the two-point correlation function of the galaxies in the field to both determine the density threshold for detection of cluster candidates and to establish their significance. This allows us to detect clusters in a self-consistent way without any assumptions about their astrophysical properties. We apply the VT to mock catalogs which extend to redshift 1.4 reproducing the ACDM cosmology and the clustering properties observed in the Sloan Digital Sky Survey data. An objective estimate of the cluster selection function in terms of the completeness and purity as a function of mass and redshift is as important as having a reliable cluster finder. We measure these quantities by matching the VT cluster catalog with the mock truth table. We show that the VT can produce a cluster catalog with completeness and purity > 80% for the redshift range up to similar to 1 and mass range down to similar to 10(13.5) solar masses.

Fases magnéticas de sistemas nanoestruturados de terras-raras

There is presently a worldwide interest in artificial magnetic systems which guide research activities in universities and companies. Thin films and multilayers have a central role, revealing new magnetic phases which often lead to breakthroughs and new technology standards, never thought otherwise. Surface and confinement effects cause large impact in the magnetic phases of magnetic materials with bulk spatially periodic patterns. New magnetic phases are expected to form in thin film thicknesses comparable to the length of the intrinsic bulk magnetic unit cell. Helimagnetic materials are prototypes in this respect, since the bulk magnetic phases consist in periodic patterns with the length of the helical pitch. In this thesis we study the magnetic phases of thin rare-earth films, with surfaces oriented along the (002) direction. The thesis includes the investigation of the magnetic phases of thin Dy and Ho films, as well as the thermal hysteresis cycles of Dy thin films. The investigation of the thermal hysteresis cycles of thin Dy films has been done in collaboration with the Laboratory of Magnetic Materials of the University of Texas, at Arlington. The theoretical modeling is based on a self-consistent theory developed by the Group of Magnetism of UFRN. Contributions from the first and second neighbors exchange energy, from the anisotropy energy and the Zeeman energy are calculated in a set of nonequivalent magnetic ions, and the equilibrium magnetic phases, from the Curie temperature up to the Nèel temperature, are determined in a self-consistent manner, resulting in a vanishing torque in the magnetic ions at all planes across the thin film. Our results reproduce the known isothermal and iso-field curves of bulk Dy and Ho, and the known spin-slip phases of Ho, and indicate that: (i) the confinement in thin films leads to a new magnetic phase, with alternate helicity, which leads to the measured thermal hysteresis of Dy ultrathin films, with thicknesses ranging from 4 nm to 16 nm; (ii) thin Dy films have anisotropy dominated surface lock-in phases, with alignment of surface spins along the anisotropy easy axis directions, similar to the known spin-slip phases of Ho ( which form in the bulk and are commensurate to the crystal lattice); and (iii) the confinement in thin films change considerably the spin-slip patterns of Ho.

A theoretical analysis on the intramolecular proton transfer of alpha-alanine in an aqueous medium

Intramolecular proton transfer from oxygen to nitrogen atoms in the alpha-alanine amino acid has been studied by ab initio methods at the HF/6-31G*, HF/6-31 ++ G** and MP2/6-31 ++ G** levels of calculation including the solvent effects by means of self-consistent reaction field theory. An analysis of the results based on the natural bond orbital charges shows that the transition structure presents an imbalance in the sense that the charge shift lags behind the proton transfer and that the bond formation is always in advance with respect to the bond cleavage. All calculation levels show that the barrier height associated with the conformational change on alpha-alanine is larger than the proton transfer process. (C) 1998 Elsevier B.V. B.V. All rights reserved.

A theoretical study on cytosine tautomers in aqueous media by using continuum models

B3LYP/6-31 + + G** and MP2/6-31 + + G** calculations have been carried out to study six tautomers of the nucleic acid base cytosine in aqueous media. Solvent effects have been analyzed using the self-consistent reaction field theory with two continuum methods. Relative stabilities and optimized geometries have been calculated for the tautomers and compared with experimental data. The present results show the importance of electrostatic solvent effects in determining observable properties of the cytosine tautomers. The amino-oxo form (C1) is the most abundant tautomer in aqueous media while the other amino-oxo form (C4) is the most energetically favored when solvent effects are included. These results can be justified by the larger values of the dipole moments for both C1 and C4 tautomers. Theoretical and experimental results of the harmonic vibrational frequencies and rotational constants show good agreement. (C) 2000 Elsevier B.V. B.V. All rights reserved.

Density functional study of the 5-methylcytosine tautomers

B3LYP/6-31++G** calculations to study seven tautomers of 5-methylcytosine in aqueous media have been carried out. Optimized geometries and relative stabilities for the different tautomers have been calculated in the gas phase, including interaction with two discrete water molecules and taking into account the solvent effects by using the self-consistent reaction field theory. The role of specific and bulk contributions of solvent effect on the observable properties of the 5-methylcytosine is clarified. The amino-oxo form is the most abundant tautomer in aqueous media. A reaction pathway connecting amino-oxo and amino-hydroxy forms along the corresponding transition structures has been characterized. Good agreement between theoretical and available experimental results of harmonic vibration frequencies is found. (C) 2001 Elsevier B.V. B.V. All rights reserved.

Electromagnetic energy-momentum and forces in matter

We discuss the electromagnetic energy-momentum distribution and the mechanical forces of the electromagnetic field in material media. There is a long-standing controversy on these notions. The Minkowski and the Abraham energy-momentum tensors are the most well-known ones. We propose a solution of this problem which appears to be natural and self-consistent from both a theoretical and an experimental point of view. (C) 2003 Elsevier B.V. B.V. All rights reserved.

Rho-nucleon tensor coupling and charge-exchange resonances

The Gamow-Teller resonance in Pb-208 is discussed in the context of a self-consistent RPA, based on the relativistic mean field theory. We inquire on the possibility of substituting the phenomenological Landau-Migdal force by a microscopic nucleon-nucleon interaction, generated from the rho-nucleon tensor coupling. The effect of this coupling turns out to be very small when the short range correlations are not taken into account, but too large when these correlations are simulated by the simple extraction of the contact terms from the resulting nucleon-nucleon interaction. (C) 2000 Elsevier B.V. B.V. All rights reserved.

Quantum equivalence between the self-dual and the Maxwell-Chern-Simons models nonlinearly coupled to U(1) scalar fields

The use of master actions to prove duality at quantum level becomes cumbersome if one of the dual fields interacts nonlinearly with other fields. This is the case of the theory considered here consisting of U(1) scalar fields coupled to a self-dual field through a linear and a quadratic term in the self-dual field. Integrating perturbatively over the scalar fields and deriving effective actions for the self-dual and the gauge field we are able to consistently neglect awkward extra terms generated via master action and establish quantum duality up to cubic terms in the coupling constant. The duality holds for the partition function and some correlation functions. The absence of ghosts imposes restrictions on the coupling with the scalar fields.