We performed the non-local thermodynamic equilibrium (non-LTE) calculations for Ti I-II with the updated model atom that includes quantum-mechanical rate coefficients for inelastic collisions with hydrogen atoms. We have calculated for the first time the rate coefficients for bound-bound transitions in inelastic collisions of titanium atoms and ions with hydrogen atoms and for the charge-exchange processes: Ti I + H <-> Ti II + H- and Ti II + H <-> Ti III + H-. The influence of these data on non-LTE abundance determinations has been tested for the Sun and four metal-poor stars. For Ti I and Ti II, the application of the derived rate coefficients has led to an increase in the departures from LTE and an increase in the titanium abundance compared to that, obtained with approximate formulas for the rate coefficients. In metal-poor stars, we have failed to achieve consistent non-LTE abundances from lines of two ionization stages. The known in the literature discrepancy in the non-LTE abundances from Ti I and Ti II lines in metal-poor stars cannot be solved by improvement of the rates of inelastic processes in collisions with hydrogen atoms in non-LTE calculations with classical model atmospheres.

Context. The pursuit of more realistic spectroscopic modelling and consistent abundances has led us to begin a new series of papers designed to improve current solar and stellar abundances of various atomic species. To achieve this, we have began updating the three-dimensional (3D) non-local thermodynamic equilibrium (non-LTE) radiative transfer code, Multi3D, and the equivalent one-dimensional (1D) non-LTE radiative transfer code, MULTI. Aims. We examine our improvements to these codes by redetermining the solar barium abundance. Barium was chosen for this test as it is an important diagnostic element of the s-process in the context of galactic chemical evolution. New Ba II + H collisional data for excitation and charge exchange reactions computed from first principles had recently become available and were included in the model atom. The atom also includes the effects of isotopic line shifts and hyperfine splitting. Method. A grid of 1D LTE barium lines were constructed with MULTI and fit to the four Ba II lines available to us in the optical region of the solar spectrum. Abundance corrections were then determined in 1D non-LTE, 3D LTE, and 3D non-LTE. A new 3D non-LTE solar barium abundance was computed from these corrections. Results. We present for the first time the full 3D non-LTE barium abundance of $A({\rm Ba})=2.27\pm0.02\pm0.01$, which was derived from four individual fully consistent barium lines. Errors here represent the systematic and random errors, respectively.

Iron plays a crucial role in studies of late-type stars. In their atmospheres, Fe I is the minority species and its lines are subject to the departures from LTE. In contrast, one believes that LTE is a realistic approximation for Fe II lines. The main source of the uncertainties in the non-LTE (NLTE) calculations for cool atmospheres is a treatment of inelastic collisions with hydrogen atoms. We investigate the effect of Fe I + H I and Fe II + H I collisions and their different treatment on the Fe I/Fe II ionisation equilibrium and Fe abundance determinations for Galactic halo benchmark stars HD84937, HD122563, and HD140283 and a sample of 38 very metal-poor giants in the dwarf galaxies with well known distances. We perform the NLTE calculations for Fe I-II with using quantum-mechanical (QM) rate coefficients for collisions with H I from Barklem (2018, B18), Yakovleva, Belyaev, and Kraemer (2018, YBK18), and Yakovleva, Belyaev, and Kraemer (2019). We find that collisions with H I serve as efficient thermalisation processes for Fe II and the NLTE abundance corrections for Fe II lines do not exceed 0.02 dex at [Fe/H] > -3 and reach +0.06~dex at [Fe/H] ~ -4. For given star, the B18 and YBK18 treatments of Fe I + H I collisions lead to similar average NLTE abundances from Fe I lines, although there exist discrepancies in the NLTE corrections for individual lines. With using QM collisional data and the Gaia based surface gravity, we obtain consistent abundances from Fe I and Fe II for a red giant HD122563. For HD84937 and HD140283, we study the Fe lines in the visible and the UV (1968-2990 A) range. For both Fe I and Fe II, abundances from the visible and UV lines are consistent. The abundances from Fe I and Fe II agree within 0.10 and 0.13 dex in the YBK18 and B18 cases. The Fe I/Fe II ionisation equilibrium is achieved for each [Fe/H] > -3.5 star of our dwarf galaxy sample.

Manganese (Mn) is a key Fe-group elements, commonly employed in stellar population and nucleosynthesis studies to explore the role of SN Ia. We have developed a new non-local thermodynamic equilibrium (NLTE) model of Mn, including new photo-ionisation cross-sections and new transition rates caused by collisions with H and H- atoms. We applied the model in combination with 1-dimensional (1D) LTE model atmospheres and 3D hydrodynamical simulations of stellar convection to quantify the impact of NLTE and convection on the line formation. We show that the effects of NLTE are present in Mn I and, to a lesser degree, in Mn II lines, and these increase with metallicity and with effective temperature of a model. Employing 3D NLTE radiative transfer, we derive new abundance of Mn in the Sun, A(Mn)=5.52 +/- 0.03 dex, consistent with the element abundance in C I meteorites. We also apply our methods to the analysis of three metal-poor benchmark stars. We find that 3D NLTE abundances are significantly higher than 1D LTE. For dwarfs, the differences between 1D NLTE and 3D NLTE abundances are typically within 0.15 dex, however, the effects are much larger in the atmospheres of giants owing to their more vigorous convection. We show that 3D NLTE successfully solves the ionisation and excitation balance for the RGB star HD 122563 that cannot be achieved by 1D LTE or 1D NLTE modelling. For HD 84937 and HD 140283, the ionisation balance is satisfied, however, the resonance Mn I triplet lines still show somewhat lower abundances compared to the high-excitation lines. Our results for the benchmark stars confirm that 1D LTE modelling leads to significant systematic biases in Mn abundances across the full wavelength range from the blue to the IR. We also produce a list of Mn lines that are not significantly biased by 3D and can be reliably, within the 0.1 dex uncertainty, modelled in 1D NLTE.

We have updated the Ba II model atom by taking into account the H I impact excitation with the rate coefficients from the quantum-mechanical calculations of Belyaev and Yakovleva (2018). Using high-resolution stellar spectra and the non-local thermodynamic equilibrium (non-LTE) line formation for Ba II, we have determined the fraction of barium isotopes with an odd mass number (fodd}) in four Galactic halo giants with well-known atmospheric parameters. We use a method based on the requirement that the abundances from the Ba II 4554 A resonance and Ba II 5853, 6496 A subordinate lines be equal. An accuracy of 0.04 dex in determining the barium abundance from individual lines has been achieved. In three stars (HD 2796, HD 108317, and HD 122563) fodd >= 0.4. This suggests that >= 80 % of the barium observed in these stars was synthesized in the r-process. In HD 128279 fodd = 0.27 exceeds the fraction of odd barium isotopes in the Solar system, but only slightly. The dominance of the r-process at the formation epoch of our sample stars is confirmed by the presence of a europium overabundance relative barium, with [Eu/Ba] > 0.3. We have calculated the non-LTE abundance corrections for five Ba II lines and investigated their dependence on atmospheric parameters in the ranges of effective temperatures from 4500 to 6500 K, surface gravities log g from 0.5 to 4.5, and metallicities [Fe/H] from 0 to -3.

We perform the non-local thermodynamic equilibrium (NLTE) calculations for Ca I-II with the updated model atom that includes new quantum-mechanical rate coefficients for Ca I + H I collisions from two recent studies, that is, by Barklem and by Mitrushchenkov, Guitou, Belyaev, Yakovleva, Spielfiedel, and Feautrier, and investigate the accuracy of calcium abundance determinations using the Sun, Procyon, and five metal-poor (MP) stars with well-determined stellar parameters. We show that both collisional recipes lead to very similar NLTE results. When using the subordinate lines of Ca I and the high-excitation lines of Ca II, NLTE provides the smaller line-to-line scatter compared with the LTE case for each star. For Procyon, NLTE removes a steep trend with line strength among strong Ca I lines seen in LTE and leads to consistent [Ca/H] abundances from the two ionisation stages. In the MP stars, the NLTE abundance from Ca II 8498 A agrees well with that from the Ca I subordinate lines. NLTE largely removes abundance discrepancies between the high-excitation lines of Ca I and Ca II 8498 A obtained for our four [Fe/H] < -2 stars under the LTE assumption. We investigate the formation of the Ca I resonance line in the [Fe/H] < -2 stars. Consistent NLTE abundances from the Ca I resonance line and the Ca II lines are found for two hyper metal-poor stars HE0107-5240 and HE1327-2326. We provide the NLTE abundance corrections for 28 lines of Ca I in a grid of model atmospheres suitable for abundance analysis of FGK-type dwarfs and subgiants.

The non-LTE line formation for Al I and Si I was calculated with model atmospheres corresponding to F-G-K type stars of different metallicity. To account for inelastic collisions with H I, for the first time we applied the cross sections calculated by Belyaev et al. using model approaches within the formalism of the Born-Oppenheimer quantum theory. For Al I non-LTE leads to overionization in the line formation layers and to weakened spectral lines, in line with earlier non-LTE studies. However, in contrast to the previuos studies, our results predict smaller magnitude of the non-LTE effects for the subordinate lines. Owing to large cross sections, the ion-pair production and mutual neutralization processes Al I(nl) + H I(1s) $\leftrightarrow$ Al~II(3s^2) + H^- provide a close coupling of high-excitation Al I levels to the Al II ground state, which causes smaller deviations from the TE populations compared to the case of pure electron collisions. For three metal-poor stars, the Al abundance was determined from seven lines in different models of their formation. In LTE and in non-LTE calculations with pure electron collisions, the Al I 3961 A resonance line gives a 0.25-0.45 dex lower abundance than the mean from the subordinate lines. The difference is removed by accounting for the collisions with H I, and the rms error of the abundance from all seven lines decreases by a factor of 1.5-3 compared to LTE. We calculated the non-LTE abundance corrections for six Al I lines depending on Teff (4500-6500 K), log g (3.0-4.5), and [Fe/H] (0, -1, -2, -3). For Si I including the collisions with H I leads to TE populations in the line formation layers and to minor departures from LTE in spectral lines.

We have explored Natural Supersymmetry (NSUSY) scenarios with low values of the $\mu$ parameter which are characterised by higgsino-like Dark Matter (DM) and compressed spectra for the lightest MSSM particles, $\chi^0_1$, $\chi^0_2$ and $\chi^\pm_1$. This scenario could be probed via monojet signatures, but as the signal-to-background ratio (S/B) is low we demonstrate that the 8 TeV LHC cannot obtain limits on the DM mass beyond those of LEP2. On the other hand, we have found, for the 13 TeV run of the LHC, that by optimising kinematical cuts we can bring the S/B ratio up to the 5(3)% level which would allow the exclusion of the DM mass up to 200(250) GeV respectively, significantly extending LEP2 limits. Moreover, we have found that LUX/XENON1T and LHC do play very complementary roles in exploring the parameter space of NSUSY, as the LHC has the capability to access regions where DM is quasi-degenerate with other higgsinos, which are challenging for direct detection experiments.

Mg is often traced in late-type stars using lines of neutral magnesium, which is expected to be subject to departures from LTE. The astrophysical importance of Mg as well as its relative simplicity from an atomic physics point of view, makes it a prime target and test bed for detailed ab initio non-LTE modelling in stellar atmospheres. For the low-lying states of Mg i, electron collision data were calculated using the R-matrix method. Calculations for collisional broadening by neutral hydrogen were also performed where data were missing. These calculations, together with data from the literature, were used to build a model atom. First, the modelling was tested by comparisons with observed spectra of benchmark stars with well-known parameters. Second, the spectral line behaviour and uncertainties were explored by extensive experiments in which sets of collisional data were changed or removed. The modelled spectra agree well with observed spectra. The line-to-line scatter in the derived abundances shows improvements compared to LTE. The observed Mg emission features at 7 and 12 microns in the spectra of the Sun and Arcturus are reasonably well reproduced. The modelling predicts non-LTE abundance corrections in dwarfs, both solar metallicity and metal-poor, to be very small (of order 0.01 dex), even smaller than found in previous studies. In giants, corrections vary greatly between lines, but can be as large as 0.4 dex. Our results emphasise the need for accurate data of Mg collisions with both electrons and H atoms for precise non-LTE predictions of stellar spectra, but demonstrate that such data can be calculated and that ab initio non-LTE modelling without resort to free parameters is possible. Grids of departure coefficients and abundance corrections for a range of stellar parameters are planned for a forthcoming paper.

Rate coefficients for inelastic Mg+H collisions are calculated for all transitions between the lowest seven levels and the ionic state (charge transfer). The rate coefficients are based on cross-sections from full quantum scattering calculations, which are themselves based on detailed quantum chemical calculations for the MgH molecule. The data are needed for non-LTE applications in cool astrophysical environments, especially cool stellar atmospheres, and are presented for a temperature range of 500-8000 K. From consideration of the sensitivity of the cross-sections to various uncertainties in the calculations, most importantly input quantum chemical data and the numerical accuracy of the scattering calculations, a measure of the possible uncertainties in the rate coefficients is estimated.

One of the main channels which allows for a large rate of neutralino dark matter annihilation in the early Universe is via the pseudoscalar Higgs A-resonance. In this case, the measured dark matter abundance can be obtained in the minimal supergravity (mSUGRA) model when tan(beta)\sim 50 and 2m_{\tz_1}\sim m_A. We investigate the reaction pp\to b\phi\to b\mu^+\mu^- +X (where \phi =A or H) at the CERN LHC where requiring the tag of a single b-jet allows for amplification of the signal-to-background ratio. The rare but observable Higgs decay to muon pairs allows for a precise measurement of the Higgs boson mass and decay width. We evaluate signal and background using CalcHEP, with muon energy smearing according to the CMS detector. We find that the Higgs width (\Gamma_A) can typically be determined with the accuracy up to \sim 8% (\sim 17%) for m_A\sim 400 (600) GeV assuming 10^3 fb^{-1} of integrated luminosity. Therefore, the pp\to b\phi\to b\mu^+\mu^- +X process provides a unique possibility for direct \Gamma_A measurement at the LHC. While the Higgs width is correlated with the parameter \tan\beta for a given value of m_A, extracting \tan\beta is complicated by an overlap of the A and H peaks, radiative corrections to the b and \tau Yukawa couplings, and the possibility that SUSY decay modes of the Higgs may be open. In the case where a dilepton mass edge from \tz_2\to\ell^+\ell^-\tz_1 is visible, it should be possible to test the relation that 2m_{\tz_1}\sim m_A.

The influence of inelastic hydrogen atom collisions on non-LTE spectral line formation has been, and remains to be, a significant source of uncertainty for stellar abundance analyses, due to the difficulty in obtaining accurate data for low-energy atomic collisions either experimentally or theoretically. For lack of a better alternative, the classical "Drawin formula" is often used. Over recent decades, our understanding of these collisions has improved markedly, predominantly through a number of detailed quantum mechanical calculations. In this paper, the Drawin formula is compared with the quantum mechanical calculations both in terms of the underlying physics and the resulting rate coefficients. It is shown that the Drawin formula does not contain the essential physics behind direct excitation by H atom collisions, the important physical mechanism being quantum mechanical in character. Quantitatively, the Drawin formula compares poorly with the results of the available quantum mechanical calculations, usually significantly overestimating the collision rates by amounts that vary markedly between transitions.

We study the potential of International Linear Collider (ILC) at sqrt(s) = 500 GeV to probe new dark matter motivated scenario where the bottom squark (sbottom) is the next to lightest supersymmetric particle. For this scenario, which is virtually impossible for the LHC to test, the ILC has a potential to cover a large fraction of the parameter space. The challenge is due to a very low energy of jets, below 20-30 GeV, which pushes the jet clustering and flavour tagging algorithms to their limits. The process of sbottom pair production was studied within the SiD detector concept. We demonstrate that ILC offers a unique opportunity to test the SUSY parameter space motivated by the sbottom neutralino co-annihilation scenario in cases when the sbottom production is kinematically accessible. The study was done with the full SiD simulation and reconstruction chain including all Standard Model and beam backgrounds.

Neutral sodium is a minority species in the atmospheres of late-type stars, and line formation in local thermodynamic equilibrium (LTE) is often a poor assumption, in particular for strong lines. We present an extensive grid of non-LTE calculations for several NaI lines in cool stellar atmospheres, including metal-rich and metal-poor dwarfs and giants. For the first time, we constructed a Na model atom that incorporates accurate quantum mechanical calculations for collisional excitation and ionisation by electrons as well as collisional excitation and charge exchange reactions with neutral hydrogen. Similar to LiI, the new rates for hydrogen impact excitation do not affect the statistical equilibrium calculations, while charge exchange reactions have a small but non-negligible influence. The presented LTE and non-LTE curves-of-growth can be interpolated to obtain non-LTE abundances and abundance corrections for arbitrary stellar parameter combinations and line strengths. The typical corrections for weak lines are -0.1...-0.2dex, whereas saturated lines may overestimate the abundance in LTE by more than 0.5dex. The non-LTE Na abundances appear very robust with respect to uncertainties in the input collisional data.

Rate coefficients for inelastic Na+H collisions are calculated for all transitions between the ten levels up to and including the ionic state (ion-pair production), namely Na(3s,3p,4s,3d,4p,5s,4d,4f,5p)+H(1s) and Na+ + H-. The calculations are based on recent full quantum scattering cross-section calculations. The data are needed for non-LTE applications in cool astrophysical environments, especially cool stellar atmospheres, and are presented for a temperature range of 500-8000 K. From consideration of the sensitivity of the cross-sections to input quantum chemical data and the results of different methods for the scattering calculations, a measure of the possible uncertainties in the rate coefficients is estimated.

Within the minimal supergravity model (mSUGRA) framework, the expectation for the relic density of neutralinos exceeds the WMAP determination unless neutralinos a) have a significant higgsino component, b) have a mass close to half that of a heavy Higgs boson, or c)can efficiently co-annihilate with a charged or colored particle. Within a 1-parameter extension of the mSUGRA model which includes non-universal Higgs masses, we show that agreement with the WMAP data can be obtained over a wide range of mSUGRA parameters for scenarios a) and b), so that the phenomenological implications may be much more diverse than in mSUGRA. We show that direct and/or indirect detection of neutralino dark matter should be possible at various current and planned facilities.

We examine the prospects for supersymmetry discovery in the minimal supergravity (mSUGRA) model via indirect detection of neutralino dark matter. We investigate rates for muon detection in neutrino telescopes, and detection of photons, positrons and anti-protons by balloon and space based experiments. We compare the discovery reach in these channels with the reach for direct detection of dark matter, and also with the reach of collider experiments such as Fermilab Tevatron, CERN LHC and a $\sqrt{s}=0.5-1$ TeV linear $e^+e^-$ collider. We pay particular attention to regions of model parameter space in accord with recent WMAP results on the dark matter density of the universe. We find that the IceCube neutrino telescope as well as 3rd generation direct dark matter detection experiments should be able to cover the {\it entire} WMAP allowed portion of the hyperbolic branch/focus point (HB/FP) region of parameter space. This is in contrast to the case of the CERN LHC or a linear $e^+e^-$ collider, where only a fraction of the HB/FP region can be accessed. In addition, we show that detection of $\gamma$s, $e^+$s and $\bar{p}$s should occur in much of the HB/FP region, as well as in the low $m_{1/2}$ portion of the $A$ annihilation funnel, and will be complementary to searches via colliders in these regions.

Rate coefficients for inelastic collisions between Li and H atoms covering all transitions between the asymptotic states Li(2s,2p,3s,3p,3d,4s,4p,4d,4f)+H(1s) and Li^+ +H^- are presented for the temperature range 2000-8000 K based on recent cross-section calculations. The data are of sufficient completeness for non-LTE modelling of the Li I 670.8 nm and 610.4 nm features in late-type stellar atmospheres. Non-LTE radiative transfer calculations in both 1D and 3D model atmospheres have been carried out for test cases of particular interest. Our detailed calculations show that the classical modified Drawin-formula for collisional excitation and de-excitation (Li*+H <-> Li*'+H) over-estimates the cross-sections by typically several orders of magnitude and consequently that these reactions are negligible for the line formation process. However, the charge transfer reactions collisional ion-pair production and mutual neutralization (Li*+H <-> Li^+ +H^-) are of importance in thermalizing Li. In particular, 3D non-LTE calculations of the Li I 670.8 nm line in metal-poor halo stars suggest that 1D non-LTE results over-estimate the Li abundance by up to about 0.1 dex, aggrevating the discrepancy between the observed Li abundances and the primordial Li abundance as inferred by the WMAP analysis of the cosmic microwave background.

We evaluate neutralino-nucleon scattering rates in several well-motivated supersymmetric models, and compare against constraints on the neutralino relic density, BF( b\to s\gamma ) as well as the muon anomalous magnetic moment a_\mu . In the mSUGRA model, the indirect constraints favor the hyperbolic branch/focus point (HB/FP) region of parameter space, and in fact this region is just where neutralino-nucleon scattering rates are high enough to be detected in direct dark matter search experiments! In Yukawa unified SUSY SO(10) models with scalar mass non-universality, the relic density of neutralinos is almost always above experimental bounds, while the corresponding direct detection rates are below experimental levels. Conversely, in five dimensional SO(10) models where gauge symmetry breaking is the result of compactification of the extra dimension, and supersymmetry breaking is communicated via gaugino mediation, the relic density is quite low, while direct detection rates can be substantial.

We evaluate the relic density of neutralinos in the minimal supergravity (mSUGRA) model. All 2-->2 neutralino annihilation diagrams, as well as all processes involving sleptons, charginos, neutralinos and third generation squarks are included. Relativistic thermal averaging of the velocity times cross sections is performed. We find that co-annihilation effects are only important on the edges of the model parameter space, where some amount of fine-tuning is necessary to obtain a reasonable relic density. Alternatively, at high tan(beta), annihilation through very broad Higgs resonances gives rise to an acceptable neutralino relic density over broad regions of parameter space where little or no fine-tuning is needed. Finally, we compare our results against the reach of various e+e- and hadron colliders for supersymmetric matter.