New Langevin and Gradient Thermostats for Rigid Body Dynamics

We introduce two new thermostats, one of Langevin type and one of gradient (Brownian) type, for rigid body dynamics. We formulate rotation using the quaternion representation of angular coordinates; both thermostats preserve the unit length of quaternions. The Langevin thermostat also ensures that the conjugate angular momenta stay within the tangent space of the quaternion coordinates, as required by the Hamiltonian dynamics of rigid bodies. We have constructed three geometric numerical integrators for the Langevin thermostat and one for the gradient thermostat. The numerical integrators reflect key properties of the thermostats themselves. Namely, they all preserve the unit length of quaternions, automatically, without the need of a projection onto the unit sphere. The Langevin integrators also ensure that the angular momenta remain within the tangent space of the quaternion coordinates. The Langevin integrators are quasi-symplectic and of weak order two. The numerical method for the gradient thermostat is of weak order one. Its construction exploits ideas of Lie-group type integrators for differential equations on manifolds. We numerically compare the discretization errors of the Langevin integrators, as well as the efficiency of the gradient integrator compared to the Langevin ones when used in the simulation of rigid TIP4P water model with smoothly truncated electrostatic interactions. We observe that the gradient integrator is computationally less efficient than the Langevin integrators. We also compare the relative accuracy of the Langevin integrators in evaluating various static quantities and give recommendations as to the choice of an appropriate integrator.Comment: 16 pages, 4 figure

Irreducible Representations of Diperiodic Groups

The irreducible representations of all of the 80 diperiodic groups, being the symmetries of the systems translationally periodical in two directions, are calculated. To this end, each of these groups is factorized as the product of a generalized translational group and an axial point group. The results are presented in the form of the tables, containing the matrices of the irreducible representations of the generators of the groups. General properties and some physical applications (degeneracy and topology of the energy bands, selection rules, etc.) are discussed.Comment: 30 pages, 5 figures, 28 tables, 18 refs, LaTex2.0

Implications of recent solar neutrino observations: an analysis of charged current data

We have analysed the recent results from the observation of charged current \nu_e d \to e^- p p events from solar neutrinos by the Sudbury Neutrino Observatory SNO assuming neutrino oscillations with three active flavours. The data seem to prefer a low mass-squared difference and large mixing angle solution (the so-called LOW solution) in (12) parameter space. However, when combined with the Gallium charged current interaction data from Gallex and GNO, distinct (1\sigma) allowed regions corresponding to the large mixing angle (LMA) and small mixing angle (SMA) appear while the LOW solution is disfavoured upto 3\sigma standard deviation. The physical electron neutrino survival probability corresponding to these best fit solutions are then determined and analysed for their energy dependence.Comment: 16 pages Latex file, with 5 epsf figures; one reference adde

Synthesis of multi-loop automatic control systems by the nonlinear programming method

The article deals with the problem of calculation of the multi-loop control systems optimal tuning parameters by numerical methods and nonlinear programming methods. For this purpose, in the paper the Optimization Toolbox of Matlab is used

On Hilberg's Law and Its Links with Guiraud's Law

Hilberg (1990) supposed that finite-order excess entropy of a random human text is proportional to the square root of the text length. Assuming that Hilberg's hypothesis is true, we derive Guiraud's law, which states that the number of word types in a text is greater than proportional to the square root of the text length. Our derivation is based on some mathematical conjecture in coding theory and on several experiments suggesting that words can be defined approximately as the nonterminals of the shortest context-free grammar for the text. Such operational definition of words can be applied even to texts deprived of spaces, which do not allow for Mandelbrot's ``intermittent silence'' explanation of Zipf's and Guiraud's laws. In contrast to Mandelbrot's, our model assumes some probabilistic long-memory effects in human narration and might be capable of explaining Menzerath's law.Comment: To appear in Journal of Quantitative Linguistic

A New 76Ge Double Beta Decay Experiment at LNGS

This Letter of Intent has been submitted to the Scientific Committee of the INFN Laboratori Nazionali del Gran Sasso (LNGS) in March 2004. It describes a novel facility at the LNGS to study the double beta decay of 76Ge using an (optionally active) cryogenic fluid shield. The setup will allow to scrutinize with high significance on a short time scale the current evidence for neutrinoless double beta decay of 76Ge using the existing 76Ge diodes from the previous Heidelberg-Moscow and IGEX experiments. An increase in the lifetime limit can be achieved by adding more enriched detectors, remaining thereby background-free up to a few 100 kg-years of exposure.Comment: 67 pages, 19 eps figures, 17 tables, gzipped tar fil

On the massless "just-so" solution to the solar neutrino problem

We study the effect of the non-resonant, vacuum oscillation-like neutrino flavor conversion induced by non-standard flavor changing and non-universal flavor diagonal neutrino interactions with electrons in the sun. We have found an acceptable fit for the combined analysis for the solar experiments total rates, the Super-Kamiokande (SK) energy spectrum and zenith angle dependence. Phenomenological constraints on non-standard flavor changing and non-universal flavor diagonal neutrino interactions are considered.Comment: 4 pages, Latex, uses eps

Continuous, Semi-discrete, and Fully Discretized Navier-Stokes Equations

The Navier--Stokes equations are commonly used to model and to simulate flow phenomena. We introduce the basic equations and discuss the standard methods for the spatial and temporal discretization. We analyse the semi-discrete equations -- a semi-explicit nonlinear DAE -- in terms of the strangeness index and quantify the numerical difficulties in the fully discrete schemes, that are induced by the strangeness of the system. By analyzing the Kronecker index of the difference-algebraic equations, that represent commonly and successfully used time stepping schemes for the Navier--Stokes equations, we show that those time-integration schemes factually remove the strangeness. The theoretical considerations are backed and illustrated by numerical examples.Comment: 28 pages, 2 figure, code available under DOI: 10.5281/zenodo.998909, https://doi.org/10.5281/zenodo.99890

The luminosity constraint on solar neutrino fluxes

A specific linear combination of the total solar neutrino fluxes must equal the measured solar photon luminosity if nuclear fusion reactions among light elements are responsible for solar energy generation. This luminosity constraint, previously used in a limited form in testing the no neutrino oscillation hypothesis, is derived in a generality that includes all of the relevant solar neutrino fluxes and which is suitable for analyzing the results of many different solar neutrino experiments. With or without allowing for neutrino oscillations, the generalized luminosity constraint can be used in future analyses of solar neutrino data. Accurate numerical values for the linear coefficients are provided.Comment: related material at http://www.sns.ias.edu/~jn