
We report on our study of asteroidal breakups, i.e. fragmentations of
targets, subsequent gravitational reaccumulation and formation of small
asteroid families. We focused on parent bodies with diameters $D_{\rm pb} = 10$
km. Simulations were performed with a smoothedparticle hydrodynamics (SPH)
code combined with an efficient Nbody integrator. We assumed various
projectile sizes, impact velocities and impact angles (125 runs in total).
Resulting sizefrequency distributions are significantly different from
scaleddown simulations with $D_{\rm pb} = 100$ km targets (Durda et al.,
2007). We derive new parametric relations describing fragment distributions,
suitable for MonteCarlo collisional models. We also characterize velocity
fields and angular distributions of fragments, which can be used as initial
conditions for Nbody simulations of small asteroid families. Finally, we
discuss a number of uncertainties related to SPH simulations.

V746 Cas is known to be a triple system composed of a close binary with an
alternatively reported period of either 25.4d or 27.8d and a third component in
a 62000d orbit. The object was also reported to exhibit multiperiodic light
variations with periods from 0.83d to 2.50d, on the basis of which it was
classified as a slowly pulsating B star. Interest in further investigation of
this system was raised by the detection of a variable magnetic field. Analysing
spectra from four instruments, earlier published radial velocities, and several
sets of photometric observations, we arrived at the following conclusions: (1)
The optical spectrum is dominated by the lines of the Btype primary
(Teff1~16500(100) K), contributing 70% of the light in the optical region, and
a slightly cooler B tertiary (Teff3~13620(150) K). The lines of the lowmass
secondary are below our detection threshold; we estimate that it could be a
normal A or F star. (2) We resolved the ambiguity in the value of the inner
binary period and arrived at a linear ephemeris of T_super.conj.=HJD
2443838.78(81)+25.41569(42)xE. (3) The intensity of the magnetic field
undergoes a~sinusoidal variation in phase with one of the known photometric
periods, namely 2.503867(19)d, which we identify with the rotational period of
the tertiary. (4) The second photometric 1.0649524(40)d period is identified
with the rotational period of the Btype primary, but this interpretation is
much less certain and needs further verification. (5) If our interpretation of
photometric periods is confirmed, the classification of the object as a slowly
pulsating B star should be revised. (6) Applying an Nbody model to different
types of available observational data, we constrain the orbital inclination of
the inner orbit to ~60 deg to 85 deg even in the absence of eclipses, and
estimate the probable properties of the triple system and its components.

We construct an advanced model for interacting multiple stellar systems in
which we compute all trajectories with a numerical Nbody integrator, namely
the BulirschStoer from the SWIFT package. We can then derive various
observables: astrometric positions, radial velocities, minima timings (TTVs),
eclipse durations, interferometric visibilities, closure phases, synthetic
spectra, spectralenergy distribution, and even complete light curves. We use a
modified version of the WilsonDevinney code for the latter, in which the
instantaneous true phase and inclination of the eclipsing binary are governed
by the Nbody integration. If one has all kinds of observations at disposal, a
joint $\chi^2$ metric and an optimisation algorithm (a~simplex or simulated
annealing) allows to search for a global minimum and construct very robust
models of stellar systems. At the same time, our Nbody model is free from
artefacts which may arise if mutual gravitational interactions among all
components are not selfconsistently accounted for. Finally, we present a
number of examples showing dynamical effects that can be studied with our code
and we discuss how systematic errors may affect the results (and how to prevent
this from happening).

The goal of this report is to give a comprehensive overview of the rich field
of forward physics, with a special attention to the topics that can be studied
at the LHC. The report starts presenting a selection of the Monte Carlo
simulation tools currently available, chapter 2, then enters the rich
phenomenology of QCD at low, chapter 3, and high, chapter 4, momentum transfer,
while the unique scattering conditions of central exclusive production are
analyzed in chapter 5. The last two experimental topics, Cosmic Ray and Heavy
Ion physics are presented in the chapter 6 and 7 respectively. Chapter 8 is
dedicated to the BFKL dynamics, multiparton interactions, and saturation. The
report ends with an overview of the forward detectors at LHC. Each chapter is
correlated with a comprehensive bibliography, attempting to provide to the
interested reader with a wide opportunity for further studies.

The 2:1 meanmotion resonance with Jupiter harbours two distinct groups of
asteroids. The shortlived population is known to be a transient group
sustained in steady state by the Yarkovsky semimajor axis drift. The longlived
asteroids, however, can exhibit dynamical lifetimes comparable to
$4\,\mathrm{Gyr}$. They reside near two isolated islands of the phase space
denoted $\mathrm{A}$ and $\mathrm{B}$, with an uneven population ratio
$\mathrm{B}/\mathrm{A} \simeq 10$. The orbits of $\mathrm{A}$island asteroids
are predominantly highly inclined, compared to island $\mathrm{B}$. The
sizefrequency distribution is steep but the orbital distribution lacks any
evidence of a collisional cluster. These observational constraints are somewhat
puzzling and therefore the origin of the longlived asteroids has not been
explained so far.
With the aim to provide a viable explanation, we first update the resonant
population and revisit its physical properties. Using an $N$body model with
seven planets and the Yarkovsky effect included, we demonstrate that the
dynamical depletion of island $\mathrm{A}$ is faster, in comparison with island
$\mathrm{B}$. Then we investigate (i) the survivability of primordial resonant
asteroids and (ii) capture of the population during planetary migration,
following a recently described scenario with an escaping fifth giant planet and
a jumpingJupiter instability. We also model the collisional evolution of the
resonant population over past $4\,\mathrm{Gyr}$. Our conclusion is that the
longlived group was created by resonant capture from a narrow part of
hypothetical outer mainbelt family during planetary migration. Primordial
asteroids surviving the migration were probably not numerous enough to
substantially contribute to the observed population.

Infrared radiation emitted from an asteroid surface causes a torque that can
significantly affect rotational state of the asteroid. The influence of small
topographic features on this phenomenon, called the YORP effect, seems to be of
utmost importance. In this work, we show that a lateral heat diffusion in
boulders of suitable sizes leads to an emergence of a local YORP effect which
magnitude is comparable to the YORP effect due to the global shape. We solve a
threedimensional heat diffusion equation in a boulder and its surroundings by
the finite element method, using the FreeFem++ code. The contribution to the
total torque is inferred from the computed temperature distribution. Our
general approach allows us to compute the torque induced by a realistic
irregular boulder. For an idealized boulder, our result is consistent with an
existing onedimensional model. We also estimated (and extrapolated) a size
distribution of boulders on (25143) Itokawa from closeup images of its
surface. We realized that topographic features on Itokawa can potentially
induce a torque corresponding to a rotational acceleration of the order 10^7
rad day^2 and can therefore explain the observed phase shift in light curves.

Asteroids formed in a dynamically quiescent disk but their orbits became
gravitationally stirred enough by Jupiter to lead to highspeed collisions. As
a result, many dozen large asteroids have been disrupted by impacts over the
age of the Solar System, producing groups of fragments known as asteroid
families. Here we explain how the asteroid families are identified, review
their current inventory, and discuss how they can be used to get insights into
longterm dynamics of main belt asteroids. Electronic tables of the membership
for 122 notable families are reported on the Planetary Data System node.

Current amount of ~500 asteroid models derived from the diskintegrated
photometry by the lightcurve inversion method allows us to study not only the
spinvector properties of the whole population of MBAs, but also of several
individual collisional families. We create a data set of 152 asteroids that
were identified by the HCM method as members of ten collisional families, among
them are 31 newly derived unique models and 24 new models with wellconstrained
poleecliptic latitudes of the spin axes. The remaining models are adopted from
the DAMIT database or the literature. We revise the preliminary family
membership identification by the HCM method according to several additional
criteria  taxonomic type, color, albedo, maximum Yarkovsky semimajor axis
drift and the consistency with the sizefrequency distribution of each family,
and consequently we remove interlopers. We then present the spinvector
distributions for eight asteroidal families. We use a combined orbital and
spinevolution model to explain the observed spinvector properties of objects
among collisional families. In general, we observe for studied families similar
trends in the (a_p, \beta) space: (i) larger asteroids are situated in the
proximity of the center of the family; (ii) asteroids with \beta>0{\deg} are
usually found to the right from the family center; (iii) on the other hand,
asteroids with \beta<0{\deg} to the left from the center; (iv) majority of
asteroids have large poleecliptic latitudes (\beta\gtrsim 30{\deg}); and
finally (v) some families have a statistically significant excess of asteroids
with \beta>0{\deg} or \beta<0{\deg}. Our numerical simulation of the longterm
evolution of a collisional family is capable of reproducing well the observed
spinvector properties. Using this simulation, we also independently constrain
the age of families Flora (1.0\pm0.5 Gyr) and Koronis (2.54 Gyr).

The larger number of models of asteroid shapes and their rotational states
derived by the lightcurve inversion give us better insight into both the nature
of individual objects and the whole asteroid population. With a larger
statistical sample we can study the physical properties of asteroid
populations, such as mainbelt asteroids or individual asteroid families, in
more detail. Shape models can also be used in combination with other types of
observational data (IR, adaptive optics images, stellar occultations), e.g., to
determine sizes and thermal properties. We use all available photometric data
of asteroids to derive their physical models by the lightcurve inversion method
and compare the observed pole latitude distributions of all asteroids with
known convex shape models with the simulated pole latitude distributions. We
used classical dense photometric lightcurves from several sources and
sparseintime photometry from the U.S. Naval Observatory in Flagstaff,
Catalina Sky Survey, and La Palma surveys (IAU codes 689, 703, 950) in the
lightcurve inversion method to determine asteroid convex models and their
rotational states. We also extended a simple dynamical model for the spin
evolution of asteroids used in our previous paper. We present 119 new asteroid
models derived from combined dense and sparseintime photometry. We discuss
the reliability of asteroid shape models derived only from Catalina Sky Survey
data (IAU code 703) and present 20 such models. By using different values for a
scaling parameter cYORP (corresponds to the magnitude of the YORP momentum) in
the dynamical model for the spin evolution and by comparing synthetics and
observed polelatitude distributions, we were able to constrain the typical
values of the cYORP parameter as between 0.05 and 0.6.

In the Nice model, the late heavy bombardment (LHB) is related to an orbital
instability of giant planets which causes a fast dynamical dispersion of a
transneptunian cometary disk. We study effects produced by these hypothetical
cometary projectiles on mainbelt asteroids. In particular, we want to check
whether the observed collisional families provide a lower or an upper limit for
the cometary flux during the LHB.
We present an updated list of observed asteroid families as identified in the
space of synthetic proper elements by the hierarchical clustering method,
colour data, albedo data and dynamical considerations and we estimate their
physical parameters. We selected 12 families which may be related to the LHB
according to their dynamical ages. We then used collisional models and Nbody
orbital simulations to gain insight into the longterm dynamical evolution of
synthetic LHB families over 4 Gyr. We account for the mutual collisions, the
physical disruptions of comets, the Yarkovsky/YORP drift, chaotic diffusion, or
possible perturbations by the giantplanet migration.
Assuming a "standard" sizefrequency distribution of primordial comets, we
predict the number of families with parentbody sizes D_PB >= 200 km which
seems consistent with observations. However, more than 100 asteroid families
with D_PB >= 100 km should be created at the same time which are not observed.
This discrepancy can be nevertheless explained by the following processes: i)
asteroid families are efficiently destroyed by comminution (via collisional
cascade), ii) disruptions of comets below some critical perihelion distance (q
<~ 1.5 AU) are common.
Given the freedom in the cometarydisruption law, we cannot provide stringent
limits on the cometary flux, but we can conclude that the observed distribution
of asteroid families does not contradict with a cometary LHB.

The ALICE Collaboration has measured inclusive J/psi production in pp
collisions at a center of mass energy sqrt(s)=2.76 TeV at the LHC. The results
presented in this Letter refer to the rapidity ranges y<0.9 and 2.5<y<4 and
have been obtained by measuring the electron and muon pair decay channels,
respectively. The integrated luminosities for the two channels are L^e_int=1.1
nb^1 and L^mu_int=19.9 nb^1, and the corresponding signal statistics are
N_J/psi^e+e=59 +/ 14 and N_J/psi^mu+mu=1364 +/ 53. We present
dsigma_J/psi/dy for the two rapidity regions under study and, for the forwardy
range, d^2sigma_J/psi/dydp_t in the transverse momentum domain 0<p_t<8 GeV/c.
The results are compared with previously published results at sqrt(s)=7 TeV and
with theoretical calculations.

We model the longterm evolution of the Hilda collisional family located in
the 3/2 meanmotion resonance with Jupiter. Its eccentricity distribution
evolves mostly due to the Yarkovsky/YORP effect and assuming that: (i) impact
disruption was isotropic, and (ii) albedo distribution of small asteroids is
the same as for large ones, we can estimate the age of the Hilda family to be
$4_{1}^{+0}\,{\rm Gyr}$. We also calculate collisional activity in the J3/2
region. Our results indicate that current collisional rates are very low for a
200\,km parent body such that the number of expected events over Gyrs is much
smaller than one.
The large age and the low probability of the collisional disruption lead us
to the conclusion that the Hilda family might have been created during the Late
Heavy Bombardment when the collisions were much more frequent. The Hilda family
may thus serve as a test of orbital behavior of planets during the LHB. We
tested the influence of the giantplanet migration on the distribution of the
family members. The scenarios that are consistent with the observed Hilda
family are those with fast migration time scales $\simeq 0.3\,{\rm Myr}$ to
$3\,{\rm Myr}$, because longer time scales produce a family that is depleted
and too much spread in eccentricity. Moreover, there is an indication that
Jupiter and Saturn were no longer in a compact configuration (with period ratio
$P_{\rm S}/P_{\rm J} > 2.09$) at the time when the Hilda family was created.

The Atype star HR 6412 = V2368 Oph was used by several investigators as a
photometric comparison star for the known eclipsing binary U Oph but was found
to be variable by three independent groups, including us. By analysing series
of new spectral and photometric observations and a critical compilation of
available radial velocities, we were able to find the correct period of light
and radialvelocity variations and demonstrate that the object is an eclipsing
and doublelined spectroscopic binary moving in a highly eccentric orbit. We
derived a linear ephemeris T min.I = HJD (2454294.67 +/ 0.01) + (38.32712 +/
0.00004)d x E and estimated preliminary basic physical properties of the
binary. The dereddened UBV magnitudes and effective temperatures of the primary
and secondary, based on our light and velocitycurve solutions, led to
distance estimates that agree with the Hipparcos distance within the errors. We
find that the mass ratio must be close to one, but the limited number and
wavelength range of our current spectra does not allow a truly precise
determination of the binary masses. Nevertheless, our results show convincingly
that both binary components are evolved away from the main sequence, which
makes this system astrophysically very important. There are only a few
similarly evolved Atype stars among known eclipsing binaries. Future
systematic observations and careful analyses can provide very stringent tests
for the stellar evolutionary theory.

The multiple system V505 Sagittarii is composed of at least three stars: a
compact eclipsing pair and a distant component, which orbit is measured
directly using speckle interferometry. In order to explain the observed orbit
of the third body in V505 Sagittarii and also other observable quantities,
namely the minima timings of the eclipsing binary and two different radial
velocities in the spectrum, we thoroughly test a fourthbody hypothesis  a
perturbation by a dim, yetunobserved object. We use an Nbody numerical
integrator to simulate future and past orbital evolution of 3 or 4 components
in this system. We construct a suitable chi^2 metric from all available
speckleinterferometry, minimatimings and radialvelocity data and we scan a
part of a parameter space to get at least some of allowed solutions. In
principle, we are able to explain all observable quantities by a presence of a
fourth body, but the resulting likelihood of this hypothesis is very low. We
also discuss other theoretical explanations of the minima timings variations.
Further observations of the minima timings during the next decade or
highresolution spectroscopic data can significantly constrain the model.