
The first observing run of Advanced LIGO spanned 4 months, from September 12,
2015 to January 19, 2016, during which gravitational waves were directly
detected from two binary black hole systems, namely GW150914 and GW151226.
Confident detection of gravitational waves requires an understanding of
instrumental transients and artifacts that can reduce the sensitivity of a
search. Studies of the quality of the detector data yield insights into the
cause of instrumental artifacts and data quality vetoes specific to a search
are produced to mitigate the effects of problematic data. In this paper, the
systematic removal of noisy data from analysis time is shown to improve the
sensitivity of searches for compact binary coalescences. The output of the
PyCBC pipeline, which is a pythonbased code package used to search for
gravitational wave signals from compact binary coalescences, is used as a
metric for improvement. GW150914 was a loud enough signal that removing noisy
data did not improve its significance. However, the removal of data with excess
noise decreased the false alarm rate of GW151226 by more than two orders of
magnitude, from 1 in 770 years to less than 1 in 186000 years.

We study lithium depletion in lowmass and solarlike stars as a function of
time, using a new diffusion coefficient describing extramixing taking place at
the bottom of a convective envelope. This new form is motivated by
multidimensional fully compressible, time implicit hydrodynamic simulations
performed with the MUSIC code. Intermittent convective mixing at the convective
boundary in a star can be modeled using extreme value theory, a statistical
analysis frequently used for finance, meteorology, and environmental science.
In this letter, we implement this statistical diffusion coefficient in a
onedimensional stellar evolution code, using parameters calibrated from
multidimensional hydrodynamic simulations of a young lowmass star. We propose
a new scenario that can explain observations of the surface abundance of
lithium in the Sun and in clusters covering a wide range of ages, from $\sim$
50 Myr to $\sim$ 4 Gyr. Because it relies on our physical model of convective
penetration, this scenario has a limited number of assumptions. It can explain
the observed trend between rotation and depletion, based on a single additional
assumption, namely that rotation affects the mixing efficiency at the
convective boundary. We suggest the existence of a threshold in stellar
rotation rate above which rotation strongly prevents the vertical penetration
of plumes and below which rotation has small effects. In addition to providing
a possible explanation for the long standing problem of lithium depletion in
premain sequence and main sequence stars, the strength of our scenario is that
its basic assumptions can be tested by future hydrodynamic simulations.

We examine a penetration layer formed between a central radiative zone and a
large convection zone in the deep interior of a young lowmass star. Using the
Multidimensional Stellar Implicit Code (MUSIC) to simulate twodimensional
compressible stellar convection in a spherical geometry over long times, we
produce statistics that characterize the extent and impact of convective
penetration in this layer. We apply extreme value theory to the maximal extent
of convective penetration at any time. We compare statistical results from
simulations which treat nonlocal convection, throughout a large portion of the
stellar radius, with simulations designed to treat local convection in a small
region surrounding the penetration layer. For each of these situations, we
compare simulations of different resolution, which have different velocity
magnitudes. We also compare statistical results between simulations that
radiate energy at a constant rate to those that allow energy to radiate from
the stellar surface according to the local surface temperature. Based on the
frequency and depth of penetrating convective structures, we observe two
distinct layers that form between the convection zone and the stable radiative
zone. We show that the probability density function of the maximal depth of
convective penetration at any time corresponds closely in space with the radial
position where internal waves are excited. We find that the maximal penetration
depth can be modeled by a Weibull distribution with a small shape parameter.
Using these results, and building on established scalings for diffusion
enhanced by largescale convective motions, we propose a new form for the
diffusion coefficient that may be used for onedimensional stellar evolution
calculations in the large P\'eclet number regime. These results should
contribute to the 321D link.

We investigate the utility of the convex hull of many Lagrangian tracers to
analyze transport properties of turbulent flows with different anisotropy. In
direct numerical simulations of statistically homogeneous and stationary
NavierStokes turbulence, neutral fluid Boussinesq convection, and MHD
Boussinesq convection a comparison with Lagrangian pair dispersion shows that
convex hull statistics capture the asymptotic dispersive behavior of a large
group of passive tracer particles. Moreover, convex hull analysis provides
additional information on the subensemble of tracers that on average disperse
most efficiently in the form of extreme value statistics and flow anisotropy
via the geometric properties of the convex hulls. We use the convex hull
surface geometry to examine the anisotropy that occurs in turbulent convection.
Applying extreme value theory, we show that the maximal square extensions of
convex hull vertices are well described by a classic extreme value
distribution, the Gumbel distribution. During turbulent convection,
intermittent convective plumes grow and accelerate the dispersion of Lagrangian
tracers. Convex hull analysis yields information that supplements standard
Lagrangian analysis of coherent turbulent structures and their influence on the
global statistics of the flow.

The first direct gravitationalwave detection was made by the Advanced Laser
Interferometer Gravitational Wave Observatory on September 14, 2015. The
GW150914 signal was strong enough to be apparent, without using any waveform
model, in the filtered detector strain data. Here, features of the signal
visible in the data are analyzed using concepts from Newtonian physics and
general relativity, accessible to anyone with a general physics background. The
simple analysis presented here is consistent with the fully
generalrelativistic analyses published elsewhere,in showing that the signal
was produced by the inspiral and subsequent merger of two black holes. The
black holes were each of approximately 35 Msun, still orbited each other as
close as ~350 km apart, and subsequently merged to form a single black hole.
Similar reasoning, directly from the data, is used to roughly estimate how far
these black holes were from the Earth, and the energy that they radiated in
gravitational waves.

We present the results of a numerical benchmark study for the
MUltidimensional Stellar Implicit Code (MUSIC) based on widely applicable two
and threedimensional compressible hydrodynamics problems relevant to stellar
interiors. MUSIC is an implicit large eddy simulation code that uses implicit
time integration, implemented as a Jacobianfree Newton Krylov method. A
physics based preconditioning technique which can be adjusted to target varying
physics is used to improve the performance of the solver. The problems used for
this benchmark study include the RayleighTaylor and KelvinHelmholtz
instabilities, and the decay of the TaylorGreen vortex. Additionally we show a
test of hydrostatic equilibrium, in a stellar environment which is dominated by
radiative effects. In this setting the flexibility of the preconditioning
technique is demonstrated. This work aims to bridge the gap between the
hydrodynamic test problems typically used during development of numerical
methods and the complex flows of stellar interiors. A series of
multidimensional tests are performed and analysed. Each of these test cases is
analysed with a simple, scalar diagnostic, with the aim of enabling direct code
comparisons. As the tests performed do not have analytic solutions we verify
MUSIC by comparing to established codes including ATHENA and the PENCIL code.
MUSIC is able to both reproduce behaviour from established and widelyused
codes as well as results expected from theoretical predictions. This
benchmarking study concludes a series of papers describing the development of
the MUSIC code and provides confidence in the future applications.

When electron cyclotron (EC) driven current is first applied to the inside of
a magnetic island, the current spreads throughout the island and after a short
period achieves a steady level. Using a two equation fluid model for the EC
current that allows us to examine this early evolution in detail, we analyze
highresolution simulations of a 2/1 classical tearing mode in a lowbeta large
aspectratio circular tokamak. These simulations use a nonlinear 3D reducedMHD
fluid model and the JOREK code. During the initial period where the EC driven
current grows and spreads throughout the magnetic island, it is not a function
of the magnetic flux. However, once it has reached a steadystate, it should be
a flux function. We demonstrate numerically that if sufficiently resolved
toroidally, the steadystate EC driven current becomes approximately a flux
function. We discuss the physics of this early period of EC evolution and its
impact on the size of the magnetic island.

The secondgeneration of gravitationalwave detectors are just starting
operation, and have already yielding their first detections. Research is now
concentrated on how to maximize the scientific potential of gravitationalwave
astronomy. To support this effort, we present here design targets for a new
generation of detectors, which will be capable of observing compact binary
sources with high signaltonoise ratio throughout the Universe.

We report here the nondetection of gravitational waves from the merger of
binary neutron star systems and neutronstarblackhole systems during the
first observing run of Advanced LIGO. In particular we searched for
gravitational wave signals from binary neutron star systems with component
masses $\in [1,3] M_{\odot}$ and component dimensionless spins $< 0.05$. We
also searched for neutronstarblackhole systems with the same neutron star
parameters, black hole mass $\in [2,99] M_{\odot}$ and no restriction on the
black hole spin magnitude. We assess the sensitivity of the two LIGO detectors
to these systems, and find that they could have detected the merger of binary
neutron star systems with component mass distributions of $1.35\pm0.13
M_{\odot}$ at a volumeweighted average distance of $\sim$ 70Mpc, and for
neutronstarblackhole systems with neutron star masses of $1.4M_\odot$ and
black hole masses of at least $5M_\odot$, a volumeweighted average distance of
at least $\sim$ 110Mpc. From this we constrain with 90% confidence the merger
rate to be less than 12,600 Gpc$^{3}$yr$^{1}$ for binaryneutron star systems
and less than 3,600 Gpc$^{3}$yr$^{1}$ for neutronstarblackhole systems.
We find that if no detection of neutronstar binary mergers is made in the next
two Advanced LIGO and Advanced Virgo observing runs we would place significant
constraints on the merger rates. Finally, assuming a rate of
$10^{+20}_{7}$Gpc$^{3}$yr$^{1}$ short gamma ray bursts beamed towards the
Earth and assuming that all short gammaray bursts have binaryneutronstar
(neutronstarblackhole) progenitors we can use our 90% confidence rate upper
limits to constrain the beaming angle of the gammaray burst to be greater than
${2.3^{+1.7}_{1.1}}^{\circ}$ (${4.3^{+3.1}_{1.9}}^{\circ}$).

Context: We study the impact of twodimensional spherical shells on
compressible convection. Realistic profiles for density and temperature from a
onedimensional stellar evolution code are used to produce a model of a large
stellar convection zone representative of a young lowmass star.
Methods: We perform hydrodynamic implicit largeeddy simulations of
compressible convection using the MUltidimensional Stellar Implicit Code
(MUSIC). Because MUSIC has been designed to use realistic stellar models
produced from onedimensional stellar evolution calculations, MUSIC simulations
are capable of seamlessly modeling a whole star. Simulations in twodimensional
spherical shells that have different radial extents are performed over hundreds
of convective turnover times, permitting the collection of wellconverged
statistics.
Results: We evaluate basic statistics of the convective turnover time, the
convective velocity, and the overshooting layer. These quantities are selected
for their relevance to onedimensional stellar evolution calculations, so that
our results are focused toward the 321D link. The inclusion in the spherical
shell of the boundary between the radiative and convection zones decreases the
amplitude of convective velocities in the convection zone. The inclusion of
nearsurface layers in the spherical shell can increase the amplitude of
convective velocities, although the radial structure of the velocity profile
established by deep convection is unchanged. The impact from including the
nearsurface layers depends on the speed and structure of smallscale
convection in the nearsurface layers. Larger convective velocities in the
convection zone result in a commensurate increase in the overshooting layer
width and decrease in the convective turnover time. These results provide
support for nonlocal aspects of convection.

This work is the first attempt to describe the multidimensional structure of
accreting young stars based on fully compressible time implicit
multidimensional hydrodynamics simulations. One major motivation is to analyse
the validity of accretion treatment used in previous 1D stellar evolution
studies. We analyse the effect of accretion on the structure of a realistic
stellar model of the young Sun. Our work is inspired by the numerical work of
Kley \& Lin (1996, ApJ, 461, 933) devoted to the structure of the boundary
layer in accretion disks. We analyse the redistribution of accreted material
with a range of values of specific entropy relative to the bulk specific
entropy of the material in the accreting object's convective envelope. A
primary goal is to understand whether and how accreted energy deposited onto a
stellar surface is redistributed in the interior. This study focusses on the
high accretion rates characteristic of FU Ori systems. We find that the highest
entropy cases produce a distinctive behaviour in the mass redistribution, rms
velocities, and enthalpy flux in the convective envelope. This change in
behaviour is characterised by the formation of a hot layer on the surface of
the accreting object, which tends to suppress convection in the envelope. We
analyse the longterm effect of such a hot buffer zone on the structure and
evolution of the accreting object with 1D stellar evolution calculations. We
study the relevance of the assumption of redistribution of accreted energy into
the stellar interior used in the literature. One conclusion is that, for a
given amount of accreted energy transferred to the accreting object, a
treatment assuming accretion energy redistribution throughout the stellar
interior could significantly overestimate the effects on the stellar structure,
in particular, on the resulting expansion.

Local regions of anomalous particle dispersion, and intermittent events that
occur in turbulent flows can greatly influence the global statistical
description of the flow. These local behaviors can be identified and analyzed
by comparing the growth of neighboring convex hulls of Lagrangian tracer
particles. Although in our simulations of homogeneous turbulence the convex
hulls generally grow in size, after the Lagrangian particles that define the
convex hulls begin to disperse, our analysis reveals short periods when the
convex hulls of the Lagrangian particles shrink, evidence that particles are
not dispersing simply. Shrinkage can be associated with anisotropic flows,
since it occurs most frequently in the presence of a mean magnetic field or
thermal convection. We compare dispersion between a wide range of statistically
homogeneous and stationary turbulent flows ranging from homogeneous isotropic
NavierStokes turbulence over different configurations of magnetohydrodynamic
turbulence and Boussinesq convection.

Intermittent largescale highshear flows are found to occur frequently and
spontaneously in direct numerical simulations of statistically stationary
turbulent Boussinesq magnetohydrodynamic (MHD) convection. The energetic
steadystate of the system is sustained by convective driving of the velocity
field and smallscale dynamo action. The intermittent emergence of flow
structures with strong velocity and magnetic shearing generates magnetic energy
at an elevated rate over timescales longer than the characteristic time of the
largescale convective motion. The resilience of magnetic energy amplification
suggests that intermittent shearbursts are a significant driver of dynamo
action in turbulent magnetoconvection.

We quantify possible differences between turbulent dynamo action in the Sun
and the dynamo action studied in idealized simulations. For this purpose we
compare Fourierspace shelltoshell energy transfer rates of three
incrementally more complex dynamo simulations: an incompressible, periodic
simulation driven by random flow, a simulation of Boussinesq convection, and a
simulation of fully compressible convection that includes physics relevant to
the nearsurface layers of the Sun. For each of the simulations studied, we
find that the dynamo mechanism is universal in the kinematic regime because
energy is transferred from the turbulent flow to the magnetic field from
wavenumbers in the inertial range of the energy spectrum. The addition of
physical effects relevant to the solar nearsurface layers, including
stratification, compressibility, partial ionization, and radiative energy
transport, does not appear to affect the nature of the dynamo mechanism. The
role of inertialrange shear stresses in magnetic field amplification is
independent from outerscale circumstances, including forcing and
stratification. Although the shelltoshell energy transfer functions have
similar properties to those seen in meanflow driven dynamos in each simulation
studied, the saturated states of these simulations are not universal because
the flow at the driving wavenumbers is a significant source of energy for the
magnetic field.

A kinetically stabilized axially symmetric tandem mirror (KSTM) uses the
momentum flux of lowenergy, unconfined particles that sample only the outer
endregions of the mirror plugs, where large favorable fieldline curvature
exists. The window of operation is determined for achieving MHD stability with
tolerable energy drain from the kinetic stabilizer. Then MHD stable systems are
analyzed for stability of the trapped particle mode. This mode is characterized
by the detachment of the centralcell plasma from the kinetic stabilizer region
without inducing fieldline bending. Stability of the trapped particle mode is
sensitive to the electron connection between the stabilizer and the end plug.
It is found that the stability condition for the trapped particle mode is more
constraining than the stability condition for the MHD mode, and it is
challenging to satisfy the required power constraint. Furthermore a severe
power drain may arise from the necessary connection of lowenergy electrons in
the kinetic stabilizer to the central region.