Rules for transition rates in nonequilibrium steady states

Just as transition rates in a canonical ensemble must respect the principle of detailed balance, constraints exist on transition rates in driven steady states. I derive those constraints, by maximum information-entropy inference, and apply them to the steady states of driven diffusion and a sheared lattice fluid. The resulting ensemble can potentially explain nonequilibrium phase behaviour and, for steady shear, gives rise to stress-mediated long-range interactions.Comment: 4 pages. To appear in Physical Review Letter

Langevin Dynamics simulations of a 2-dimensional colloidal crystal under confinement and shear

Langevin Dynamics simulations are used to study the effect of shear on a two-dimensional colloidal crystal confined by structured parallel walls. When walls are sheared very slowly, only two or three crystalline layers next to the walls move along with them, while the inner layers of the crystal are only slightly tilted. At higher shear velocities, this inner part of the crystal breaks into several pieces with different orientations. The velocity profile across the slit is reminiscent of shear-banding in flowing soft materials, where liquid and solid regions coexist; the difference, however, is that in the latter case the solid regions are glassy while here they are crystalline. At even higher shear velocities, the effect of the shearing becomes smaller again. Also the effective temperature near the walls (deduced from the velocity distributions of the particles) decreases again when the wall velocity gets very large. When the walls are placed closer together, thereby introducing a misfit, a structure containing a soliton staircase arises in simulations without shear. Introducing shear increases the disorder in these systems until no solitons are visible any more. Instead, similar structures like in the case without misfit result. At high shear rates, configurations where the incommensurability of the crystalline structure is compensated by the creation of holes become relevant

A model colloidal fluid with competing interactions: bulk and interfacial properties

Using a simple mean-field density functional theory theory (DFT), we investigate the structure and phase behaviour of a model colloidal fluid composed of particles interacting via a pair potential which has a hard core of diameter $\sigma$, is attractive Yukawa at intermediate separations and repulsive Yukawa at large separations. We analyse the form of the asymptotic decay of the bulk fluid correlation functions, comparing results from our DFT with those from the self consistent Ornstein-Zernike approximation (SCOZA). In both theories we find rich crossover behaviour, whereby the ultimate decay of correlation functions changes from monotonic to long-wavelength damped oscillatory decay on crossing certain lines in the phase diagram, or sometimes from oscillatory to oscillatory with a longer wavelength. For some choices of potential parameters we find, within the DFT, a $\lambda$-line at which the fluid becomes unstable with respect to periodic density fluctuations. SCOZA fails to yield solutions for state points near such a $\lambda$-line. The propensity to clustering of particles, which is reflected by the presence of a long wavelength $\gg \sigma$, slowly decaying oscillatory pair correlation function, and a structure factor that exhibits a very sharp maximum at small but non zero wavenumbers, is enhanced in states near the $\lambda$-line. We present density profiles for the planar liquid-gas interface and for fluids adsorbed at a planar hard wall. The presence of a nearby $\lambda$-transition gives rise to pronounced long-wavelength oscillations in the one-body densities at both types of interface.Comment: 14 pages, 11 figure

Critical behaviors of sheared frictionless granular materials near jamming transition

Critical behaviors of sheared dense and frictionless granular materials in the vicinity of the jamming transition are numerically investigated. From the extensive molecular dynamics simulation, we verify the validity of the scaling theory near the jamming transition proposed by Otsuki and Hayakawa (Prog. Theor. Phys., 121, 647 (2009)). We also clarify the critical behaviors of the shear viscosity and the pair correlation function based on both a phenomenology and the simulation.Comment: 13pages, 26 figure

An XMM-Newton observation of the young open cluster NGC 2547: coronal activity at 30 Myr

We report XMM-Newton observations of the young open cluster NGC 2547 which allow us to characterise coronal activity in solar-type stars at an age of 30 Myr. X-ray emission peaks among G-stars at luminosities (0.3-3keV) of Lx~10^{30.5} erg/s and declines to Lx<=10^{29.0} erg/s among M-stars. Coronal spectra show evidence for multi-temperature differential emission measures and low coronal metal abundances (Z~0.3). The G- and K-type stars follow the same relationship between X-ray activity and Rossby number established in older clusters and field stars, although most solar-type stars in NGC 2547 exhibit saturated/super-saturated X-ray activity levels. Median levels of Lx and Lx/Lbol in the solar-type stars of NGC 2547 are similar to T-Tauri stars of the Orion Nebula cluster (ONC), but an order of magnitude higher than in the older Pleiades. The spread in X-ray activity levels among solar-type stars in NGC 2547 is much smaller than in older or younger clusters. Coronal temperatures increase with Lx, Lx/Lbol and surface X-ray flux. Active solar-type stars in NGC 2547 have coronal temperatures between those in the ONC and the most active older ZAMS stars. A flaring rate (for total flare energies [0.3-3keV] >10^{34} erg) of 1 every 350^{+350}_{-120} ks was found for solar-type stars, similar to rates found in the ONC and Pleiades. Comparison with ROSAT HRI data taken 7 years previously reveals that only 10-15 percent of solar-type stars or stars with Lx>3x10^{29} erg/s exhibit X-ray variability by more than a factor of two. The similar levels of X-ray activity and rate of occurrence for large flares in NGC 2547 and the ONC demonstrate that the X-ray radiation environment around young solar-type stars remains relatively constant over their first 30 Myr (abridged).Comment: Accepted for publication in MNRAS. Electronic tables available from the autho

Fluctuation theorem for stochastic systems

The fluctuation theorem describes the probability ratio of observing trajectories that satisfy or violate the second law of thermodynamics. It has been proved in a number of different ways for thermostatted deterministic nonequilibrium systems. In the present paper we show that the fluctuation theorem is also valid for a class of stochastic nonequilibrium systems. The theorem is therefore not reliant on the reversibility or the determinism of the underlying dynamics. Numerical tests verify the theoretical result

Steady shear flow thermodynamics based on a canonical distribution approach

A non-equilibrium steady state thermodynamics to describe shear flows is developed using a canonical distribution approach. We construct a canonical distribution for shear flow based on the energy in the moving frame using the Lagrangian formalism of the classical mechanics. From this distribution we derive the Evans-Hanley shear flow thermodynamics, which is characterized by the first law of thermodynamics $dE = T dS - Q d\gamma$ relating infinitesimal changes in energy $E$, entropy $S$ and shear rate $\gamma$ with kinetic temperature $T$. Our central result is that the coefficient $Q$ is given by Helfand's moment for viscosity. This approach leads to thermodynamic stability conditions for shear flow, one of which is equivalent to the positivity of the correlation function of $Q$. We emphasize the role of the external work required to sustain the steady shear flow in this approach, and show theoretically that the ensemble average of its power $\dot{W}$ must be non-negative. A non-equilibrium entropy, increasing in time, is introduced, so that the amount of heat based on this entropy is equal to the average of $\dot{W}$. Numerical results from non-equilibrium molecular dynamics simulation of two-dimensional many-particle systems with soft-core interactions are presented which support our interpretation.Comment: 23 pages, 7 figure

A census of massive stars in NGC 346. Stellar parameters and rotational velocities

Spectroscopy for 247 stars towards the young cluster NGC 346 in the Small Magellanic Cloud has been combined with that for 116 targets from the VLT-FLAMES Survey of Massive Stars. Spectral classification yields a sample of 47 O-type and 287 B-type spectra, while radial-velocity variations and/or spectral multiplicity have been used to identify 45 candidate single-lined systems, 17 double-lined systems, and one triple-lined system. Atmospheric parameters (T$_eff$ and log$g$) and projected rotational velocities ($v_e$sin$i$) have been estimated using TLUSTY model atmospheres; independent estimates of $v_e$sin$i$ were also obtained using a Fourier Transform method. Luminosities have been inferred from stellar apparent magnitudes and used in conjunction with the T$_eff$ and $v_e$sin$i$ estimates to constrain stellar masses and ages using the BONNSAI package. We find that targets towards the inner region of NGC 346 have higher median masses and projected rotational velocities, together with smaller median ages than the rest of the sample. There appears to be a population of very young targets with ages of less than 2 Myr, which have presumably all formed within the cluster. The more massive targets are found to have lower $v_e$sin$i$ consistent with previous studies. No significant evidence is found for differences with metallicity in the stellar rotational velocities of early-type stars, although the targets in the SMC may rotate faster than those in young Galactic clusters. The rotational velocity distribution for single non-supergiant B-type stars is inferred and implies that a significant number have low rotational velocity ($\simeq$10\% with $v_e$<40 km/s), together with a peak in the probability distribution at $v_e \simeq$300 km/s. Larger projected rotational velocity estimates have been found for our Be-type sample and imply that most have rotational velocities between 200-450 km/s.Comment: Accepted by A&