Cyclic variability of the circumstellar disk of the Be star zeta Tauri II. Testing the 2D global disk oscillation model

Context. About 2/3 of the Be stars present the so-called V/R variations, a phenomenon characterized by the quasi-cyclic variation in the ratio between the violet and red emission peaks of the HI emission lines. These variations are generally explained by global oscillations in the circumstellar disk forming a one-armed spiral density pattern that precesses around the star with a period of a few years. Aims. This paper presents self-consistent models of polarimetric, photometric, spectrophotometric, and interferometric observations of the classical Be star zeta Tauri. The primary goal is to conduct a critical quantitative test of the global oscillation scenario. Methods. Detailed three-dimensional, NLTE radiative transfer calculations were carried out using the radiative transfer code HDUST. The most up-to-date research on Be stars was used as input for the code in order to include a physically realistic description for the central star and the circumstellar disk. The model adopts a rotationally deformed, gravity darkened central star, surrounded by a disk whose unperturbed state is given by a steady-state viscous decretion disk model. It is further assumed that this disk is in vertical hydrostatic equilibrium. Results. By adopting a viscous decretion disk model for zeta Tauri and a rigorous solution of the radiative transfer, a very good fit of the time-average properties of the disk was obtained. This provides strong theoretical evidence that the viscous decretion disk model is the mechanism responsible for disk formation. The global oscillation model successfully fitted spatially resolved VLTI/AMBER observations and the temporal V/R variations in the H alpha and Br gamma lines. This result convincingly demonstrates that the oscillation pattern in the disk is a one-armed spiral. Possible model shortcomings, as well as suggestions for future improvements, are also discussed.

The merging cluster Abell 1758 revisited: multi-wavelength observations and numerical simulations

Context. Cluster properties can be more distinctly studied in pairs of clusters, where we expect the effects of interactions to be strong. Aims. We here discuss the properties of the double cluster Abell 1758 at a redshift z similar to 0.279. These clusters show strong evidence for merging. Methods. We analyse the optical properties of the North and South cluster of Abell 1758 based on deep imaging obtained with the Canada-France-Hawaii Telescope (CFHT) archive Megaprime/Megacam camera in the g' and r' bands, covering a total region of about 1.05 x 1.16 deg(2), or 16.1 x 17.6 Mpc(2). Our X-ray analysis is based on archive XMM-Newton images. Numerical simulations were performed using an N-body algorithm to treat the dark-matter component, a semi-analytical galaxy-formation model for the evolution of the galaxies and a grid-based hydrodynamic code with a parts per million (PPM) scheme for the dynamics of the intra-cluster medium. We computed galaxy luminosity functions (GLFs) and 2D temperature and metallicity maps of the X-ray gas, which we then compared to the results of our numerical simulations. Results. The GLFs of Abell 1758 North are well fit by Schechter functions in the g' and r' bands, but with a small excess of bright galaxies, particularly in the r' band; their faint-end slopes are similar in both bands. In contrast, the GLFs of Abell 1758 South are not well fit by Schechter functions: excesses of bright galaxies are seen in both bands; the faint-end of the GLF is not very well defined in g'. The GLF computed from our numerical simulations assuming a halo mass-luminosity relation agrees with those derived from the observations. From the X-ray analysis, the most striking features are structures in the metal distribution. We found two elongated regions of high metallicity in Abell 1758 North with two peaks towards the centre. In contrast, Abell 1758 South shows a deficit of metals in its central regions. Comparing observational results to those derived from numerical simulations, we could mimic the most prominent features present in the metallicity map and propose an explanation for the dynamical history of the cluster. We found in particular that in the metal-rich elongated regions of the North cluster, winds had been more efficient than ram-pressure stripping in transporting metal-enriched gas to the outskirts. Conclusions. We confirm the merging structure of the North and South clusters, both at optical and X-ray wavelengths.

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.

Numerical study of a model for nonequilibrium wetting

We revisit the scaling properties of a model for nonequilibrium wetting [Phys. Rev. Lett. 79, 2710 (1997)], correcting previous estimates of the critical exponents and providing a complete scaling scheme. Moreover, we investigate a special point in the phase diagram, where the model exhibits a roughening transition related to directed percolation. We argue that in the vicinity of this point evaporation from the middle of plateaus can be interpreted as an external field in the language of directed percolation. This analogy allows us to compute the crossover exponent and to predict the form of the phase transition line close to its terminal point.

Edge effects in bilayer graphene nanoribbons: Ab initio total-energy density functional theory calculations

We show that the ground state of zigzag bilayer graphene nanoribbons is nonmagnetic. It also possesses a finite gap, which has a nonmonotonic dependence with the width as a consequence of the competition between bulk and strongly attractive edge interactions. All results were obtained using ab initio total-energy density functional theory calculations with the inclusion of parametrized van der Waals interactions.

Spin filtering and disorder-induced magnetoresistance in carbon nanotubes: Ab initio calculations

Nitrogen-doped carbon nanotubes can provide reactive sites on the porphyrin-like defects. It is well known that many porphyrins have transition-metal atoms, and we have explored transition-metal atoms bonded to those porphyrin-like defects inN-doped carbon nanotubes. The electronic structure and transport are analyzed by means of a combination of density functional theory and recursive Green's function methods. The results determined the heme B-like defect (an iron atom bonded to four nitrogens) is the most stable and has a higher polarization current for a single defect. With randomly positioned heme B defects in nanotubes a few hundred nanometers long, the polarization reaches near 100%, meaning they are effective spin filters. A disorder-induced magnetoresistance effect is also observed in those long nanotubes, and values as high as 20 000% are calculated with nonmagnectic eletrodes.

Barrier-free substitutional doping of graphene sheets with boron atoms: Ab initio calculations

Using ab initio methods, we propose a simple and effective way to substitutionally dope graphene sheets with boron. The method consists of selectively exposing each side of the graphene sheet to different elements. We first expose one side of the membrane to boron while the other side is exposed to nitrogen. Proceeding this way, the B atoms will be spontaneously incorporated into the graphene membrane without any activation barrier. In a second step, the system should be exposed to a H-rich environment, which will remove the CN radical from the layer and form HCN, leading to a perfect substitutional doping.

Optical absorption and electronic band structure first-principles calculations of alpha-glycine crystals

Light absorption of alpha-glycine crystals grown by slow evaporation at room temperature was measured, indicating a 5.11 +/- 0.02 eV energy band gap. Structural, electronic, and optical absorption properties of alpha-glycine crystals were obtained by first-principles quantum mechanical calculations using density functional theory within the generalized gradient approximation in order to understand this result. To take into account the contribution of core electrons, ultrasoft and norm-conserving pseudopotentials, as well as an all electron approach were considered to compute the electronic density of states and band structure of alpha-glycine crystals. They exhibit three indirect energy band gaps and one direct Gamma-Gamma energy gap around 4.95 eV. The optical absorption related to transitions between the top of the valence band and the bottom of the conduction band involves O 2p valence states and C, O 2p conduction states, with the carboxyl group contributing significantly to the origin of the energy band gap. The calculated optical absorption is highly dependent on the polarization of the incident radiation due to the spatial arrangement of the dipolar glycine molecules; in the case of a polycrystalline sample, the first-principles calculated optical absorption is in good agreement with the measurement when a rigid energy shift is applied.

Understanding the ""Antikick"" in the Merger of Binary Black Holes

The generation of a large recoil velocity from the inspiral and merger of binary black holes represents one of the most exciting results of numerical-relativity calculations. While many aspects of this process have been investigated and explained, the ""antikick,"" namely, the sudden deceleration after the merger, has not yet found a simple explanation. We show that the antikick can be understood in terms of the radiation from a deformed black hole where the anisotropic curvature distribution on the horizon correlates with the direction and intensity of the recoil. Our analysis is focused on Robinson-Trautman spacetimes and allows us to measure both the energies and momenta radiated in a gauge-invariant manner. At the same time, this simpler setup provides the qualitative and quantitative features of merging black holes, opening the way to a deeper understanding of the nonlinear dynamics of black-hole spacetimes.

Numerical evidence against a conjecture on the cover time of planar graphs

We investigate a conjecture on the cover times of planar graphs by means of large Monte Carlo simulations. The conjecture states that the cover time tau (G(N)) of a planar graph G(N) of N vertices and maximal degree d is lower bounded by tau (G(N)) >= C(d)N(lnN)(2) with C(d) = (d/4 pi) tan(pi/d), with equality holding for some geometries. We tested this conjecture on the regular honeycomb (d = 3), regular square (d = 4), regular elongated triangular (d = 5), and regular triangular (d = 6) lattices, as well as on the nonregular Union Jack lattice (d(min) = 4, d(max) = 8). Indeed, the Monte Carlo data suggest that the rigorous lower bound may hold as an equality for most of these lattices, with an interesting issue in the case of the Union Jack lattice. The data for the honeycomb lattice, however, violate the bound with the conjectured constant. The empirical probability distribution function of the cover time for the square lattice is also briefly presented, since very little is known about cover time probability distribution functions in general.

Precise determination of quantum critical points by the violation of the entropic area law

Finite-size scaling analysis turns out to be a powerful tool to calculate the phase diagram as well as the critical properties of two-dimensional classical statistical mechanics models and quantum Hamiltonians in one dimension. The most used method to locate quantum critical points is the so-called crossing method, where the estimates are obtained by comparing the mass gaps of two distinct lattice sizes. The success of this method is due to its simplicity and the ability to provide accurate results even considering relatively small lattice sizes. In this paper, we introduce an estimator that locates quantum critical points by exploring the known distinct behavior of the entanglement entropy in critical and noncritical systems. As a benchmark test, we use this new estimator to locate the critical point of the quantum Ising chain and the critical line of the spin-1 Blume-Capel quantum chain. The tricritical point of this last model is also obtained. Comparison with the standard crossing method is also presented. The method we propose is simple to implement in practice, particularly in density matrix renormalization group calculations, and provides us, like the crossing method, amazingly accurate results for quite small lattice sizes. Our applications show that the proposed method has several advantages, as compared with the standard crossing method, and we believe it will become popular in future numerical studies.

Accurate band gaps of AlGaN, InGaN, and AlInN alloys calculations based on LDA-1/2 approach

We present parameter-free calculations of electronic properties of InGaN, InAlN, and AlGaN alloys. The calculations are based on a generalized quasichemical approach, to account for disorder and composition effects, and first-principles calculations within the density functional theory with the LDA-1/2 approach, to accurately determine the band gaps. We provide precise results for AlGaN, InGaN, and AlInN band gaps for the entire range of compositions, and their respective bowing parameters. (C) 2011 American Institute of Physics. [doi:10.1063/1.3576570]

Structure, stability, depolarized light scattering, and vibrational spectra of fullerenols from all-electron density-functional-theory calculations

We have investigated the stability, electronic properties, Rayleigh (elastic), and Raman (inelastic) depolarization ratios, infrared and Raman absorption vibrational spectra of fullerenols [C(60)(OH)(n)] with different degrees of hydroxylation by using all-electron density-functional-theory (DFT) methods. Stable arrangements of these molecules were found by means of full geometry optimizations using Becke's three-parameter exchange functional with the Lee, Yang, and Parr correlation functional. This DFT level has been combined with the 6-31G(d,p) Gaussian-type basis set, as a compromise between accuracy and capability to treat highly hydroxylated fullerenes, e.g., C(60)(OH)(36). Thus, the molecular properties of fullerenols were systematically analyzed for structures with n=1, 2, 3, 4, 8, 10, 16, 18, 24, 32, and 36. From the electronic structure analysis of these molecules, we have evidenced an important effect related to the weak chemical reactivity of a possible C(60)(OH)(24) isomer. To investigate Raman scattering and the vibrational spectra of the different fullerenols, frequency calculations are carried out within the harmonic approximation. In this case a systematic study is only performed for n=1-4, 8, 10, 16, 18, and 24. Our results give good agreements with the expected changes in the spectral absorptions due to the hydroxylation of fullerenes.

Bulk structures of PtO and PtO(2) from density functional calculations

Platinum plays an important role in catalysis and electrochemistry, and it is known that the direct interaction of oxygen with Pt surfaces can lead to the formation of platinum oxides (PtO(x)), which can affect the reactivity. To contribute to the atomistic understanding of the atomic structure of PtO(x), we report a density functional theory study of the atomic structure of bulk PtO(x) (1 <= x <= 2). From our calculations, we identified a lowest-energy structure (GeS type, space group Pnma) for PtO, which is 0.181 eV lower in energy than the structure suggested by W. J. Moore and L. Pauling [J. Am. Chem. Soc. 63, 1392 (1941)] (PtS type). Furthermore, two atomic structures were identified for PtO(2), which are almost degenerate in energy with the lowest-energy structure reported so far for PtO(2) (CaCl(2) type). Based on our results and analysis, we suggest that Pt and O atoms tend to form octahedron motifs in PtO(x) even at lower O composition by the formation of Pt-Pt bonds.

Failure analysis of low velocity impact on thin composite laminates: Experimental and numerical approaches

The dynamic behavior of composite laminates is very complex because there are many concurrent phenomena during composite laminate failure under impact load. Fiber breakage, delaminations, matrix cracking, plastic deformations due to contact and large displacements are some effects which should be considered when a structure made from composite material is impacted by a foreign object. Thus, an investigation of the low velocity impact on laminated composite thin disks of epoxy resin reinforced by carbon fiber is presented. The influence of stacking sequence and energy impact was investigated using load-time histories, displacement-time histories and energy-time histories as well as images from NDE. Indentation tests results were compared to dynamic results, verifying the inertia effects when thin composite laminate was impacted by foreign object with low velocity. Finite element analysis (FEA) was developed, using Hill`s model and material models implemented by UMAT (User Material Subroutine) into software ABAQUS (TM), in order to simulate the failure mechanisms under indentation tests. (C) 2007 Elsevier Ltd. All rights reserved.