Tilted and warped discs inside tilted dark matter haloes are predicted from numerical and semi-analytical studies. In this paper, we use deep imaging to demonstrate the likely existence of tilted outer structures in real galaxies. We consider two SB0 edge-on galaxies, NGC4469 and NGC4452, which exhibit apparent tilted outer discs with respect to the inner structure. In NGC4469, this structure has a boxy shape, inclined by $\Delta$PA$\approx$3$^{\circ}$ with respect to the inner disc, whereas NGC4452 harbours a discy outer structure with $\Delta$PA$\approx$6$^{\circ}$. In spite of the different shapes, both structures have surface brightness profiles close to exponential and make a large contribution ($\sim30$%) to the total galaxy luminosity. In the case of NGC4452, we propose that its tilted disc likely originates from a former fast tidal encounter (probably with IC3381). For NGC4469, a plausible explanation may also be galaxy harassment, which resulted in a tilted or even a tumbling dark matter halo. A less likely possibility is accretion of gas-rich satellites several Gyr ago. New deep observations may potentially reveal more such galaxies with tilted outer structures, especially in clusters. We also consider galaxies, mentioned in the literature, where a central component (a bar or a bulge) is tilted with respect to the stellar disc. According to our numerical simulations, one of the plausible explanations of such observed "tilts" of the bulge/bar is a projection effect due to a not exactly edge-on orientation of the galaxy coupled with a skew angle of the triaxial bulge/bar.

The HERON project is aimed at studying halos and low surface brightness details near galaxies. In this second HERON paper we consider in detail deep imaging (down to surface brightness of ~28 mag/arcsec$^2$ in the r band) for 35 galaxies, viewed edge-on. We confirm a range of low surface brightness features previously described in the literature but also report new ones. We classify the observed outer shapes of the galaxies into three main types (and their prototypes): disc/diamond-like (NGC891), oval (NGC4302), and boxy (NGC3628). We show that the shape of the outer disc in galaxies does not often follow the general 3D model of an exponential disc: 17 galaxies in our sample exhibit oval or even boxy isophotes at the periphery. Also, we show that the less flattened the outer disc is, the more oval or boxy its structure. Many galaxies in our sample have an asymmetric outer structure. We propose that the observed diversity of the galaxy outer shapes is defined by the merger history and its intensity: if no recent multiple minor or single major merging took place, the outer shape is diamond-like or discy. On the contrary, interacting galaxies show oval outer shapes, whereas recent merging appears to transform the outer shape to boxy.

We report an ab initio study on the spectroscopy of the open-shell diatomic molecule yttrium oxide, YO. The study considers the six lowest doublet states, $X\,{}^{2}\Sigma^{+}$, $A\,{}^{2}\Pi$, $C\,{}^{2}\Pi$, $A'\,{}^{2}\Delta$, $B\,{}^{2}\Sigma^{+}$, $D\,{}^{2}\Sigma^{+}$ and a few higher-lying quartet states using high levels of electronic structure theory and accurate nuclear motion calculations. The coupled cluster singles, doubles, and perturbative triples, CCSD(T), and multireference configuration interaction (MRCI) methods are employed in conjunction with a relativistic pseudopotential on the yttrium atom and a series of correlation-consistent basis sets ranging in size from triple-$\zeta$ to quintuple-$\zeta$ quality. Core-valence correlation effects are taken into account and complete basis set limit extrapolation is performed for CCSD(T). Spin-orbit coupling is included through the use of both MRCI state-interaction with spin-orbit (SI-SO) approach and four-component relativistic equation-of-motion CCSD calculations. Using the ab initio data for bond lengths ranging from 1.0 to 2.5 A, we compute 6 potential energy, 12 spin-orbit, 8 electronic angular momentum, 6 electric dipole moment and 12 transition dipole moment (4 parallel and 8 perpendicular) curves which provide a complete description of the spectroscopy of the system of six lowest doublet states. The Duo nuclear motion program is used to solve the coupled nuclear motion Schr\"{o}dinger equation for these six electronic states. The spectra of $^{89}$Y$^{16}$O simulated for different temperatures are compared with several available high resolution experimental studies; good agreement is found once minor adjustments are made to the electronic excitation energies.

X-structures are often observed in galaxies hosting the so-called B/PS (boxy/peanuts) bulges and are visible from the edge-on view. They are the most notable features of B/PS bulges and appear as four rays protruding from the disk of the host galaxy and distinguishable against the B/PS bulge background. In some works their origin is thought to be connected with the so-called banana-shaped orbits with a vertical resonance 2:1. A star in such an orbit performs two oscillations in the vertical direction per one revolution in the bar frame. Several recent studies that analyzed ensembles of orbits arising in different $N$-body models do not confirm the dominance of the resonant 2:1 orbits in X-structures. In our work we analyze two $N$-body models and show how the X-structure in our models is gradually assembled from the center to the periphery from orbits with less than 2:1 frequency ratio. The most number of such orbits is enclosed in a 'farfalle'-shape (Italian pasta) form and turns out to be non-periodic. We conclude that the X-structure is only the envelope of regions of high density caused by the crossing or folding of different types of orbits at their highest points, and does not have a "backbone" similar to that of the in-plane bar. Comparing the orbital structure of two different numerical models, we show that the dominance of one or another family of orbits with a certain ratio of the vertical oscillations frequency to the in-plane frequency depends on the parameters of the underlying galaxy and ultimately determines the morphology of the X-structure and the opening angle of its rays.

In supernovae neutrinos are emitted from a region with a width $r_{\rm eff}$ of a few kilometers (rather than from a surface of infinitesimal width). We study the effect of integration (averaging) over such an extended emission region on collective oscillations. The averaging leads to additional suppression of the correlation (off-diagonal element of the density matrix) by a factor $ \sim 1/r_{\rm eff} V_e \sim 10^{-10}$, where $V_e$ is the matter potential. This factor enters the initial condition for further collective oscillations and, consequently, leads to a delay of the strong flavour transitions. We justify and quantify this picture using a simple example of collective effects in two intersecting fluxes. We have derived the evolution equation for the density matrix elements integrated over the emission region and solved it both numerically and analytically. For the analytic solution we have used linearized equations. We show that the delay of the development of the instability and the collective oscillations depends on the suppression factor due to the averaging (integration) logarithmically. If the instability develops inside the production region, the integration leads not only to a delay but also to a modification of the exponential grow.

By means of self-consistent numerical simulations we investigated the dynamical impact of classical bulges on the growth of the secondary buckling of a bar. Overall we considered 14 models with different disc and bulge parameters. We obtained that a bulge with a quite modest mass $B/D=0.1$ leads to completely symmetrical evolution of the bar almost independently of the initial stellar disc parameters and even can damp the first bending. At the same time, the bars in all our bulgeless models suffer from the short primary and prolonged secondary buckling. Given the smallness of the mass suppressing secondary buckling, we conclude that a classical bulge along with the gas central concentration may be the main culprits for the rarity of bars with ongoing buckling in the local Universe.

Recent observational studies of X-shaped structures revealed that values of their opening angles lie in a narrow range: from 20 deg to 43 deg with smaller X-shaped structures appearing to follow a characteristic opening angle ~ 27 deg - 31 deg. We use self-consistent numerical simulations to uncover what parameters of host galaxies govern the opening angle spread. We constructed a series of equilibrium models of galaxies with high spatial resolution, varying the dark halo and bulge contribution in the overall gravitational potential, the initial disc thickness of models and the Toomre parameter Q and followed their evolution for almost 8 Gyrs. Each model demonstrated the formation of clear X-structures with different flatness. We have obtained that opening angles lie in the range from 25 deg to 42 deg throughout the entire evolution. These values are roughly consistent with observational data. The greatest variation in the opening angles is obtained by varying the mass of the dark halo. The initial thickness of the disc and the Toomre parameter are responsible for smaller variations of the angle and shorter X-structures. An increase of both parameters changes the morphology of in-plane bars and X-structures. In some cases we observed even double X-structures. The main effect of the bulge is to prevent buckling at intermediate and late stages of the disc evolution. Comparison of models with different halo masses indicates that the smallest observable values of opening angles can be associated with the presence of a heavy dark halo (more than 3 masses of the disc within the optical radius).

The Cluster mission, launched in 2000, has produced a large database of electron flux intensity measurements in the Earths magnetosphere by the Research with Adaptive Particle Imaging Detector (RAPID)/ Imaging Electron Spectrometer (IES) instrument. However, due to background contamination of the data with high-energy electrons (>400 keV) and inner- zone protons (230-630 keV) in the radiation belts and ring current, the data have been rarely used for inner-magnetospheric science. The current paper presents two algorithms for background correction. The first algorithm is based on the empirical contamination percentages by both protons and electrons. The second algorithm uses simultaneous proton observations. The efficiencies of these algorithms are demonstrated by comparison of the corrected Cluster/RAPID/IES data with Van Allen Probes/Magnetic Electron Ion Spectrometer (MagEIS) measurements for 2012-2015. Both techniques improved the IES electron data in the radiation belts and ring current.Yearly averaged flux intensities of the two missions show the ratio of measurements close to 1. IES corrected measurements were also compared with Arase Medium-Energy Particle Experiments-Electron Analyzer (MEP-e) electron data during two conjunction periods in 2017 and also exhibited ratio close to 1. We demonstrate a scientific application of the corrected IES electron data analyzing its evolution during solar cycle. Spin-averaged yearly mean IES electron intensities in the outer belt for energies 40-400 keV at L-shell between 4 and 6 showed high positive correlation with AE index and solar wind dynamic pressure during 2001- 2016. Relationship between solar wind dynamic pressure and IES electron measurements in the outer radiation belt was derived as a uniform linear-logarithmic equation.

Building on previous results on the quadratic helicity in magnetohydrodynamics (MHD) we investigate particular minimum helicity states. Those are eigenfunctions of the curl operator and are shown to constitute solutions of the quasi-stationary incompressible ideal MHD equations. We then show that these states have indeed minimum quadratic helicity.

Collective oscillations of supernova neutrinos above the neutrino sphere can be completely described by the propagation of individual neutrinos in external potentials and are in this sense a linear phenomenon. An effective theory of collective oscillations can be developed based on certain assumptions about time dependence of these potentials. General conditions for strong flavor transformations are formulated and these transformations can be interpreted as parametric resonance effects induced by periodic modulations of the potentials. We study a simplified and solvable example, where a probe neutrino is propagating in a flux of collinear neutrinos, such that $\nu \nu-$ interactions in the flux are absent. Still, this example retains the main feature - the coherent flavor exchange. Properties of the parametric resonance are studied, and it is shown that integrations over energies and emission points of the flux neutrinos suppress modulations of the potentials and therefore strong transformations. The transformations are also suppressed by changes in densities of background neutrinos and electrons.

We present an approach to testing the gravitational redshift effect using the RadioAstron satellite. The experiment is based on a modification of the Gravity Probe A scheme of nonrelativistic Doppler compensation and benefits from the highly eccentric orbit and ultra-stable atomic hydrogen maser frequency standard of the RadioAstron satellite. Using the presented techniques we expect to reach an accuracy of the gravitational redshift test of order $10^{-5}$, a magnitude better than that of Gravity Probe A. Data processing is ongoing, our preliminary results agree with the validity of the Einstein Equivalence Principle.

The Nobel prize in physics 2015 has been awarded "... for the discovery of neutrino oscillations which show that neutrinos have mass". While SuperKamiokande (SK), indeed, has discovered oscillations, SNO observed effect of the adiabatic (almost non-oscillatory) flavor conversion of neutrinos in the matter of the Sun. Oscillations are irrelevant for solar neutrinos apart from small $\nu_e$ regeneration inside the Earth. Both oscillations and adiabatic conversion do not imply masses uniquely and further studies were required to show that non-zero neutrino masses are behind the SNO results. Phenomena of oscillations (phase effect) and adiabatic conversion (the MSW effect driven by the change of mixing in matter) are described in pedagogical way.

Solar neutrino studies triggered and largely motivated the major developments in neutrino physics in the last 50 years. Theory of neutrino propagation in different media with matter and fields has been elaborated. It includes oscillations in vacuum and matter, resonance flavor conversion and resonance oscillations, spin and spin-flavor precession, etc. LMA MSW has been established as the true solution of the solar neutrino problem. Parameters theta12 and Delta_m21^2 have been measured; theta13 extracted from the solar data is in agreement with results from reactor experiments. Solar neutrino studies provide a sensitive way to test theory of neutrino oscillations and conversion. Characterized by long baseline, huge fluxes and low energies they are a powerful set-up to search for new physics beyond the standard 3nu paradigm: new neutrino states, sterile neutrinos, non-standard neutrino interactions, effects of violation of fundamental symmetries, new dynamics of neutrino propagation, probes of space and time. These searches allow us to get stringent, and in some cases unique bounds on new physics. We summarize the results on physics of propagation, neutrino properties and physics beyond the standard model obtained from studies of solar neutrinos.

Supernova neutrinos have several exceptional features which can lead to interesting physical consequences. At the production point their wave packets have an extremely small size $\sigma_x \sim 10^{-11}$ cm; hence the energy uncertainty can be as large as the energy itself, $\sigma_E \sim E$, and the coherence length is short. On the way to the Earth the wave packets of mass eigenstates spread to macroscopic sizes and separate. Inside the Earth the mass eigenstates split into eigenstates in matter and oscillate again. The coherence length in the Earth is comparable with the radius of the Earth. We explore these features and their consequences. (i) We present new estimates of the wave packet size. (ii) We consider the decoherence condition for the case of wave packets with spatial spread and show that it is not modified by the spread. (iii) We study the coherence of neutrinos propagating in a multi-layer medium with density jumps at the borders of layers. In this case coherence can be partially restored due to a "catch-up effect", increasing the coherence length beyond the usual estimate. This catch-up effect can occur for supernova neutrinos as they cross the shock wave fronts in the exploding star or the core of the Earth.

We consider several aspects related to the form, derivation and applications of the Liouville equation (LE) for flavour evolution of neutrinos. To take into account the quantum nature of neutrinos we derive the evolution equation for the matrix of densities using wave packets instead of Wigner functions. The obtained equation differs from the standard LE by an additional term which is proportional to the difference of group velocities. We show that this term describes loss of the propagation coherence in the system. In absence of inelastic collisions, the LE can be reduced to a single derivative equation over a trajectory coordinate. Additional time and spacial dependence may steam from initial (production) conditions. The transition from single neutrino evolution to the evolution of a neutrino gas is considered.

A test of a cornerstone of general relativity, the gravitational redshift effect, is currently being conducted with the RadioAstron spacecraft, which is on a highly eccentric orbit around Earth. Using ground radio telescopes to record the spacecraft signal, synchronized to its ultra-stable on-board H-maser, we can probe the varying flow of time on board with unprecedented accuracy. The observations performed so far, currently being analyzed, have already allowed us to measure the effect with a relative accuracy of $4\times10^{-4}$. We expect to reach $2.5\times10^{-5}$ with additional observations in 2016, an improvement of almost a magnitude over the 40-year old result of the GP-A mission.

We study hypermagnetic helicity and lepton asymmetry evolution in plasma of the early Universe before the electroweak phase transition (EWPT) accounting for chirality flip processes via inverse Higgs decays and sphaleron transitions which violate the left lepton number and wash out the baryon asymmetry of the Universe (BAU). In the scenario where the right electron asymmetry supports the BAU alone through the conservation law $B/3 - L_{eR}=const$ at temperatures $T>T_{RL}\simeq 10~TeV$ the following universe cooling leads to the production of a non-zero left lepton (electrons and neutrinos) asymmetry. This is due to the Higgs decays becoming more faster when entering the equilibrium at $T=T_{RL}$ with the universe expansion, $\Gamma_{RL}\sim T> H\sim T^2$ , resulting in the parallel evolution of the right and the left electron asymmetries at $T<T_{RL}$ through the corresponding Abelian anomalies in SM in the presence of a seed hypermagnetic field. The hypermagnetic helicity evolution proceeds in a self-consistent way with the lepton asymmetry growth. The role of sphaleron transitions decreasing the left lepton number turns out to be negligible in given scenario. The hypermagnetic helicity plays a key role in lepto/baryogenesis in our scenario and the more hypermagnetic field is close to the maximum helical one the faster BAU grows up the observable value , $B_{obs}\sim 10^{-10}$.

We investigate the evolution of the baryon asymmetry of the Universe (BAU) in its symmetric phase before the electroweak phase transition (EWPT) induced by leptogenesis in the hypermagnetic field of an arbitrary structure and with a maximum hypermagnetic helicity density. The novelty of this work is that the BAU has been calculated for a continuous hypermagnetic helicity spectrum. The observed BAU $B_{obs} = 10^{-10}$ that can be in large-scale hypermagnetic fields satisfying the wave number inequality $k \leq k_{max}$ grows with increasing $k_{max}$. We will also show that the initial right-handed electron asymmetry $\xi_{eR}(\eta_0)$ used in our leptogenesis model as a free parameter cannot take too large values, $\xi_{eR}(\eta_0) = 10^{-4}$, because this leads to a negative BAU by the EWPT time. In contrast, a sufficiently small initial right-handed electron asymmetry, $\xi_{eR}(\eta_0)$, provides its further growth and the corresponding BAU growth from zero to some positive value,including the observed $B_{obs} = 10^{-10}$.

This paper describes outstanding issues in astrophysics and cosmology that can be solved by astronomical observations in a broad spectral range from far infrared to millimeter wavelengths. The discussed problems related to the formation of stars and planets, galaxies and the interstellar medium, studies of black holes and the development of the cosmological model can be addressed by the planned space observatory Millimetron (the "Spectr-M" project) equipped with a cooled 10-m mirror. Millimetron can operate both as a single-dish telescope and as a part of a space-ground interferometer with very long baseline.

We study oscillations of the high energy atmospheric neutrinos in the Earth into sterile neutrinos with the eV-scale mass. The MSW resonance and parametric enhancement of the $\bar{\nu}_\mu\to\bar{\nu}_s$ oscillations lead to distortion of the zenith angle distribution of the muon-track events which can be observed by IceCube. Due to matter effect, the IceCube signal depends not only on the mixing element $U_{\mu 4}$ relevant for LSND and MiniBooNE but also on $U_{\tau 4}$ and the CP-violating phase $\delta_{24}$. We show that the case with $U_{\tau 4} = \delta_{24} = 0$ leads to the weakest IceCube signal and therefore should be used to bound $U_{\mu 4}$. We compute the zenith angle distributions of the $\nu_\mu-$events for different energy intervals in the range (0.1 - 10) TeV and find that inclusion of the energy information (binning in energy) improves the sensitivity to $\nu_s$ drastically. We estimate that with already collected (during 3 - 4 years) IceCube statistics the bound $|U_{\mu 4}|^2 < 0.01$ ($99\%$ C.L.) can be established and the mixing required by LSND and MiniBooNE can be excluded at $(4 - 6) \sigma$ confidence level.

We study hypermagnetic helicity and lepton asymmetry evolution in plasma of the early Universe before the electroweak phase transition (EWPT) accounting for chirality flip processes via inverse Higgs decays and sphaleron transitions which violate the left lepton number and wash out the baryon asymmetry of the Universe (BAU). In the scenario where the right electron asymmetry supports the BAU alone through the conservation law $B/3 - L_{eR}=const$ at temperatures $T>T_{RL}\simeq 10 TeV$ the following universe cooling leads to the production of a non-zero left lepton (electrons and neutrinos) asymmetry. This is due to the Higgs decays becoming more faster when entering the equilibrium at $T=T_{RL}$ with the universe expansion, $\Gamma_{RL}\sim T> H\sim T^2$, resulting in the parallel evolution of both the right and the left electron asymmetries at $T<T_{RL}$ through the corresponding Abelian anomalies in SM in the presence of a seed hypermagnetic field. The hypermagnetic helicity evolution proceeds in a self-consistent way with the lepton asymmetry growth. The role of sphaleron transitions decreasing the left lepton number turns out to be negligible in given scenario. The hypermagnetic helicity can be a supply for the magnetic one in Higgs phase assuming a strong seed hypermagnetic field in symmetric phase.

We explore the physics potential of multi-megaton scale ice or water Cherenkov detectors with low ($\sim 1$ GeV) threshold. Using some proposed characteristics of the PINGU detector setup we compute the distributions of events versus neutrino energy $E_\nu$ and zenith angle $\theta_z$, and study their dependence on yet unknown neutrino parameters. The $(E_\nu - \theta_z)$ regions are identified where the distributions have the highest sensitivity to the neutrino mass hierarchy, to the deviation of the 2-3 mixing from the maximal one and to the CP-phase. We evaluate significance of the measurements of the neutrino parameters and explore dependence of this significance on the accuracy of reconstruction of the neutrino energy and direction. The effect of degeneracy of the parameters on the sensitivities is also discussed. We estimate the characteristics of future detectors (energy and angle resolution, volume, etc.) required for establishing the neutrino mass hierarchy with high confidence level. We find that the hierarchy can be identified at $3\sigma$ -- $10\sigma$ level (depending on the reconstruction accuracies) after 5 years of PINGU operation.

We show that study of the atmospheric neutrinos in the 10--100 GeV energy range by DeepCore sub-array of the IceCube Neutrino Observatory can substantially constrain the mixing of sterile neutrinos of mass $\sim 1$ eV with active neutrinos. In the scheme with one sterile neutrino we calculate $\nu_\mu-$ and $\bar{\nu}_\mu-$ oscillation probabilities as well as zenith angle distributions of $\nu_\mu^{CC}$ events in different energy intervals in DeepCore. The distributions depend on the mass hierarchy of active neutrinos. Therefore, in principle, the hierarchy can be identified, if $\nu_s$ exists. After a few years of exposure the DeepCore data will allow to exclude the mixing $|U_{\mu 4}|^2 \geq 0.02$ indicated by the LSND/MiniBooNE results. Combination of the DeepCore and high energy IceCube data will further improve sensitivity to $\nu_s$ mixing parameters.

Neutrino oscillations are only observable when the neutrino production, propagation and detection coherence conditions are satisfied. In this paper we consider in detail neutrino production coherence, taking \pi\to \mu \nu \ decay as an example. We compare the oscillation probabilities obtained in two different ways: (1) coherent summation of the amplitudes of neutrino production at different points along the trajectory of the parent pion; (2) averaging of the standard oscillation probability over the neutrino production coordinate in the source. We demonstrate that the results of these two different approaches exactly coincide, provided that the parent pion is considered as pointlike and the detection process is perfectly localized. In this case the standard averaging of the oscillation probability over the finite spatial extensions of the neutrino source (and detector) properly takes possible decoherence effects into account. We analyze the reason for this equivalence of the two approaches and demonstrate that for pion wave packets of finite width \sigma_{x\pi} the equivalence is broken. The leading order correction to the oscillation probability due to \sigma_{x\pi}\ne 0 is shown to be \sim [v_g/(v_g-v_\pi)]\sigma_{x\pi}/l_{osc}, where v_g and v_\pi \ are the group velocities of the neutrino and pion wave packets, and l_{osc} is the neutrino oscillation length.

Oscillation interpretation of the results from the LSND, MiniBooNE and some other experiments requires existence of sterile neutrino with mass $\sim 1$ eV and mixing with the active neutrinos $|U_{\mu 0}|^2 \sim (0.02 - 0.04)$. It has been realized some time ago that existence of such a neutrino affects significantly the fluxes of atmospheric neutrinos in the TeV range which can be tested by the IceCube Neutrino Observatory. In view of the first IceCube data release we have revisited the oscillations of high energy atmospheric neutrinos in the presence of one sterile neutrino. Properties of the oscillation probabilities are studied in details for various mixing schemes both analytically and numerically. The energy spectra and angular distributions of the $\nu_\mu-$events have been computed for the simplest $\nu_s-$mass, and $\nu_s - \nu_\mu$ mixing schemes and confronted with the IceCube data. An illustrative statistical analysis of the present data shows that in the $\nu_s-$mass mixing case the sterile neutrinos with parameters required by LSND/MiniBooNE can be excluded at about $3\sigma$ level. The $\nu_s- \nu_\mu$ mixing scheme, however, can not be ruled out with currently available IceCube data.

Recent results from the SNO, Super-Kamiokande and Borexino experiments do not show the expected upturn of the energy spectrum of events (the ratio $R \equiv N_{obs}/N_{SSM}$) at low energies. At the same time, cosmological observations testify for possible existence of additional relativistic degrees of freedom in the early Universe: $\Delta N_{eff} = 1 - 2$. These facts strengthen the case of very light sterile neutrino, $\nu_s$, with $\Delta m^2_{01} \sim (0.7 - 2) \cdot 10^{-5}$ eV$^2$, which mixes weakly with the active neutrinos. The $\nu_s$ mixing in the mass eigenstate $\nu_1$ characterized by $\sin^2 2\alpha \sim 10^{-3}$ can explain an absence of the upturn. The mixing of $\nu_s$ in the eigenstate $\nu_3$ with $\sin^2 \beta \sim 0.1$ leads to production of $\nu_s$ via oscillations in the Universe and to additional contribution $\Delta N_{eff} \approx 0.7 - 1$ before the big bang nucleosynthesis and later. Such a mixing can be tested in forthcoming experiments with the atmospheric neutrinos as well as in future accelerator long baseline experiments. It has substantial impact on conversion of the supernova neutrinos.

High energy cosmic neutrino fluxes can be produced inside relativistic jets under the envelopes of collapsing stars. In the energy range E ~ (0.3 - 1e5) GeV, flavor conversion of these neutrinos is modified by various matter effects inside the star and the Earth. We present a comprehensive (both analytic and numerical) description of the flavor conversion of these neutrinos which includes: (i) oscillations inside jets, (ii) flavor-to-mass state transitions in an envelope, (iii) loss of coherence on the way to observer, and (iv) oscillations of the mass states inside the Earth. We show that conversion has several new features which are not realized in other objects, in particular interference effects ("L- and H- wiggles") induced by the adiabaticity violation. The neutrino-neutrino scattering inside jet and inelastic neutrino interactions in the envelope may produce some additional features at E > 1e4 GeV. We study dependence of the probabilities and flavor ratios in the matter-affected region on angles theta13 and theta23, on the CP-phase delta, as well as on the initial flavor content and density profile of the star. We show that measurements of the energy dependence of the flavor ratios will, in principle, allow to determine independently the neutrino and astrophysical parameters.

Collective oscillations of supernova neutrinos swap the electron neutrino and antineutrino spectra with those of another flavor in certain energy intervals bounded by sharp spectral splits. This phenomenon is far more general than previously appreciated: typically one finds one or more swaps and accompanying splits in the neutrino and antineutrino channels for both inverted and normal neutrino mass hierarchies. Depending on an instability condition, swaps develop around spectral crossings (energies where the electron neutrino or antineutrino fluxes are equal to that of another flavor, as well as infinite E where all fluxes vanish), and the widths of swaps are determined by the spectra and fluxes. Wash-out by multi-angle decoherence varies across the spectrum and splits can survive as sharp spectral features.

Despite the theory of neutrino oscillations being rather old, some of its basic issues are still being debated in the literature. We discuss, in the framework of the wave packet approach, a number of such issues, including the relevance of the "same energy" and "same momentum" assumptions, the role of quantum-mechanical uncertainty relations in neutrino oscillations, the dependence of the production/detection and propagation coherence conditions that ensure the observability of neutrino oscillations on neutrino energy and momentum uncertainties, the question of (in)dependence of the oscillation probabilities on the neutrino production and detection processes, the applicability limits of the stationary source approximation, and Lorentz invariance of the oscillation probability. We also develop a novel approach to calculation of the oscillation probability in the wave packet picture, based on the summation/integration conventions different from the standard one, which gives a new insight into the oscillation phenomenology. We discuss a number of apparently paradoxical features of the theory of neutrino oscillations.

There are several new features in production, oscillations and detection of the atmospheric neutrinos of low energies, E<100 MeV. The flavor ratio, r, of muon to electron neutrino fluxes is substantially smaller than 2 and decreases with energy, significant part of events is due to the decay of invisible muons at rest, etc. Oscillations in two-layer medium (atmosphere - earth) should be taken into account. We derive analytical and semi-analytical expressions for the oscillation probabilities of these "sub-subGeV" neutrinos. The energy spectra of the e-like events in water cherenkov detectors are computed and dependence of the spectra on the 2-3 mixing angle, the 1-3 mixing and CP-violation phase are studied. We find that variations of 2-3 mixing angle in the presently allowed region change the number of e-like events by about 15 - 20 % as well as to distortion of the energy spectrum. The 1-3 mixing and CP-violation can lead to ~10% effects. Detailed study of the sub-subGeV neutrinos will be possible in future Megaton-scale detectors.

We develop a comprehensive description of three flavor neutrino oscillations inside the Earth in terms of neutrino oscillograms in the whole range of nadir angles and for energies above 0.1 GeV. The effects of the 1-2 mass splitting and mixing as well the interference of the 1-2 and 1-3 modes of oscillations are quantified. The 1-2 mass splitting and mixing lead to the appearance, apart from the resonance MSW peaks, of the parametric resonance peak for core-crossing trajectories at E_nu ~= 0.2 GeV. We show that the interference effects, in particular CP violation, have a domain structure with borders determined by the solar and atmospheric magic lines and the lines of the interference phase condition. The dependence of the oscillograms on the Dirac CP-violating phase is studied. We show that for sin^2(2 theta13) < 0.1 the strongest dependence of the oscillograms on delta is in the 1-2 and 1-3 resonance regions.

In the dense-neutrino region above the neutrino sphere of a supernova (r > 400 km), neutrino-neutrino refraction causes collective flavor transformations. They can lead to "spectral splits" where an energy E_split splits the transformed spectrum sharply into parts of almost pure but different flavors. Unless there is an ordinary MSW resonance in the dense-neutrino region, E_split is determined by flavor-lepton number conservation alone. Spectral splits are created by an adiabatic transition between regions of large and small neutrino density. We solve the equations of motion in the adiabatic limit explicitly and provide analytic expressions for a generic example.

We develop a detailed and comprehensive description of neutrino oscillations driven by the 1-3 mixing in the matter of the Earth. The description is valid for the realistic (PREM) Earth density profile in the whole range of nadir angles and for neutrino energies above 1 GeV. It can be applied to oscillations of atmospheric and accelerator neutrinos. The results are presented in the form of neutrino oscillograms of the Earth, i.e. the contours of equal oscillation probabilities in the neutrino energy--nadir angle plane. A detailed physics interpretation of the oscillograms, which includes the MSW peaks, parametric ridges, local maxima, zeros and saddle points, is given in terms of the amplitude and phase conditions. Precise analytic formulas for the probabilities are obtained. We study the dependence of the oscillation pattern on theta_13 and find, in particular, that the transition probability P > 1/2 appears for sin^2(2theta_13) as small as 0.009. We consider the dependence of the oscillation pattern on the matter density profile and comment on the possibility of the oscillation tomography of the Earth.

We present a formalism for precise description of oscillation phenomena in matter at high energies or high densities, V > \Delta m^2/2E, where V is the matter-induced potential of neutrinos. The accuracy of the approximation is determined by the quantity \sin^2 2\theta_m \Delta V/2\pi V, where \theta_m is the mixing angle in matter and \Delta V is a typical change of the potential over the oscillation length (l \sim 2\pi/V). We derive simple and physically transparent formulas for the oscillation probabilities, which are valid for arbitrary matter density profiles. They can be applied to oscillations of high energy (E > 10 GeV) accelerator, atmospheric and cosmic neutrinos in the matter of the Earth, substantially simplifying numerical calculations and providing an insight into the physics of neutrino oscillations in matter. The effect of parametric enhancement of the oscillations of high energy neutrinos is considered. Future high statistics experiments can provide an unambiguous evidence for this effect.

Neutrinos may possibly violate the spin-statistics theorem, and hence obey Bose statistics or mixed statistics despite having spin half. We find the generalized equilibrium distribution function of neutrinos which depends on a single fermi-bose parameter, \kappa, and interpolates continuously between the bosonic and fermionic distributions when \kappa changes from -1 to +1. We consider modification of the Big Bang Nucleosynthesis (BBN) in the presence of bosonic or partly bosonic neutrinos. For pure bosonic neutrinos the abundances change (in comparison with the usual Fermi-Dirac case) by -3.2% for 4He (which is equivalent to a decrease of the effective number of neutrinos by \Delta N_\nu = - 0.6), +2.6% for 2H and -7% for 7Li. These changes provide a better fit to the BBN data. Future BBN studies will be able to constrain the fermi-bose parameter to \kappa > 0.5, if no deviation from fermionic nature of neutrinos is found. We also evaluate the sensitivity of future CMB and LSS observations to the fermi-bose parameter.

We assume that the Pauli exclusion principle is violated for neutrinos, and consequently, neutrinos obey the Bose-Einstein statistics. Cosmological and astrophysical consequences of this assumption are considered. Neutrinos may form cosmological Bose condensate which accounts for all (or a part of) the dark matter in the universe. ``Wrong'' statistics of neutrinos could modify big bang nucleosynthesis, leading to the effective number of neutrino species smaller than three. Dynamics of the supernova collapse would be influenced and spectra of the supernova neutrinos may change. The presence of neutrino condensate would enhance contributions of the Z-bursts to the flux of the UHE cosmic rays and lead to substantial refraction effects for neutrinos from remote sources. The Pauli principle violation for neutrinos can be tested in the two-neutrino double beta decay.

Solar neutrino physics enters a stage of precision measurements. In this connection we present a precise analytic description of the neutrino conversion in the context of LMA MSW solution of the solar neutrino problem. Using the adiabatic perturbation theory we derive an analytic formula for the $\nu_e$ survival probability which takes into account the non-adiabatic corrections and the regeneration effect inside the Earth. The probability is averaged over the neutrino production region. We find that the non-adiabatic corrections are of the order $10^{-9}-10^{-7}$. Using the formula for the Earth regeneration effect we discuss features of the zenith angle dependence of the $\nu_e$ flux. In particular, we show that effects of small structures at the surface of the Earth can be important.

The Homestake result is about ~ 2 \sigma lower than the Ar-production rate, Q_{Ar}, predicted by the LMA MSW solution of the solar neutrino problem. Also there is no apparent upturn of the energy spectrum (R \equiv N_{obs}/N_{SSM}) at low energies in SNO and Super-Kamiokande. Both these facts can be explained if a light, \Delta m^2_{01} ~ (0.2 - 2) \cdot 10^{-5} eV^2, sterile neutrino exists which mixes very weakly with active neutrinos: \sin^2 2\alpha ~ (10^{-5} - 10^{-3}). We perform both the analytical and numerical study of the conversion effects in the system of two active neutrinos with the LMA parameters and one weakly mixed sterile neutrino. The presence of sterile neutrino leads to a dip in the survival probability in the intermediate energy range E = (0.5 - 5) MeV thus suppressing the Be, or/and pep, CNO as well as B electron neutrino fluxes. Apart from diminishing Q_{Ar} it leads to decrease of the Ge-production rate and may lead to decrease of the BOREXINO signal and CC/NC ratio at SNO. Future studies of the solar neutrinos by SNO, SK, BOREXINO and KamLAND as well as by the new low energy experiments will allow us to check this possibility. We present a general analysis of modifications of the LMA energy profile due to mixing with new neutrino states.

After recent results from solar neutrino experiments and KamLAND we know that neutrinos from SN1987A underwent flavor conversion. This must be taken into account in the analysis of the data. Assuming the normal mass hierarchy of neutrinos we calculate the permutation factors p for the Kamiokande-2, IMB and Baksan detectors. The conversion inside the star leads to p = 0.28 - 0.32; regeneration effects inside the Earth reduce this factor down to 0.15 - 0.20 (at E = 40 MeV). We study in details the influence of conversion on the observed signal depending on the parameters of the original neutrino spectra. For a given set of these parameters, the conversion could increase the average energy of the observed events by up to 50 % and the number of events by a factor of 2 at Kamiokande-2 and of 3 - 5 at IMB. Inversely, we find that neglecting the conversion effects can lead up to 50 % error in the determination of the average energy of the original antinu_e flux and about 50 % error on its luminosity. By comparison with experimental data we conclude that the Kamiokande-2 data alone do not favor strong conversion effect. In contrast, the combined analysis of the Kamiokande and IMB results slightly favors strong conversion effects (that is, large difference of the original spectra). In comparison with the no oscillation case, the latter require lower average energy and higher luminosity of the original antinu_e flux.

We have performed analysis of the solar neutrino data including results from the SNO salt phase as well as the combined analysis of the solar and the KamLAND results. The best fit values of neutrino parameters are Delta m^2 = 7.1e-5 eV^2, tan^2\theta = 0.40 with the boron flux f_B = 1.04. New SNO results strongly disfavor maximal mixing and the h-LMA region (Delta m^2 > 1e-4 eV^2) which is accepted now at the 3-sigma level. We find the 3-sigma upper bounds: Delta m^2 < 1.7e-4$ eV^2 and tan^2\theta < 0.64, and the lower bound Delta m^2 > 4.8e-5 eV^2. Non-zero 13-mixing does not change these results significantly. The present data determine quantitatively the physical picture of the solar neutrino conversion. At high energies relevant for SNO and Super-Kamiokande the deviation of the effective survival probability from the non-oscillatory value is about 10 - 14%. The oscillation effect contribution to this difference about 10% and the Earth regeneration is about 3 - 4%. At low energies (E < 1 MeV) the matter corrections to vacuum oscillation effect are below 5%. The predictions for the forthcoming measurements are given which include the spectral distortion and CC/NC ratio at SNO, the Day-Night asymmetry, the KamLAND spectrum and rate.

In the framework of the seesaw mechanism with three heavy right-handed Majorana neutrinos and no Higgs triplets we carry out a systematic study of the structure of the right-handed neutrino sector. Using the current low-energy neutrino data as an input and assuming hierarchical Dirac-type neutrino masses $m_{Di}$, we calculate the masses $M_i$ and the mixing of the heavy neutrinos. We confront the inferred properties of these neutrinos with the constraints coming from the requirement of a successful baryogenesis via leptogenesis. In the generic case the masses of the right-handed neutrinos are highly hierarchical: $M_i \propto m_{Di}^2$; the lightest mass is $M_1 \approx 10^3 - 10^6$ GeV and the generated baryon-to-photon ratio $\eta_B\lesssim 10^{-14}$ is much smaller than the observed value. We find the special cases which correspond to the level crossing points, with maximal mixing between two quasi-degenerate right-handed neutrinos. Two level crossing conditions are obtained: ${m}_{ee}\approx 0$ (1-2 crossing) and $d_{12}\approx 0$ (2-3 crossing), where ${m}_{ee}$ and $d_{12}$ are respectively the 11-entry and the 12-subdeterminant of the light neutrino mass matrix in the basis where the neutrino Yukawa couplings are diagonal. We show that sufficient lepton asymmetry can be produced only in the 1-2 crossing where $M_1 \approx M_2 \approx 10^{8}$ GeV, $M_3 \approx 10^{14}$ GeV and $(M_2 - M_1)/ M_2 \lesssim 10^{-5}$.

We consider in details the effects of the 13-mixing (sin^2 theta_{13}) and of the type of mass hierarchy/ordering (sign[ Delta m^2_{13}]) on neutrino signals from the gravitational collapses of stars. The observables (characteristics of the energy spectra of nu_e and antinu_e events) sensitive to sin^2 theta_{13} and sign[Delta m^2_{13}] have been calculated. They include the ratio of average energies of the spectra, r_E = <E>/<\bar E >, the ratio of widths of the energy distributions, r_Gamma, the ratios of total numbers of nu_e and antinu_e events at low energies, S, and in the high energy tails, R_{tail}. We construct and analyze scatter plots which show the predictions for the observables for different intervals of sin^2 theta_{13} and signs of Delta m^2_{13}, taking into account uncertainties in the original neutrino spectra, the star density profile, etc.. Regions in the space of observables r_E, r_Gamma, S, R_{tail} exist in which certain mass hierarchy and intervals of sin^2 theta_{13} can be identified or discriminated. We elaborate on the method of the high energy tails in the spectra of events. The conditions are formulated for which sin^2 theta_{13} can be (i) measured, (ii) restricted from below, (iii) restricted from above. We comment on the possibility to determine sin^2 theta_{13} using the time dependence of the signals due to the propagation of the shock wave through the resonance layers of the star. We show that the appearance of the delayed Earth matter effect in one of the channels (nu_e or antinu_e) in combination with the undelayed effect in the other channel will allow to identify the shock wave appeareance and determine the mass hierarchy.

The first KamLAND results are in a very good agreement with the predictions made on the basis of the solar neutrino data and the LMA realization of the MSW mechanism. We perform a combined analysis of the KamLAND (rate, spectrum) and the solar neutrino data with a free boron neutrino flux f_B. The best fit values of neutrino parameters are Delta m^2 = 7.3e-5 eV^2, tg^2 theta = 0.41 and f_B = 1.05 with the 1 sigma intervals: Delta m^2 = (6.2 - 8.4)e-5 eV^2, tg^2 theta = 0.33 - 0.54. We find the 3 sigma upper bounds: Delta m^2 < 2.8e-4 eV^2 and tg^2 theta < 0.84, and the lower bound Delta m^2 > 4e-5 eV^2. At 99% C.L. the KamLAND spectral result splits the LMA region into two parts with the preferred one at Delta m^2 < 1e-4 eV^2. The higher Delta m^2 region is accepted at about 2 sigma level. We show that effects of non-zero 13-mixing, sin^2 theta_{13} < 0.04, are small leading to a slight improvement of the fit in higher Delta m^2 region. In the best fit point we predict for SNO: CC/NC = 0.33 +0.05-0.03 and A_{DN}(SNO) = 2.8+-0.8 % (68% C.L.), and A_{DN}(SNO) < 9 % at the 3 sigma level. Further improvements in the determination of the parameters are discussed and implications of the solar neutrino and KamLAND results are considered.

The formalism of flavor conversion of supernova neutrinos is generalized to include possible differences in the fluxes of the muon and tau neutrinos produced in the star. In this case the radiatively induced difference of the nu_mu and nu_tau potentials in matter becomes important. The nu_mu and nu_tau flux differences can manifest themselves in the effects of the Earth matter on the observed nu_e (antinu_e) signal if: (i) the neutrino mass hierarchy is normal (inverted); (ii) the solution of the solar neutrino problem is in the LMA region; (iii) the mixing U_{e3} is relatively large: |U_{e3}|>10^{-3}. We find that for differences in the nu_mu - nu_tau (antinu_mu - antinu_tau) average energies and/or integrated luminosities < 20 %, the relative deviation of the observed nu_e (antinu_e) energy spectrum at E > 50 MeV from that in the case of the equal fluxes can reach 20 - 30 % (10 - 15 %) for neutrinos crossing the Earth. It could be detected in future if large detectors sensitive to the nu_e (antinu_e) energy spectrum become available. The study of this effect would allow one to test the predictions of the nu_mu, nu_tau, antinu_mu, antinu_tau fluxes from supernova models and therefore give an important insight into the properties of matter at extreme conditions. It should be taken into account in the reconstruction of the neutrino mass spectrum and mixing matrix from the supernova neutrino observations. We show that even for unequal nu_mu and nu_tau fluxes, effects of leptonic CP violation can not be studied in the supernova neutrino experiments.

We perform a detailed study of the Earth matter effects on supernova neutrinos. The dependences of these effects on the properties of the original neutrino fluxes, on the trajectory of the neutrinos inside the Earth and on the oscillation parameters are described. We show that, for a large fraction (60 %) of the possible arrival times of the signal, the neutrino flux crosses a substantial amount of the matter of the Earth at least for one of the existing detectors. For oscillation parameters from the LMA solution of the solar neutrino problem the Earth matter effect consists in an oscillatory modulation of the $\barnue$ and/or $\nu_e$ energy spectra. The relative deviation with respect to the undistorted spectra can be as large as 20-30 % for E> 20 MeV and 70-100 % for E> 40 MeV. For parameters from the SMA and LOW solutions the effect is localized at low energies (E<10 MeV) and is not larger than 10 %. The Earth matter effects can be revealed (i) by the observation of oscillatory distortions of the energy spectra in a single experiment and (ii) by the comparison between the spectra at different detectors. For a supernova at distance D=10 Kpc, comparing the results of SuperKamiokande (SK), SNO and LVD experiments one can establish the effect at (2-3) \sigma level, whereas larger statistical significance ((4-5) \sigma) is obtained if two experiments of SK-size or larger are available. Studies of the Earth matter effect will select or confirm the solution of the solar neutrino problem, probe the mixing U_{e 3} and identify the hierarchy of the neutrino mass spectrum.

We study matter effects on oscillations of high-energy neutrinos in the Universe. Substantial effect can be produced by scattering of the neutrinos from cosmological sources ($z\gta 1$) on the relic neutrino background, provided that the latter has large CP-asymmetry: $\eta\equiv (n_\nu-n_{\bar{\nu}})/n_\gamma\gta 1$, where $n_\nu$, $n_{\bar{\nu}}$ and $n_\gamma$ are the concentrations of neutrinos, antineutrinos and photons. We consider in details the dynamics of conversion in the expanding neutrino background. Applications are given to the diffuse fluxes of neutrinos from GRBs, AGN, and the decay of super-heavy relics. We find that the vacuum oscillation probability can be modified by $\sim (10-20)%$ and in extreme cases allowed by present bounds on $\eta$ the effect can reach $\sim 100%$. Signatures of matter effects would consist (i) for both active-active and active-sterile conversion, in a deviation of the numbers of events produced in a detector by neutrinos of different flavours, $N_{\alpha}~(\alpha=e,\mu,\tau)$, and of their ratios from the values given by vacuum oscillations; such deviations can reach $\sim 5-15%$, (ii) for active-sterile conversion, in a characteristic energy dependence of the ratios $N_{e}/N_{\mu},N_{e}/N_{\tau},N_{\mu}/N_{\tau}$. Searches for these matter effects will probe large CP and lepton asymmetries in the universe.

We present a global analysis of all the available solar neutrino data treating consistently the 8B and hep neutrino fluxes as free parameters. The analysis reveals at 99.7% C.L. eight currently-allowed discrete regions in two-neutrino oscillation space, five regions corresponding to active neutrinos and three corresponding to sterile neutrinos. Most of the allowed solutions are robust with respect to changes in the analysis procedure, but the traditional vacuum solution is fragile. The globally-permitted range of the 8B neutrino flux, 0.45 to 1.95 in units of the BP2000 flux, is comparable to the 3 sigma range allowed by the standard solar model. We discuss the implications for SNO of a low mass, Delta m^2 ~ 6 times 10^{-12} eV^2, vacuum oscillation solution, previously found by Raghavan, and by Krastev and Petcov, but absent in recent analyses that included Super-Kamiokande data. For the SNO experiment, we present refined predictions for the charged-current rate and the ratio of the neutral-current rate to charged-current rate. The predicted charged-current rate can be clearly distinguished from the no-oscillation rate only for the LMA solution. The predicted ratio of the neutral-current rate to charged-current rate is distinguishable from the no-oscillation ratio for the LMA, SMA, LOW, and VAC solutions for active neutrinos.

We have performed a detailed study of the zenith angle dependence of the regeneration factor and distributions of events at SNO and SK for different solutions of the solar neutrino problem. In particular, we discuss oscillatory behaviour and the synchronization effect in the distribution for the LMA solution, the parametric peak for the LOW solution, etc.. Physical interpretation of the effects is given. We suggest a new binning of events which emphasizes distinctive features of zenith angle distributions for the different solutions. We also find the correlations between the integrated day-night asymmetry and the rates of events in different zenith angle bins. Study of these correlations strengthens the identification power of the analysis.

We study properties of the oscillation effects in the matter of the Earth on antineutrino fluxes from supernovae. We show that these effects can provide explanation of the difference in the energy spectra of the events detected by Kamiokande-2 and IMB detectors from SN1987A as well as the absence of high-energy events with $E\gta 40$ MeV. This explanation requires the neutrino oscillation parameters $\Delta m^2$ and $\sin^2 2\theta $ to be in the region of the LMA solution of the solar neutrino problem and the normal mass hierarchy if $|U_{e 3}|^2\gta 10^{-3}$. The hierarchy can be inverted if $|U_{e 3}|^2\ll 10^{-3}$. The solution of the solar neutrino problem based on $\nu_e$-conversion to a pure sterile state is disfavoured by SN1987A data.

Neutrino oscillation scenarios predict correlations, and zones of avoidance, among measurable quantities such as spectral energy distortions, total fluxes, time dependences, and flavor content. The comparison of observed and predicted correlations will enhance the diagnostic power of solar neutrino experiments. A general test of all presently-allowed (two neutrino) oscillation solutions is that future measurements must yield values outside the predicted zones of avoidance. To illustrate the discriminatory power of the simultaneous analysis of multiple observables, we map currently allowed regions of neutrino masses and mixing angles onto planes of quantities measurable with the Sudbury Neutrino Observatory (SNO). We calculate the correlations that are predicted by vacuum and MSW (active and sterile) neutrino oscillation solutions that are globally consistent with all available neutrino data. We derive approximate analytic expressions for the dependence of individual observables and specific correlations upon neutrino oscillations parameters. We also discuss the prospects for identifying the correct oscillation solution using multiple SNO observables.

We study the phenomenological consequences of maximal and near-maximal mixing of the electron neutrino with other ($x$=tau and/or muon) neutrinos. We describe the deviations from maximal mixing in terms of a parameter $\epsilon\equiv1-2\sin^2\theta_{ex}$ and quantify the present experimental status for $|\epsilon|<0.3$. We find that the global analysis of solar neutrino data allows maximal mixing with confidence level better than 99% for $10^{-8}$ eV$^2\lsim\Delta m^2\lsim2\times10^{-7}$ eV$^2$. In the mass ranges $\Delta m^2\gsim 1.5\times10^{-5}$ eV$^2$ and $4\times10^{-10}$ eV$^2\lsim\Delta m^2\lsim2\times10^{-7}$ eV$^2$ the full interval $|\epsilon|<0.3$ is allowed within 4$\sigma$(99.995 % CL). We suggest ways to measure $\epsilon$ in future experiments. The observable that is most sensitive to $\epsilon$ is the rate [NC]/[CC] in combination with the Day-Night asymmetry in the SNO detector. With theoretical and statistical uncertainties, the expected accuracy after 5 years is $\Delta \epsilon\sim 0.07$. We also discuss the effects of maximal and near-maximal $\nu_e$-mixing in atmospheric neutrinos, supernova neutrinos, and neutrinoless double beta decay.

We report on the first space-VLBI observations of the OH34.26+0.15 maser in two main line OH transitions at 1665 and 1667 MHz. The observations involved the space radiotelescope on board the Japanese satellite HALCA and an array of ground radio telescopes. The map of the maser region and images of individual maser spots were produced with an angular resolution of 1 milliarcsec which is several times higher than the angular resolution available on the ground. The maser spots were only partly resolved and a lower limit to the brightness temperature 6x10^{12} K was obtained. The maser seems to be located in the direction of low interstellar scattering, an order of magnitude lower than the scattering of a nearby extragalactic source and pulsar.

We calculate the range of predicted values for 10 quantities that will be measured by the Sudbury Neutrino Observatory (SNO). We use neutrino oscillation solutions (vacuum and MSW; active and sterile neutrinos) that are globally consistent with all available neutrino data and estimate realistic theoretical and experimental uncertainties. The neutral current to charged current double ratio is predicted to be more than 9 sigma from the no-oscillation solution for all of the currently favored neutrino oscillation solutions. The best-fit oscillation solutions predict a CC day-night rate difference between -0.1% and +12.5% and a NC day-night difference < 0.01%. We present also the predicted range for the first and the second moments of the charged current electron recoil energy spectrum, the charged current, the neutral current, and the neutrino-electron scattering rates, the seasonal dependence of the charged current rate, and the double ratio of neutrino-electron scattering rate to charged current rate.

The ratio of the measured to the predicted standard model CC event rates in SNO will be 0.47 if no oscillations occur. The best-fit active oscillation predictions for the CC ratio are: 0.35-39 (MSW) and 0.38-42 (vacuum) (all for a 5 MeV energy threshold), typically about 20% less than the no-oscillation expectation. We calculate the predicted ratios for six active and sterile neutrino oscillation solutions allowed at 99% CL and determine the dependence of the ratios on energy threshold. If the high-energy anomaly observed by SuperKamiokande is due to an enhanced hep flux, MSW active solutions predict that out of a total of 5000 CC events above 5 MeV in SNO between 49 and 54 events will be observed above 13 MeV whereas only 19 events are expected for no-oscillations and a nominal standard hep flux.

We show that the conditions for total neutrino conversion found in [1] are equivalent to the conditions of maximal depth (parametric resonance) and ($\pi/2 + \pi k$) - phase of parametric oscillations. Therefore the effects considered in [1] are a particular case of the parametric resonance in neutrino oscillations. The existence of strong enhancement peaks in transition probability P rather than the condition P=1 is of physical relevance. We comment on possible realizations and implications of the parametric enhancement of neutrino oscillations.

Recent results on solar neutrinos provide hints that the LMA MSW solution could be correct. We perform accurate calculations for potential `smoking-gun' effects of the LMA solution in the SuperKamiokande solar neutrino experiment, including: (1) an almost constant reduction of the standard recoil electron energy spectrum (with a weak, < 10%, relative increase below 6.5 MeV); (2) an integrated difference in day-night rates (2% to 14%); (3) an approximately constant zenith-angle dependence of the nighttime event rate; (4) a new test for the difference in the shape of the equally-normalized day-night energy spectra (~ 1%); and (5) annual variations of the signal due to the regeneration effect (~ 6 times smaller than the integrated day-night effect). We also establish a relation between the integrated day-night asymmetry and the seasonal asymmetry due to LMA regeneration. As a cautionary example, we simulate the effect of an absolute energy calibration error on the shape (distortion) of the recoil energy spectrum. We compare LMA predictions with SuperKamiokande data and discuss the possibilities for distinguishing experimentally between LMA and vacuum oscillations.

We consider main ingredients which determine neutrino transformations in media. Strong transformations relevant for the astrophysics can be due to large depth oscillations, resonance conversion, parametric resonance effect, interplay of oscillations and inelastic collisions. Properties of transitions are discussed using the graphic representation. The applications of the transitions to supernova neutrinos are described. The supernova neutrinos can probe whole neutrino mass spectrum. Their studies will help to identify the pattern of neutrino mass and mixing.

We determine the neutrino parameters for MSW and vacuum oscillations (active and sterile neutrinos) that are allowed by the separate, and collective, imposition of the constraints from total event rates in the chlorine, GALLEX, SAGE, and SuperKamiokande experiments (504 days), the SuperKamiokande electron energy spectrum, and the SuperKamiokande zenith-angle dependence. The small mixing angle MSW solution is acceptable at the 7% C.L. (8% for sterile nu's) and the vacuum solution is acceptable at the 6% C.L. . The best-fit global MSW solution for active neutrinos is: Delta m^2 = 5 x 10^-6 eV^2, sin^2 (2 theta) = 5.5 x 10^{-3} (and for sterile neutrinos: Delta m^2 = 4 x 10^-6 eV^2, sin^2 (2 theta) = 7 x 10^-3). For vacuum oscillations, the best-fit solution is: Delta m^2 = 6.5 x 10^-11 eV^2, sin^2 (2 theta) = 0.75 . An arbitrary combination of undistorted (no oscillations) pp, 7Be, 8B, and CNO neutrino fluxes is inconsistent with the combined data sets at the 3.5 sigma C.L., independent of astrophysical considerations. We use improved calculations of solar model fluxes, neutrino absorption cross sections and energy spectra, and a detailed evaluation of regeneration effects.

We propose a new mechanism of leptogenesis in which the asymmetries in lepton numbers are produced through the CP-violating oscillations of ``sterile'' (electroweak singlet) neutrinos. The asymmetry is communicated from singlet neutrinos to ordinary leptons through their Yukawa couplings. The lepton asymmetry is then reprocessed into baryon asymmetry by electroweak sphalerons. We show that the observed value of baryon asymmetry can be generated in this way, and the masses of ordinary neutrinos induced by the seesaw mechanism are in the astrophysically and cosmologically interesting range. Except for singlet neutrinos, no physics beyond the Standard Model is required.

We consider the Grand Unification (GU) scenario for neutrino masses which is based on the see-saw mechanism with the mass of the heaviest right handed (RH) neutrino at the GU-scale: $M_3 \sim \Lambda_{GU}$, and on the quark-lepton symmetry for fermions from the third generation. The scenario predicts for the light neutrinos: $m_3 \sim (2 - 4) \cdot 10^{-3}$ eV and $m_2 \sim (0.3 - 3) \cdot 10^{-5}$ eV (in the case of a linear mass hierarchy of the RH neutrinos or/and in presence of the Planck scale suppressed non-renormalizable operators). It also predicts large $\nu_e - \nu_{\mu}$ mixing: $~\sin^2 2\theta_{e\mu} \stackrel{_>}{_\sim} 0.2$. In this scenario the solar neutrinos ($\nu_{\odot}$) undergo both the \nue $\to$ \nutau resonance conversion in the Sun and substantial \nue $\to$ \numu vacuum oscillations on the way from the Sun to the Earth. The interplay of both effects enlarges the range of neutrino parameters which solve the $\nu_{\odot}$-problem. In particular, $\nu_e - \nu_{\tau}$ mixing angle can be as small as the corresponding quark mixing: $\sin^2 2\theta_{e\tau} \geq (2~-~5) \cdot 10^{-4}$. The scenario predicts peculiar (oscillatory) distortion of the boron neutrino energy spectrum and seasonal variations of signals. Manifestations of these effects in the Super-Kamiokande and SNO experiments are studied.

Electron polarization induced by magnetic fields can modify the potentials relevant for describing neutrino conversions in media with magnetic fields. The magnitudes of polarization potentials are determined for different conditions. We show that variations of the electron polarization along the neutrino trajectory can induce resonant conversions in the active-sterile neutrino system, but cannot lead to level crossing in the active-active neutrino system. For neutrino flavour conversions the polarisation leads only to a shift of the standard MSW resonance. For polarizations $\lambda \lsim 0.04$ the direct modifications of the potential (density) due to the magnetic field pressure are smaller than the modifications due to the polarization effect. We estimate that indeed the typical magnitude of the polarization in the sun or in a supernova are not expected to exceed $10^{-2}$. However even such a small polarization may lead to interesting consequences for supernova physics and for properties of neutrino signals from collapsing stars.

We study the spin-flavor transitions of neutrinos in the magnetic fields of a collapsing star. For the neutrino mass squared difference (10^{-10} - 10) eV^2 the transitions take place in almost isotopically neutral region of the star, where the effective matter density is suppressed up to 3 - 4 orders of magnitude. This suppression is shown to increase the sensitivity of the neutrino bursts studies to the magnetic moment of neutrino by 1.5 - 2 orders of magnitude, and for realistic magnetic field the observable effects may exist for (2 - 3)10^{-14} Bohr magneton. In the isotopically neutral region the jumps of the effective potential exist which influence the probabilities of transitions. The experimental signatures of the spin-flavor transitions are discussed. In particular, in the case of direct mass hierarchy, the spin-flip effects result in a variety of modifications of the electron antineutrino spectrum. Taking this into account, we estimated the upper bounds on the magnetic moment from the SN1987A data. In the isotopically neutral region the effects of possible twist of the magnetic field on the way of neutrinos can be important, inducing distortion of the neutrino energy spectra and further increasing the sensitivity. However, if the total rotation angle is restricted the absolute change of probabilities is small.

There are large uncertainties in the predictions of the boron neutrino flux from the Sun which cannot be considered as being of purely statistical origin. We treat the magnitude of this flux, $\Phi_B$, as a parameter to be found from experiment. The properties of the MSW solution to the solar neutrino problem for different values of $\Phi_B$ are studied. Present data give the bounds: $0.38 < \Phi_B/ \Phi_B^0 < 3.1$ (2$\sigma$), where $\Phi_B^0 \equiv 5.7 \cdot 10^6$ cm$^{-2}$s$^{-1}$ is the flux in the reference SSM. The variations of the flux in this interval enlarge the allowed region of mixing angles: $\sin^2 2\theta = 8\cdot 10^{-4} \div 2\cdot 10^{-2}$ (small mixing solutions) and $\sin^2 2\theta = 0.2 \div 0.85$ (large mixing solution). If the value of the original boron neutrino flux is about that measured by Kamiokande, a consistent description of the data is achieved for $\sin^2 2\theta \sim (0.8 \div 2)\cdot 10^{-3}$ (``very small mixing solution"). The solution is characterized by a strong suppression of the beryllium neutrino line, a weak distortion of the high energy part of the boron neutrino spectrum and a value of the double ratio $(CC/NC)^{exp}/(CC/NC)^{SSM}$ at $E > 5$ MeV close to 1. We comment on the possibility to measure the neutrino parameters and the original boron neutrino flux in future experiments.

Transitions in a system of neutrinos with vacuum mixing and magnetic moments, propagating in matter and transverse magnetic field, are considered. It is shown that in the realistic case of magnetic field direction varying along the neutrino path qualitatively new phenomena become possible: permutation of neutrino conversion resonances, appearance of resonances in the neutrino-antineutrino ($\nu_{lL}\leftrightarrow\bar{\nu}_{lR}$) transition channels, neutrino-antineutrino resonant conversion, large amplitude $\nu_{lL}\leftrightarrow\bar{\nu}_{lR}$ oscillations, merging of different resonances (triple resonances). Possible phenomenological implications of these effects are briefly discussed.