Orientational order in dipolar fluids consisting of nonspherical hard particles

We investigate fluids of dipolar hard particles by a certain variant of density-functional theory. The proper treatment of the long range of the dipolar interactions yields a contribution to the free energy which favors ferromagnetic order. This corrects previous theoretical analyses. We determine phase diagrams for dipolar ellipsoids and spherocylinders as a function of the aspect ratio of the particles and their dipole moment. In the nonpolar limit the results for the phase boundary between the isotropic and nematic phase agree well with simulation data. Adding a longitudinal dipole moment favors the nematic phase. For oblate or slightly elongated particles we find a ferromagnetic liquid phase, which has also been detected in computer simulations of fluids consisting of spherical dipolar particles. The detailed structure of the phase diagram and its evolution upon changing the aspect ratio are discussed in detail.Comment: 35 pages LaTeX with epsf style, 11 figures in eps format, submitted to Phys. Rev.

Multi-Wavelength Implications of the Companion Star in Eta Carinae

Eta Carinae is considered to be a massive colliding wind binary system with a highly eccentric (e \sim 0.9), 5.54-yr orbit. However, the companion star continues to evade direct detection as the primary dwarfs its emission at most wavelengths. Using three-dimensional (3-D) SPH simulations of Eta Car's colliding winds and radiative transfer codes, we are able to compute synthetic observables across multiple wavebands for comparison to the observations. The models show that the presence of a companion star has a profound influence on the observed HST/STIS UV spectrum and H-alpha line profiles, as well as the ground-based photometric monitoring. Here, we focus on the Bore Hole effect, wherein the fast wind from the hot secondary star carves a cavity in the dense primary wind, allowing increased escape of radiation from the hotter/deeper layers of the primary's extended wind photosphere. The results have important implications for interpretations of Eta Car's observables at multiple wavelengths.Comment: 5 pages, 4 figures, To be published in the proceedings of the meeting 'Four Decades of Research on Massive Stars' in honor of Tony Moffat, 11-15 July 2011, Saint-Michel-des-Saints, Quebe

On the influence of the companion star in Eta Carinae: 2D radiative transfer modeling of the ultraviolet and optical spectra

We present 2D radiative transfer modeling of the Eta Carinae binary system accounting for the presence of a wind-wind collision (WWC) cavity carved in the optically-thick wind of the primary star. By comparing synthetic line profiles with HST/STIS spectra obtained near apastron, we show that the WWC cavity has a strong influence on multi-wavelength diagnostics. This influence is regulated by the modification of the optical depth in the continuum and spectral lines. We find that H-alpha, H-beta, and Fe II lines are the most affected by the WWC cavity, since they form over a large volume of the primary wind. These spectral lines depend on latitude and azimuth since, according to the orientation of the cavity, different velocity regions of a spectral line are affected. For 2D models with orientation corresponding to orbital inclination angle 110deg < i < 140deg and longitude of periastron 210deg < omega < 330deg, the blueshifted and zero-velocity regions of the line profiles are the most affected. These orbital orientations are required to simultaneously fit the UV and optical spectrum of Eta Car, for a half-opening angle of the cavity in the range 50-70deg. We find that the excess P-Cygni absorption seen in H-alpha, H-beta and optical Fe II lines in spherical models becomes much weaker or absent in the 2D models, in agreement with the observations. The observed UV spectrum of Eta Car, dominated by Fe II absorption lines, is superbly reproduced by our 2D cavity models. Small discrepancies still remain, as H-gamma and H-delta absorptions are overestimated by our models. We suggest that photoionization of the wind of the primary by the hot companion star is responsible for the weak absorption seen in these lines. Our CMFGEN models indicate that the primary star has a mass-loss rate of 8.5x10e-4 Msun/yr and wind terminal velocity of 420 km/s around the 2000 apastron.Comment: 20 pages, 14 figures, accepted for publication in MNRA

Inhomogeneous magnetization in dipolar ferromagnetic liquids

At high densities fluids of strongly dipolar spherical particles exhibit spontaneous long-ranged orientational order. Typically, due to demagnetization effects induced by the long range of the dipolar interactions, the magnetization structure is spatially inhomogeneous and depends on the shape of the sample. We determine this structure for a cubic sample by the free minimization of an appropriate microscopic density functional using simulated annealing. We find a vortex structure resembling four domains separated by four domain walls whose thickness increases proportional to the system size L. There are indications that for large L the whole configuration scales with the system size. Near the axis of the mainly planar vortex structure the direction of the magnetization escapes into the third dimension or, at higher temperatures, the absolute value of the magnetization is strongly reduced. Thus the orientational order is characterized by two point defects at the top and the bottom of the sample, respectively. The equilibrium structure in an external field and the transition to a homogeneous magnetization for strong fields are analyzed, too.Comment: 17 postscript figures included, submitted to Phys. Rev.

Spectroscopic and photometric oscillatory envelope variability during the S Doradus outburst of the Luminous Blue Variable R71

To better understand the LBV phenomenon, we analyze multi-epoch and multi-wavelength spectra and photometry of R71. Pre-outburst spectra are analyzed with the radiative transfer code CMFGEN to determine the star's fundamental stellar parameters. During quiescence, R71 has an effective temperature of $T_\mathrm{{eff}} = 15\,500~K$ and a luminosity of log$(L_*/L_{\odot})$ = 5.78 and is thus a classical LBV, but at the lower luminosity end of this group. We determine its mass-loss rate to $4.0 \times 10^{-6}~M_{\odot}~$yr$^{-1}$. We present R71's spectral energy distribution from the near-ultraviolet to the mid-infrared during its present outburst. Mid-infrared observations suggest that we are witnessing dust formation and grain evolution. Semi-regular oscillatory variability in the star's light curve is observed during the current outburst. Absorption lines develop a second blue component on a timescale twice that length. The variability may consist of one (quasi-)periodic component with P ~ 425/850 d with additional variations superimposed. During its current S Doradus outburst, R71 occupies a region in the HR diagram at the high-luminosity extension of the Cepheid instability strip and exhibits similar irregular variations as RV Tau variables. LBVs do not pass the Cepheid instability strip because of core evolution, but they develop comparable cool, low-mass, extended atmospheres in which convective instabilities may occur. As in the case of RV Tau variables, the occurrence of double absorption lines with an apparent regular cycle may be due to shocks within the atmosphere and period doubling may explain the factor of two in the lengths of the photometric and spectroscopic cycles.Comment: 18 pages, 14 figures, submitted to A&

Particle dynamics of a cartoon dune

The spatio-temporal evolution of a downsized model for a desert dune is observed experimentally in a narrow water flow channel. A particle tracking method reveals that the migration speed of the model dune is one order of magnitude smaller than that of individual grains. In particular, the erosion rate consists of comparable contributions from creeping (low energy) and saltating (high energy) particles. The saltation flow rate is slightly larger, whereas the number of saltating particles is one order of magnitude lower than that of the creeping ones. The velocity field of the saltating particles is comparable to the velocity field of the driving fluid. It can be observed that the spatial profile of the shear stress reaches its maximum value upstream of the crest, while its minimum lies at the downstream foot of the dune. The particle tracking method reveals that the deposition of entrained particles occurs primarily in the region between these two extrema of the shear stress. Moreover, it is demonstrated that the initial triangular heap evolves to a steady state with constant mass, shape, velocity, and packing fraction after one turnover time has elapsed. Within that time the mean distance between particles initially in contact reaches a value of approximately one quarter of the dune basis length

Crystal structures and freezing of dipolar fluids

We investigate the crystal structure of classical systems of spherical particles with an embedded point dipole at T=0. The ferroelectric ground state energy is calculated using generalizations of the Ewald summation technique. Due to the reduced symmetry compared to the nonpolar case the crystals are never strictly cubic. For the Stockmayer (i.e., Lennard-Jones plus dipolar) interaction three phases are found upon increasing the dipole moment: hexagonal, body-centered orthorhombic, and body-centered tetragonal. An even richer phase diagram arises for dipolar soft spheres with a purely repulsive inverse power law potential $\sim r^{-n}$. A crossover between qualitatively different sequences of phases occurs near the exponent $n=12$. The results are applicable to electro- and magnetorheological fluids. In addition to the exact ground state analysis we study freezing of the Stockmayer fluid by density-functional theory.Comment: submitted to Phys. Rev.

Ferromagnetic Liquid Thin Films Under Applied Field

Theoretical calculations, computer simulations and experiments indicate the possible existence of a ferromagnetic liquid state, although definitive experimental evidence is lacking. Should such a state exist, demagnetization effects would force a nontrivial magnetization texture. Since liquid droplets are deformable, the droplet shape is coupled with the magnetization texture. In a thin-film geometry in zero applied field, the droplet has a circular shape and a rotating magnetization texture with a point vortex at the center. We calculate the elongation and magnetization texture of such ferromagnetic thin film liquid droplet confined between two parallel plates under a weak applied magnetic field. The vortex stretches into a domain wall and exchange forces break the reflection symmetry. This behavior contrasts qualitatively and quantitatively with the elongation of paramagnetic thin films.Comment: 10 pages, 4 figures, Submitted to Phys. Rev.