Investigating the rotational evolution of young, low mass stars using Monte Carlo simulations

We investigate the rotational evolution of young stars through Monte Carlo simulations. We simulate 280,000 stars, each of which is assigned a mass, a rotational period, and a mass accretion rate. The mass accretion rate depends on mass and time, following power-laws indices 1.4 and -1.5, respectively. A mass-dependent accretion threshold is defined below which a star is considered as diskless, which results in a distribution of disk lifetimes that matches observations. Stars are evolved at constant angular spin rate while accreting and at constant angular momentum when they become diskless. We recover the bimodal period distribution seen in several young clusters. The short period peak consists mostly of diskless stars and the long period one is mainly populated by accreting stars. Both distributions present a long tail towards long periods and a population of slowly rotating diskless stars is observed at all ages. We reproduce the observed correlations between disk fraction and spin rate, as well as between IR excess and rotational period. The period-mass relation we derive from the simulations exhibits the same global trend as observed in young clusters only if we release the disk locking assumption for the lowest mass stars. We find that the time evolution of median specific angular momentum follows a power law index of -0.65 for accreting stars and of -0.53 for diskless stars, a shallower slope that results from a wide distribution of disk lifetimes. Using observationally-documented distributions of disk lifetimes, mass accretion rates, and initial rotation periods, and evolving an initial population from 1 to 12 Myr, we reproduce the main characteristics of pre-main sequence angular momentum evolution, which supports the disk locking hypothesis. (abridged)Comment: 11 pages, 14 figures, accepted for publication in A&

Doubly-periodic array of bubbles in a Hele-Shaw cell

Exact solutions are presented for a doubly-periodic array of steadily moving bubbles in a Hele-Shaw cell when surface tension is neglected. It is assumed that the bubbles either are symmetrical with respect to the channel centreline or have fore-and-aft symmetry, or both, so that the relevant flow domain can be reduced to a simply connected region. By using conformal mapping techniques, a general solution with any number of bubbles per unit cell is obtained in integral form. Several examples are given, including solutions for multi-file arrays of bubbles in the channel geometry and doubly-periodic solutions in an unbounded cell.Comment: 15 pages, 12 figure

Multi-strange particle production in relativistic heavy ion collisions at sNN=62.4\sqrt{s_{NN}}=62.4 GeV

We present preliminary STAR results on measurements of multi-strange particles $\Xi$, $\Omega$ and their anti-particles from Au+Au and Cu+Cu at $\sqrt{s_{NN}}=62.4$ GeV collisions. In order to better understand the role of strangeness enhancement in nucleus-nucleus collisions and its scaling properties with system size, we compare the results from Au+Au and Cu+Cu reactions for different event centrality classes. Strangeness enhancement is discussed in the context of multi-strange to pion ratios. Finally, $\Omega/\phi$ ratio is shown for different systems and energies for a systematic study

Vortex motion around a circular cylinder above a plane

The study of vortex flows around solid obstacles is of considerable interest from both a theoretical and practical perspective. One geometry that has attracted renewed attention recently is that of vortex flows past a circular cylinder placed above a plane wall, where a stationary recirculating eddy can form in front of the cylinder, in contradistinction to the usual case (without the plane boundary) for which a vortex pair appears behind the cylinder. Here we analyze the problem of vortex flows past a cylinder near a wall through the lenses of the point-vortex model. By conformally mapping the fluid domain onto an annular region in an auxiliary complex plane, we compute the vortex Hamiltonian analytically in terms of certain special functions related to elliptic theta functions. A detailed analysis of the equilibria of the model is then presented. The location of the equilibrium in front of the cylinder is shown to be in qualitative agreement with the experimental findings. We also show that a topological transition occurs in phase space as the parameters of the systems are variedComment: 17 pages, 8 figure

An instance of the MIKADO migration model

In this document, we briefly describe the main contribution to the deliverable on experimenting with the implementation of most of the calculi considered in the project. First, we describe how two well known calculi for mobile processes KLAIM and Dπ have been implemented on the top of IMC. We then describe the implementation of the MiKO programming language, an instance of the parametric calculus introduced in the WP1 with the TyCO calculus as the content of the membrane itself. After this, we outline the description of the implementation of the abstract machine for an instance of the Kell Calculus that dedicates particular attention to the proof of its correctness. Our presentation ends with a discussion of the problem of implementing security membranes on the top of an execution platform

A maximum entropy approach to H-theory: Statistical mechanics of hierarchical systems

A novel formalism, called H-theory, is applied to the problem of statistical equilibrium of a hierarchical complex system with multiple time and length scales. In this approach, the system is formally treated as being composed of a small subsystem---representing the region where the measurements are made---in contact with a set of `nested heat reservoirs' corresponding to the hierarchical structure of the system. The probability distribution function (pdf) of the fluctuating temperatures at each reservoir, conditioned on the temperature of the reservoir above it, is determined from a maximum entropy principle subject to appropriate constraints that describe the thermal equilibrium properties of the system. The marginal temperature distribution of the innermost reservoir is obtained by integrating over the conditional distributions of all larger scales, and the resulting pdf is written in analytical form in terms of certain special transcendental functions, known as the Fox $H$-functions. The distribution of states of the small subsystem is then computed by averaging the quasi-equilibrium Boltzmann distribution over the temperature of the innermost reservoir. This distribution can also be written in terms of $H$-functions. The general family of distributions reported here recovers, as particular cases, the stationary distributions recently obtained by Mac\^edo {\it et al.} [Phys.~Rev.~E {\bf 95}, 032315 (2017)] from a stochastic dynamical approach to the problem.Comment: 20 pages, 2 figure