Cyclic variability of the circumstellar disk of the Be star zeta Tauri I. Long-term monitoring observations

Context. Emission lines formed in decretion disks of Be stars often undergo long-term cyclic variations, especially in the violet-to-red (V/R) ratio of their primary components. The underlying structural and dynamical variations of the disks are only partly understood. From observations of the bright Be-shell star. Tau, the possibly broadest and longest data set illustrating the prototype of this behaviour was compiled from our own and archival observations. It comprises optical and infrared spectra, broad-band polarimetry, and interferometric observations. Aims. The dense, long-time monitoring permits a better separation of repetitive and ephemeral variations. The broad wavelength coverage includes lines formed under different physical conditions, i.e. different locations in the disk, so that the dynamics can be probed throughout much of the disk. Polarimetry and interferometry constrain the spatial structure. All together, the objective is a better understand the dynamics and life cycle of decretion disks. Methods. Standard methods of data acquisition, reduction, and analysis were applied. Results. From 3 V/R cycles between 1997 and 2008, a mean cycle length in Ha of 1400-1430 days was derived. After each minimum in V/R, the shell absorption weakens and splits into two components, leading to 3 emission peaks. This phase may make the strongest contribution to the variability in cycle length. There is no obvious connection between the V/R cycle and the 133-day orbital period of the not otherwise detected companion. V/R curves of different lines are shifted in phase. Lines formed on average closer to the central star are ahead of the others. The shell absorption lines fall into 2 categories differing in line width, ionization/excitation potential, and variability of the equivalent width. They seem to form in separate regions of the disk, probably crossing the line of sight at different times. The interferometry has resolved the continuum and the line emission in Br gamma and HeI 2.06. The phasing of the Br gamma emission shows that the photocenter of the line-emitting region lies within the plane of the disk but is offset from the continuum source. The plane of the disk is constant throughout the observed V/R cycles. The observations lay the foundation for the fully self-consistent, one-armed, disk-oscillation model developed in Paper II.

Near infrared polarimetry of a sample of YSOs

Aims. Our goal is to study the physical properties of the circumstellar environment of young stellar objetcs (YSOs). In particular, the determination of the scattering mechanism can help us to constrain the optical depth of the disk and/or envelope in the near infrared. Methods. We used the IAGPOL imaging polarimeter along with the CamIV infrared camera at the LNA observatory to obtain near infrared polarimetry measurements in the H band of a sample of optically visible YSOs, namely, eleven T Tauri stars and eight Herbig Ae/Be stars. An independent determination of the disk (or jet) orientation was obtained for twelve objects from the literature. The circumstellar optical depth could then be estimated by comparing the integrated polarization position angle (PA) with the direction of the major axis of the disk projected onto the plane of the sky. Optically thin disks have, in general, a polarization PA that is perpendicular to the disk plane. In contrast, optically thick disks have polarization PAs parallel to the disks. Results. Among the T Tauri stars, three are consistent with having optically thin disks (AS 353A, RY Tau and UY Aur) and five with optically thick disks (V536 Aql, DG Tau, DO Tau, HL Tau and LkH alpha 358). Among the Herbig Ae/Be stars, two stars exhibit evidence of optically thin disks (Hen 3-1191 and VV Ser) and two of optically thick disks (PDS 453 and MWC 297). Our results seem consistent with optically thick disks at near infrared bands, which are more likely to be associated with younger YSOs. Marginal evidence of polarization reversal is found in RY Tau, RY Ori, WW Vul, and UY Aur. In the first three cases, this feature can be associated with the UXOR phenomenon. Correlations with the IRAS colors and the spectral index yielded evidence of an evolutionary segregation in which the disks tend to be optically thin when they are older.

H alpha spectropolarimetry of RY Tauri and PX Vulpeculae (Research Note)

Aims. To detect line effects using spectropolarimetry in order to find evidence of rotating disks and their respective symmetry axes in T Tauri stars. Methods. We used the IAGPOL imaging polarimeter along with the Eucalyptus-IFU to obtain spectropolarimetric measurements of the T Tauri stars RY Tau (two epochs) and PX Vul (one epoch). Evidence of line effects showing a loop in the Q-U diagram favors a compact rather than an extended source for the line photons in a rotating disk. In addition, the polarization position angle (PA) obtained using the line effect can constrain the symmetry axis of the disk. Results. RY Tau shows a variable H alpha double peak in 2004-2005 data. A polarization line effect is evident in the Q-U diagram for both epochs confirming a clockwise rotating disk. A single loop is evident in 2004 changing to a linear excursion plus a loop in 2005. Interestingly, the intrinsic PA calculated using the line effect is consistent between our two epochs (similar to 167 degrees). An alternative intrinsic PA computed from the interstellar polarization-corrected continuum and averaged between 2001-2005 yielded a PA similar to 137 degrees. This last value is closer to perpendicular to the observed disk direction (similar to 25 degrees), as expected from single scattering in an optically thin disk. For PX Vul, we detected spectral variability in H alpha along with non-variable continuum polarization when compared with previous data. The Q-U diagram shows a well-defined loop in H alpha associated with a counter-clockwise rotating disk. The symmetry axis inferred from the line effect has a PA similar to 91 degrees (with an ambiguity of 90 degrees). Our results confirm previous evidence that the emission line in T Tauri stars has its origin in a compact source scattered off a rotating accretion disk.

Remote sensing the vertical profile of cloud droplet effective radius, thermodynamic phase, and temperature

Cloud-aerosol interaction is a key issue in the climate system, affecting the water cycle, the weather, and the total energy balance including the spatial and temporal distribution of latent heat release. Information on the vertical distribution of cloud droplet microphysics and thermodynamic phase as a function of temperature or height, can be correlated with details of the aerosol field to provide insight on how these particles are affecting cloud properties and their consequences to cloud lifetime, precipitation, water cycle, and general energy balance. Unfortunately, today's experimental methods still lack the observational tools that can characterize the true evolution of the cloud microphysical, spatial and temporal structure in the cloud droplet scale, and then link these characteristics to environmental factors and properties of the cloud condensation nuclei. Here we propose and demonstrate a new experimental approach (the cloud scanner instrument) that provides the microphysical information missed in current experiments and remote sensing options. Cloud scanner measurements can be performed from aircraft, ground, or satellite by scanning the side of the clouds from the base to the top, providing us with the unique opportunity of obtaining snapshots of the cloud droplet microphysical and thermodynamic states as a function of height and brightness temperature in clouds at several development stages. The brightness temperature profile of the cloud side can be directly associated with the thermodynamic phase of the droplets to provide information on the glaciation temperature as a function of different ambient conditions, aerosol concentration, and type. An aircraft prototype of the cloud scanner was built and flew in a field campaign in Brazil. The CLAIM-3D (3-Dimensional Cloud Aerosol Interaction Mission) satellite concept proposed here combines several techniques to simultaneously measure the vertical profile of cloud microphysics, thermodynamic phase, brightness temperature, and aerosol amount and type in the neighborhood of the clouds. The wide wavelength range, and the use of multi-angle polarization measurements proposed for this mission allow us to estimate the availability and characteristics of aerosol particles acting as cloud condensation nuclei, and their effects on the cloud microphysical structure. These results can provide unprecedented details on the response of cloud droplet microphysics to natural and anthropogenic aerosols in the size scale where the interaction really happens.

Theory of small-signal ac response of a dielectric liquid containing two groups of ions

The analysis of Macdonald for electrolytes is generalized to the case in which two groups of ions are present. We assume that the electrolyte can be considered as a dispersion of ions in a dielectric liquid, and that the ionic recombination can be neglected. We present the differential equations governing the ionic redistribution when the liquid is subjected to an external electric field, describing the simultaneous diffusion of the two groups of ions in the presence of their own space charge fields. We investigate the influence of the ions on the impedance spectroscopy of an electrolytic cell. In the analysis, we assume that each group of ions have equal mobility, the electrodes perfectly block and that the adsorption phenomena can be neglected. In this framework, it is shown that the real part of the electrical impedance of the cell has a frequency dependence presenting two plateaux, related to a type of ambipolar and free diffusion coefficients. The importance of the considered problem on the ionic characterization performed by means of the impedance spectroscopy technique was discussed. (c) 2008 American Institute of Physics.

Numerical investigation of the quantum fluctuations of optical fields transmitted through an atomic medium

We have numerically solved the Heisenberg-Langevin equations describing the propagation of quantized fields through an optically thick sample of atoms. Two orthogonal polarization components are considered for the field, and the complete Zeeman sublevel structure of the atomic transition is taken into account. Quantum fluctuations of atomic operators are included through appropriate Langevin forces. We have considered an incident field in a linearly polarized coherent state (driving field) and vacuum in the perpendicular polarization and calculated the noise spectra of the amplitude and phase quadratures of the output field for two orthogonal polarizations. We analyze different configurations depending on the total angular momentum of the ground and excited atomic states. We examine the generation of squeezing for the driving-field polarization component and vacuum squeezing of the orthogonal polarization. Entanglement of orthogonally polarized modes is predicted. Noise spectral features specific to (Zeeman) multilevel configurations are identified.

Study of the optical and magnetic properties of pyrimidine in water combining PCM and QM/MM methodologies

The solvent effects on the low-lying absorption spectrum and on the (15)N chemical shielding of pyrimidine in water are calculated using the combined and sequential Monte Carlo simulation and quantum mechanical calculations. Special attention is devoted to the solute polarization. This is included by an iterative procedure previously developed where the solute is electrostatically equilibrated with the solvent. In addition, we verify the simple yet unexplored alternative of combining the polarizable continuum model (PCM) and the hybrid QM/MM method. We use PCM to obtain the average solute polarization and include this in the MM part of the sequential QM/MM methodology, PCM-MM/QM. These procedures are compared and further used in the discrete and the explicit solvent models. The use of the PCM polarization implemented in the MM part seems to generate a very good description of the average solute polarization leading to very good results for the n-pi* excitation energy and the (15)N nuclear chemical shield of pyrimidine in aqueous environment. The best results obtained here using the solute pyrimidine surrounded by 28 explicit water molecules embedded in the electrostatic field of the remaining 472 molecules give the statistically converged values for the low lying n-pi* absorption transition in water of 36 900 +/- 100 (PCM polarization) and 36 950 +/- 100 cm(-1) (iterative polarization), in excellent agreement among one another and with the experimental value observed with a band maximum at 36 900 cm(-1). For the nuclear shielding (15)N the corresponding gas-water chemical shift obtained using the solute pyrimidine surrounded by 9 explicit water molecules embedded in the electrostatic field of the remaining 491 molecules give the statistically converged values of 24.4 +/- 0.8 and 28.5 +/- 0.8 ppm, compared with the inferred experimental value of 19 +/- 2 ppm. Considering the simplicity of the PCM over the iterative polarization this is an important aspect and the computational savings point to the possibility of dealing with larger solute molecules. This PCM-MM/QM approach reconciles the simplicity of the PCM model with the reliability of the combined QM/MM approaches.

Electron collisions with the CH(2)O-H(2)O complex

We report cross sections for elastic collisions of low-energy electrons with the CH(2)O-H(2)O complex. We employed the Schwinger multichannel method with pseudopotentials in the static-exchange and in the static-exchange-polarization approximations for energies from 0.1 to 20 eV. We considered four different hydrogen-bonded structures for the complex that were generated by classical Monte Carlo simulations. Our aim is to investigate the effect of the water molecule on the pi* shape resonance of formaldehyde. Previous studies reported a pi* shape resonance for CH(2)O at around 1 eV. The resonance positions of the complexes appear at lower energies in all cases due to the mutual polarization between the two molecules. This indicates that the presence of water may favor dissociation by electron impact and may lead to an important effect on strand breaking in wet DNA by electron impact.

Inclusive cross section and double helicity asymmetry for pi(0) production in p plus p collisions at root s=62.4 GeV

The PHENIX experiment presents results from the RHIC 2006 run with polarized p + p collisions at root s = 62.4 GeV, for inclusive pi(0) production at midrapidity. Unpolarized cross section results are measured for transverse momenta p(T) = 0.5 to 7 GeV/c. Next-to-leading order perturbative quantum chromodynamics calculations are compared with the data, and while the calculations are consistent with the measurements, next-to-leading logarithmic corrections improve the agreement. Double helicity asymmetries A(LL) are presented for p(T) = 1 to 4 GeV/c and probe the higher range of Bjorken x of the gluon (x(g)) with better statistical precision than our previous measurements at root s = 200 GeV. These measurements are sensitive to the gluon polarization in the proton for 0.06 < x(g) < 0.4.

Combined Monte Carlo and quantum mechanics study of the solvatochromism of phenol in water. The origin of the blue shift of the lowest pi-pi* transition

A combined and sequential use of Monte Carlo simulations and quantum mechanical calculations is made to analyze the spectral shift of the lowest pi-pi* transition of phenol in water. The solute polarization is included using electrostatic embedded calculations at the MP2/aug-cc-pVDZ level giving a dipole moment of 2.25 D, corresponding to an increase of 76% compared to the calculated gas-phase value. Using statistically uncorrelated configurations sampled from the MC simulation,first-principle size-extensive calculations are performed to obtain the solvatochromic shift. Analysis is then made of the origin of the blue shift. Results both at the optimized geometry and in room-temperature liquid water show that hydrogen bonds of water with phenol promote a red shift when phenol is the proton-donor and a blue shift when phenol is the proton-acceptor. In the case of the optimized clusters the calculated shifts are in very good agreement with results obtained from mass-selected free jet expansion experiments. In the liquid case the contribution of the solute-solvent hydrogen bonds partially cancels and the total shift obtained is dominated by the contribution of the outer solvent water molecules. Our best result, including both inner and outer water molecules, is 570 +/- 35 cm(-1), in very good agreement with the small experimental shift of 460 cm(-1) for the absorption maximum.

Event structure and double helicity asymmetry in jet production from polarized p plus p collisions at root s=200 GeV

We report on the event structure and double helicity asymmetry (A(LL)) of jet production in longitudinally polarized p + p collisions at root s = 200 GeV. Photons and charged particles were measured by the PHENIX experiment at midrapidity vertical bar eta vertical bar < 0.35 with the requirement of a high-momentum (> 2 GeV/c) photon in the event. Event structure, such as multiplicity, p(T) density and thrust in the PHENIX acceptance, were measured and compared with the results from the PYTHIA event generator and the GEANT detector simulation. The shape of jets and the underlying event were well reproduced at this collision energy. For the measurement of jet A(LL), photons and charged particles were clustered with a seed-cone algorithm to obtain the cluster pT sum (p(T)(reco)). The effect of detector response and the underlying events on p(T)(reco) was evaluated with the simulation. The production rate of reconstructed jets is satisfactorily reproduced with the next-to-leading-order and perturbative quantum chromodynamics jet production cross section. For 4< p(T)(reco) < 12 GeV/c with an average beam polarization of < P > = 49% we measured Lambda(LL) = -0.0014 +/- 0.0037(stat) at the lowest p(T)(reco) bin (4-5 GeV= c) and -0.0181 +/- 0.0282(stat) at the highest p(T)(reco) bin (10-12 GeV= c) with a beam polarization scale error of 9.4% and a pT scale error of 10%. Jets in the measured p(T)(reco) range arise primarily from hard-scattered gluons with momentum fraction 0: 02 < x < 0: 3 according to PYTHIA. The measured A(LL) is compared with predictions that assume various Delta G(x) distributions based on the Gluck-Reya-Stratmann-Vogelsang parameterization. The present result imposes the limit -a.1 < integral(0.3)(0.02) dx Delta G(x, mu(2) = GeV2) < 0.4 at 95% confidence level or integral(0.3)(0.002) dx Delta G(x, mu(2) = 1 GeV2) < 0.5 at 99% confidence level.

Electronic properties of liquid ammonia: A sequential molecular dynamics/quantum mechanics approach

The electronic properties of liquid ammonia are investigated by a sequential molecular dynamics/quantum mechanics approach. Quantum mechanics calculations for the liquid phase are based on a reparametrized hybrid exchange-correlation functional that reproduces the electronic properties of ammonia clusters [(NH(3))(n); n=1-5]. For these small clusters, electron binding energies based on Green's function or electron propagator theory, coupled cluster with single, double, and perturbative triple excitations, and density functional theory (DFT) are compared. Reparametrized DFT results for the dipole moment, electron binding energies, and electronic density of states of liquid ammonia are reported. The calculated average dipole moment of liquid ammonia (2.05 +/- 0.09 D) corresponds to an increase of 27% compared to the gas phase value and it is 0.23 D above a prediction based on a polarizable model of liquid ammonia [Deng , J. Chem. Phys. 100, 7590 (1994)]. Our estimate for the ionization potential of liquid ammonia is 9.74 +/- 0.73 eV, which is approximately 1.0 eV below the gas phase value for the isolated molecule. The theoretical vertical electron affinity of liquid ammonia is predicted as 0.16 +/- 0.22 eV, in good agreement with the experimental result for the location of the bottom of the conduction band (-V(0)=0.2 eV). Vertical ionization potentials and electron affinities correlate with the total dipole moment of ammonia aggregates. (c) 2008 American Institute of Physics.

CMB in a box: Causal structure and the Fourier-Bessel expansion

This paper makes two points. First, we show that the line-of-sight solution to cosmic microwave anisotropies in Fourier space, even though formally defined for arbitrarily large wavelengths, leads to position-space solutions which only depend on the sources of anisotropies inside the past light cone of the observer. This foretold manifestation of causality in position (real) space happens order by order in a series expansion in powers of the visibility gamma = e(-mu), where mu is the optical depth to Thomson scattering. We show that the contributions of order gamma(N) to the cosmic microwave background (CMB) anisotropies are regulated by spacetime window functions which have support only inside the past light cone of the point of observation. Second, we show that the Fourier-Bessel expansion of the physical fields (including the temperature and polarization momenta) is an alternative to the usual Fourier basis as a framework to compute the anisotropies. The viability of the Fourier-Bessel series for treating the CMB is a consequence of the fact that the visibility function becomes exponentially small at redshifts z >> 10(3), effectively cutting off the past light cone and introducing a finite radius inside which initial conditions can affect physical observables measured at our position (x) over right arrow = 0 and time t(0). Hence, for each multipole l there is a discrete tower of momenta k(il) (not a continuum) which can affect physical observables, with the smallest momenta being k(1l) similar to l. The Fourier-Bessel modes take into account precisely the information from the sources of anisotropies that propagates from the initial value surface to the point of observation-no more, no less. We also show that the physical observables (the temperature and polarization maps), and hence the angular power spectra, are unaffected by that choice of basis. This implies that the Fourier-Bessel expansion is the optimal scheme with which one can compute CMB anisotropies.

Magnetoelectric effects in an organometallic quantum magnet

Metal-organic materials constitute a new field in which to search for ferroelectricity and coupling between electricity and magnetism. We observe a magnetic field-induced change in the electric polarization, Delta P(H), that reaches 50 mu C/m(2) in single crystals of NiCl(2)-4SC(NH(2))(2) (DTN). DTN forms a tetragonal structure that breaks inversion symmetry with the electrically polar thiourea molecules [SC(NH(2))] all tilted in the same direction along the c axis. The field H induces canted antiferromagnetism of the Ni S = 1 spins between 2 and 12 T and our measurements show that the electric polarization increases monotonically in this range, saturating above 12 T. By modeling the microscopic origin of this magnetoelectric effect, we find that the leading contribution to Delta P comes from the change in the crystal electric field, with a smaller contribution from magnetic exchange striction. The finite value of Delta P induced by magnetostriction results from the polar nature of the thiourea molecules bonded to the Ni atoms, and it is amplified by the softness of these organic molecules.

Optical third-harmonic spectroscopy of the magnetic semiconductor EuTe

EuTe possesses the centrosymmetric crystal structure m3m of rocksalt type in which the second-harmonic generation is forbidden in electric dipole approximation but the third-harmonic generation (THG) is allowed. We studied the THG spectra of this material and observed several resonances in the vicinity of the band gap at 2.2-2.5 eV and at higher energies up to 4 eV, which are related to four-photon THG processes. The observed resonances are assigned to specific combinations of electronic transitions between the ground 4f(7) state at the top of the valence band and excited 4f(6)5d(1) states of Eu(2+) ions, which form the lowest energy conduction band. Temperature, magnetic field, and rotational anisotropy studies allowed us to distinguish crystallographic and magnetic-field-induced contributions to the THG. A strong modification of THG intensity for the 2.4 eV band and suppression of the THG for the 3.15 eV band was observed in applied magnetic field. Two main features of the THG spectra were assigned to 5d(t(2g)) and 5d(e(g)) subbands at 2.4 eV and 3.15 eV, respectively. A microscopic quantum-mechanical model of the THG response was developed and its conclusions are in qualitative agreement with the experimental results.