abstract: In this review, we summarize the observed statistical properties of visual binaries in various young stellar clusters. The binary frequency, orbital period distribution and mass ratio distribution are considered for populations of both low-mass and high-mass stars. These properties are then compared to numerical models of the dynamical evolution of stellar clusters and to the predictions of binary formation mechanisms.
  .
  : visual binaries – statistical properties – young stellar clusters – dynamical evolution – binary formation mechanisms
author: Gaspard Duchêne & Jérôme Bouvier; Jochen Eislöffel; Theodore Simon
title: Statistical properties of visual binaries as tracers of the formation and early evolution of young stellar clusters

# Introduction

The current paradigm for star formation is that most stars form in stellar clusters rather than in Taurus-like star-forming regions. Since binary systems are more frequent than single stars, the study of multiple systems in clusters is expected to provide us with some insight on the stellar formation process itself. For instance, it has been suggested that the gravitational collapse of molecular clouds leads to a binary fraction close to 100 %, which would decrease over time to about 60 %, the value observed on the Main Sequence (Ghez *et al.* 1993). The physical processes leading to such an evolution and the timescale over which they act still need to be asserted. Alternatively, the impact of local physical conditions in the molecular cloud, such as the magnetic field strength or gas temperature, on the output of the star formation process has still to be ascertain and might lead to different binary frequencies in various environments (Durisen & Sterzik 1994).

Observing stellar clusters with a wide range of physical parameters allows to perform a direct comparison between binary properties and model predictions. Clusters with ages ranging from one to several hundred million years are known within 500 pc from the Sun with some hosting many massive stars while other containing only low-mass stars. The stellar density can vary over more than an order of magnitude from cluster to cluster, as well as the physical conditions of the parent molecular cloud. In this contribution, we review the observed properties of visual binaries in clusters and compare them with cluster properties and with model predictions regarding the formation and evolution of stellar clusters.

# Summary of observational results

## Visual binary frequency

The advent of several high-angular imaging techniques makes the binary frequency a relatively straightforward property to derive in stellar populations. For instance, $HST$ imaging (Prosser *et al.* 1994), speckle interferometry (Leinert *et al.* 1993, Patience *et al.* 1998) and adaptive optics (Bouvier *et al.* 1997) have been used, mostly in the near-infrared in order to facilitate the detection of relatively faint and cold low-mass companions. The low-mass populations of extremely young clusters (the Trapezium cluster, IC 348) as well as of more evolved ones ($\alpha$ Per, the Pleiades, the Hyades, Præsepe) have been studied, which span an age range from $\sim$`<!-- -->`{=html}1 to several 100 Myr. The range of semi-major axes probed by these surveys usually extends from a few 10 to $\sim$`<!-- -->`{=html}2000 AU. The Hyades cluster, which is the closest to the Sun, is the only one for which visual binaries could be resolved down to about 5 AU.

Most surveys miss the faintest companions so that a statistical completeness correction has to be applied in order to estimate the true binary frequency from the observed one. This correction can be large, especially at small separations for adaptive optics surveys due to the separation-dependent detection limit inherent to this technique. The correction factor also depends upon assumptions regarding the distribution of mass ratios and/or orbital periods of the binary population, which obviously increases the statistical uncertainties on the derived binary frequency.

With these limitations, all studied clusters are found to exhibit a visual binary frequency consistent with that of field G-dwarfs (Duquennoy & Mayor 1991), to within 10 % or less (e.g. Prosser *et al.* 1994, Bouvier *et al.* 1997, Duchêne *et al.* 1999). This result does not seem to depend upon the age of the cluster between $\sim$`<!-- -->`{=html}1 Myr and 600 Myr. One possible exception might be the Hyades cluster which appears to host more tight (5–50 AU) binaries than the field population does, though statistical uncertainties could account for this slight excess (Patience *et al.* 1998). It hence appears that the fraction of binaries in very young clusters, at an age of about a million years, is already similar to that measured for field dwarfs and that it does not significantly evolve as the clusters dynamically relax on the ZAMS.

## Orbital period distribution

While the overall binary frequency integrated over a range of semi-major axes already provides clues to the properties of binary populations in clusters, one would ideally wish to compare the distribution of their orbital periods. One is thus led to convert the observed angular separations into orbital periods assuming an average system mass and applying a statistical correction factor to take into account projection effects.

Over a limited range of orbital periods, the distribution of semi-major axes in cluster binaries again appears to be consistent with that of field G-dwarf binaries (Duquennoy & Mayor, 1991). The comparison between cluster and field binaries can be further investigated by considering not only visual binaries but spectroscopic binaries as well. This has been done for ZAMS clusters, where the low-mass stars have already been searched for short-period companions. Within statistical uncertainties, the results reinforce the similarity between cluster and field binaries, as shown for the Pleiades in Figure 1.

In order to improve the statistics, we combined the results obtained for several ZAMS clusters. This yields the orbital period distribution shown in Figure 2 for visual binaries. While the distribution of period in cluster binaries agrees with that of field binaries over the orbital range $10^{4.4}$–$10^{6.4}\,$days, a paucity of wide (sep $\geq$ 300 AU) cluster binaries appears relative to the field at the 2.5 $\sigma$ level. Such a deficit cannot be related to a detection bias, since the distant companions of wide binaries are the most easily detected. The uncompleteness correction which was applied here affects only the first orbital bin of the distribution.

This difference between field and cluster binaries is reinforced by the evidence that the orbital period distribution peaks at a smaller value in clusters than in the field (Patience *et al.* 2000). The mode of the distribution is 5 AU for cluster binaries and 40 AU for field binaries, the former being estimated under the assumption that all clusters, including the Hyades, can be combined into a single sample. Furthermore, Scally *et al.* (1999) have shown that there are no very wide binaries (sep $\geq$ 2000 AU) in the Trapezium cluster.

Hence, even though clusters at various evolutionary stages appear to host a roughly similar fraction of low-mass binary systems, which itself is similar to the field binary fraction, the distribution of orbital periods appears to exhibit small but significant differences between clusters and field binary populations.

## Mass ratio distribution

The mass of the primaries and secondaries of binary systems on the ZAMS can be reliably estimated from their near-infrared flux using main sequence mass-luminosity relationships. In order to increase the statistical significance of the derived mass ratio distribution, the combined results for three clusters observed with similar set-ups are shown in Figure 3. The observed distribution appears to be essentially flat over the range $q=M_b/M_a=0.1$–1. This agrees with the $q$-distribution derived by Duquennoy & Mayor (1991) for low-mass field binaries if the limited sensitivity of our observations, which prevents us from detecting tight systems with low mass ratios, is taken into account.

Contrary to main sequence binaries, young low-mass systems frequently exhibit significant continuum excesses at near-infrared wavelengths and this may affect the derivation of mass ratios from near-IR flux ratios. However, it is usually assumed on a statistical basis that both stars in a young binary system exhibit similar excesses, and that the overall flux ratio distribution can be safely transformed into a mass ratio distribution. Baraffe *et al.* (1998)'s models were used to get the mass-luminosity relationship relevant for such young stars. Both in the Taurus star-forming region (Leinert *et al.* 1993) and in the Trapezium cluster (Prosser *et al.* 1994), the observed $q$-distribution is consistent with being flat down to $q=0.1$, and therefore not different from the ZAMS and field estimates.

Current results thus suggest that the mass ratio distribution is similar for all low-mass binary populations studied so far, independent of their evolutionary status and environmental conditions (field or clusters). No peak is found in the distribution towards $q\approx1$ nor towards low $q$. Furthermore, the ZAMS clusters mass ratio distribution does not seem to vary with separation for visual binaries (Patience *et al.* 2000).

## Massive binaries

Studies conducted so far have mostly focussed on the low-mass binary population of clusters. More recently, however, the investigation of massive binaries has started, with the additional difficulty of a huge contrast between the bright primary and the low-mass companion which is a major limitation of current massive binary surveys.

The recent *HST* study of thirteen OB stars in the Trapezium cluster by Preibisch *et al.* (1999) already revealed eight visual companions with separations ranging from 15 to 450 AU. The binary frequency of high-mass stars in this cluster thus appears significantly higher than that of their low-mass counterparts (Prosser *et al.* 1994), and the difference is further enhanced when completeness corrections are taken into account. Similarly, an adaptive optics survey of about fifty OB stars in the distant cluster NGC 6611 revealed eight companions in the separation range 200–2000 AU (Duchêne 2000), indicative of a high binary fraction compared to G-type field dwarfs if the orbital period distribution is the same for high-mass and low-mass binaries. A study of young, mostly isolated, Herbig AeBe stars conducted by Bouvier & Corporon (2000) yields an uncorrected binary frequency which is significantly higher than that of lower mass Main Sequence stars for both Be and Ae PMS stars.

The mass-ratio distribution of massive binaries can be estimated from the flux ratios measured by imaging surveys. The $q$-distribution derived for massive binaries with B0–B7 primaries is shown in Figure 4. The sample includes members of the Trapezium cluster and NGC 6611 as well as isolated Herbig Be stars, with a semi-major axis ranging from 100 to 2000 AU. The mass ratio distribution is strongly peaked towards small values: only about 25 % of the binaries have $q>0.5$. Patience *et al.* (2000) obtained a similar trend for A and F stars in various ZAMS clusters. This distribution is dramatically different from that of low-mass binaries, indicating that the primary's mass is a relevant parameter when investigating the formation and evolution of binary systems.

# Interpretation

The observed properties of young binary systems hold clues to their formation mechanism and their early evolution. In the following, we consider both classes of models, binary formation and dynamical evolution of clusters, and compare their predictions with observations.

## Dynamical evolution of stellar clusters

The fact that all clusters studied so far exhibit a similar fraction of visual binaries regardless of their age rules out a long-term decrease of the binary frequency past the first million years after gravitational collapse. In addition, since similarly-aged loose associations and clusters display different binary frequencies, it is tempting to conclude that environmental effects play a key role in fixing the binary content of a star forming region.

$N$-body simulations of the dynamical evolution of stellar clusters were ran by Kroupa *et al.* (1999) and compared with observations of the Trapezium cluster. They found that close encounters between cluster members could efficiently disrupt binaries in less than a million years while the number of binaries would remain constant afterwards. This results from the fact that the number of encounters falls off dramatically after one relaxation time. If the initial low-mass binary fraction is of the order of 100 % at the end of the protostellar collapse, such a scenario would naturally explain the high binary frequency observed in T Tauri associations: in these regions, the stellar density is several orders of magnitudes too low to generate numerous enough close-by interactions so as to impact on the initial binary content (Clarke & Pringle 1991). In contrast, frequent gravitational encounters in dense clusters would rapidly erode the initially high binary fraction and lead to a binary frequency similar to that of the field after a million years of evolution or so.

The observed properties of binary populations described in Section 2 further support the hypothesis that dynamical evolution of clusters may largely drive the binary properties. In dynamical models, wide binaries are preferentially disrupted, because their binding energy is small compared to the average kinetic energy of the cluster (see the review by P. Kroupa in this volume). This prediction agrees with the observed distribution of orbital periods for binaries in ZAMS clusters which exhibits a paucity of wide systems, as shown in Figure 5. The higher binary fraction found for massive stars is also in qualitative agreement with this model since the potential well of massive primaries is deeper than that of solar-mass stars, making it more difficult to disrupt more massive binaries. The "softening boundary" which delineates hard and soft systems, the latter being prone to disruption, is thus moved towards wider separations.

A number of issues remain however to be addressed before accepting this paradigm. For instance, dynamical models predict that low-$q$ systems are preferentially disrupted, introducing a non-negligible modification in the shape of this distribution. This constrasts with the derivation of similar mass ratio distributions for low-mass binaries in clusters and associations. Quantitatively, Kroupa (1995) estimated that 70 % of the initial binaries must have a mass ratio $q<0.4$ in order to account for the observed distribution of field binaries. This is unlike the $q$-distribution observed for T Tauri associations. Another puzzling result is the fact that all clusters seem to have very similar binary frequency and orbital period distributions. Hence, binary properties appear to be weakly dependent on small differences in initial conditions. Dyamical models have focussed on reproducing the binary population of the Trapezium cluster and it remains to be seen whether the same results would be obtained for slightly less dense clusters such as IC 348.

According to these models, wide binaries found in the solar neighbourhood cannot have formed in clusters, since they would have been rapidly disrupted. One thus has to assume that they were born in loose associations that have subsequently dissipated in the field. However, the field binary population does not appear to be a simple combination of cluster and association members: a weighted mixture of these 2 populations that reproduces the observed fraction of tight and wide binaries also yields a large excess of systems in the 5–50 AU separation range compared to field binaries. The issue of the efficiency of star formation in clusters as opposed to associations is still open.

## Binary formation models

If the dynamical evolution of clusters is responsible for the observed binary properties, then the formation process itself must yield a large initial binary fraction. Furthermore, the masses of the components must be taken at random from the initial mass function, as suggested by Kroupa *et al.* (1995). We note that the observations of high-mass binaries are qualitatively consistent with random pairing, with a limited number of $q>0.5$ pairs, though this result is still preliminary due to the small sample size. We will now discuss whether binary formation mechanisms do predict such properties for initial binaries, prior to cluster's evolution.

The formation of a stellar cluster probably results from the fragmentation of the parent molecular cloud into a large number of individual cores. In this environment, wide binary systems can form either by capture or through the fragmentation of a small collapsing core. Models considering capture in small-$N$ ($N<10$) subclusters predict that the two most massive members are usually associated in a binary system, while the other stars are ejected (McDonald & Clarke 1993). This leads to a mass ratio distribution which either peaks towards unity or is roughly flat, depending on the value of $N$. However, taking into account the dissipative effect of protostellar disks during close encounters leads to a secondary mass distribution that mimics the IMF (McDonald & Clarke 1995). Furthermore, the model predicts a higher binary frequency for more massive stars. Similar simulations made by Sterzik & Durisen (1998) additionally predict that the mass ratio distribution is more heavily weighted towards low-$q$ values in higher mass systems. All these predictions are consistent with the observations, indicating that stars in clusters may form in compact subclusters. However, such a scenario does not end up with a 100 % initial binary frequency, since some single stars are ejected from each subcluster.

Capture during the relaxation of small-$N$ subgroups appears to be the model that best fits current observations, and especially the random pairing of secondary masses. The alternative model of scale-free fragmentation (Clarke 1996) predicts that binary properties do not depend on the primary mass, which conflicts with the observations. Binary properties could also result, at least partly, from the accretion of residual material onto two gravitationally bound seeds (Bate & Bonnell 1997, Bate 2000). However, these models predict a larger number of $q\sim1$ binaries for tighter systems as well as for more massive binaries, which is not verified in the observational trends. Finally, binaries formed through disk fragmentation should preferentially have small separations and low mass ratios, a prediction which contrasts with the apparent lack of a dependence of the mass ratio distribution upon separation.

# Summary

The advent of high-angular resolution imaging techniques have offered the opportunity to systematically survey binaries in young clusters, thus complementing spectroscopic studies. We reviewed some statistical properties of visual binaries such as their frequency, orbital period distribution and mass ratio distribution; both low-mass and high-mass stars were considered. Observational results can be summarized as follows: the low-mass binaries in cluster have similar properties as Main Sequence (field) binaries, except for a significant deficit of systems with a separation larger than 300 AU. High-mass stars host many more companions than do solar-type stars, and most massive binaries have a small mass ratio. These observational constraints are consistent with the results of numerical simulations describing the early dynamical evolution of stellar clusters provided that the initial binary frequency is of the order of 100 % for low-mass stars. Disk-assisted capture in small subclusters is the formation model that best approaches the initial conditions required to account for the observations. In both classes of models, however, several issues still need to be investigated and additional simulations over a wider range of physical conditions have to be run before these scenarios can be considered as offering an acceptable new paradigm for the star and binary formation process. In parallel, future observations of young embedded clusters will help to derive the properties of protobinaries at an even earlier stage of evolution than discussed in this review, which will undoubtly lead to new clues to their formation.


author: Murray Gell-Mann[^1][^2]; Pierre Ramond$^\dagger$; Richard Slansky
title: COMPLEX SPINORS AND UNIFIED THEORIES

In *Supergravity*, P. van Nieuwenhuizen & D.Z . Freedman (eds.) *North Holland Publishing Company* 1979

We were told by Frank Yang in his welcoming speech that supergravity is a phenomenon of theoretical physics. Why, at this time, is it not more than that? Self-coupled extended supergravity, especially for $N = 8$, seems very close to the overall unified theory for which all of us have yearned since the time of Einstein. There are no quanta of spin $>2$; there is just one graviton of spin $2$; there are $N$ gravitini of spin $3/2$, just right for eating the $N$ Goldstone fermions of spin $1/2$ that are needed if $N$-fold supersymmetry is to be violated spontaneously; there are $N(N -1)/ 2$ spin $1$ bosons, perfectly suited to be the gauge bosons for $SO_N$ in the theory with self-coupling. There are $N(N -1) (N-2)/ 6$ spin $1/2$ Majorana particles, and with the simplest assignments of charge and colour they include isotopic doublets of quarks and leptons. The theory is highly non-singular in perturbation theory, and the threatened divergence at the level of three loops has not even been demonstrated. The apparently arbitrary cancellation of huge contributions of opposite sign to the cosmological constant (from self-coupling on the one hand and from spontaneous violation of supersymmetry on the other) has been phrased in such an elegant way that it may be acceptable. (Of course, if we follow Hawking et al., we may not even need to cancel out the cosmological constant!)

What is wrong then? Of course the spontaneous violation of $SO_N$ and of supersymmetry is not known to happen in the supergravity theory. But what seems much worse, the spectrum of elementary particles includes too few spin $1$ and spin $1/2$ objects to agree with the list that we would like to see on the basis of our experience at energies $\lesssim~50$ $GeV$. Of course, looking up at the Planck mass of $\sim 2 \times 10^{19}$ $GeV$, we are in a position of greater inferiority than an ant staring up at a skyscraper (facing a factor of only $10^6$ or so) and it may not be reasonable to expect that what looks elementary to us should be elementary on a grand scale.

Nevertheless, we make the comparison and we find that $SO_8\supset  SU_3^C \times U_1\times U_1$ but $SO_8 \not\supset SU_3^c\times  SU_2\times U_1$, so that there is no room for the $X^\pm$ intermediate bosons of the charged current weak interaction. Among the spin $1/2$ particles, we have room for at most two flavours of lepton (say $e$ and $\nu_e$) and four flavours of quark (say $u$, $d$, $c$ and $s$). Even these numbers may be reduced if we try to locate the Goldstone fermions among the elementary spin $1/2$ particles or if we use the generator of the second $U_1$ in a restrictive way.

We would then be forced to regard all or most of the known quarks and leptons as non-elementary, as well as at least two of the intermediate bosons of the weak interaction. The broken Yang-Mills theory of the weak interaction would be only an effective gauge theory, not a fundamental one. All this may prove to be the case, and we will then have to understand the rather complicated relation existing between the elementary particles of the theory and the elementary particles as we perceive them today.

Various investigators have looked into superconformal supergravity, in which one tries to use the full $SU_N$ as a gauge group; such a theory is plagued with particles appearing as multiple poles in propagators, involving difficulties with negative probabilities or lack of causality. Ignoring these serious difficulties, we may ask about the algebraic description of the spin $1/2$ fermions in such a theory. Apparently they are again connected with third-rank tensor representations, forming part of $({\bf N} + {\bf\overline N})_A$ of $SU_N$ (where $A$ means totally antisymmetrized) instead of being assigned to $({\bf N})^3_A$ of $SO_N$.

This tendency to assign the spin $1/2$ fermions to a tensor representation, probably a third rank tensor, of $SO_N$ or $SU_N$, exists even in theories having nothing to do with supergravity. We may, for example, consider a composite model of quarks and leptons, in which they are made up of $N$ kinds of fermionic sub-units. We may think of such a scheme in algebraic terms as assigning these sub-units to the representation ${\bf N}$ of $SO_N$ or $SU_N$ and the quarks and leptons to tensors that are part of $({\bf N})^3$ of $SO_N$ or part of $({\bf N }+ {\bf\overline N})^3$ of $SU_N$, provided each known particle is made up of three sub-units. (Of course one might use five or a higher odd number and obtain fifth rank tensors and so forth, but three is much simpler.)

Now what indications come from the attempts to construct a unified Yang-Mills theory? Do they also point to such a third-rank tensor for the spin $1/2$ fermions?

We turn, then, to the program of formulating a broken Yang-Mills theory of strong and weak interactions, with an effective energy of unification between $10^{14}$ $GeV$ and the Planck mass. This program is only slightly less immodest in conception than the overall unification program of self-coupled supergravity. For the sake of expressing all the Yang-Mills coupling constants in terms of a single one, a simple group $G$ is employed. (Actually one could use $G \times G, G  \times G \times G$, etc. with discrete symmetries connecting the factors, but we shall treat here only the case of a single $G$.)

The smallest $G$ that has been used is $SU_5$; the known left-handed spin $1/2$ fermions are then assigned to three families, each belonging to the reducible representation ${\bf \bar 5} + {\bf 10}$ of $SU_5$, where ${\bf \bar 5}$ contains $\bar d$, $e^-$, and $\nu_e$ for the lowest family, while ${\bf 10}$ consists of $d,~ u,~\bar u$ and $e^+$, The combination ${\bf \bar 5} + {\bf 10}$ is anomaly-free. The violation of symmetry takes place in two stages. First the symmetry $SU_5$ is broken down to $SU^c_3\times SU_2\times  U_1$ by means of a non-zero vacuum expected value of an operator transforming like the adjoint representation ${\bf 24}$, with no direct effect on the fermion masses, and then $SU_2\times U_1$ is broken down to $U_1^{e-m}$ by means of operators transforming like ${\bf 5} + {\bf \bar 5}$, with perhaps an admixture of ${\bf 45}$ giving the masses of the quarks and leptons. The detailed work is done using explicit spinless Higgs boson fields, with various constants for mass, for self-coupling and for coupling to fermions, constants that must be delicately adjusted to make the masses of the fermions and of the intermediate boson for weak interactions tiny with respect to the unification mass. A quantity of roughly similar magnitude, the renormalization-group-invariant mass $\Lambda$ of $QCD$, is tiny with respect to the unification mass for a totally different reason, namely the smallness of the unified coupling constant, which is $\sim  10^{-2}$ near the unification mass, corresponding to the fine structure constant at low energies, and is proportional to the reciprocal of the logarithm of $(10^{14}~ GeV)/\Lambda$. Despite some successes, which we mention below, the $SU_5$ scheme seems to us a temporary expedient rather than a final theory, because of the arbitrariness associated with the Higgs bosons and also because the particles and antiparticles among the left-handed spin $1/2$ fermions have no relation to each other (i.e., there is no $C$ or $P$ operator for the theory).

The $SU_5$ scheme has at least two successes: a roughly correct prediction of the weak angle $\theta_{\rm w}$ and the prediction that after allowing for renormalization $m_b =m_\tau$, which works quite well. The violation of $SU_2\times U_1$ by ${\bf 5} + {\bf \bar 5}$ of $SU_5$ would also give $m_s\sim m_\mu$ and $m_d\sim m_e$ after renormalization. The second of these relations does not work but might be subject to large corrections because the quantities are so small; the first might work if the usual estimates of $m_s$ are in error - otherwise some admixture of a ${\bf 45}$ of $SU_5$ has been suggested, along with ${\bf 5} + {\bf \bar 5}$, but affecting mainly the two lower families.

We have studied various complex spinor schemes that reduce to the $SU_5$ system after some symmetry violation. work on such schemes has also been done by Georgi et al. at Harvard, Susskind and collaborators at Stanford, Wilczek and Zee, Gürsey et al. in the case of $E_6$, and no doubt by many others. Early investigations of complex spinor assignments were carried out by Fritzsch and Minkowski.

First, let us restrict our attention to a single family, say the third one, assuming that the $t$ quark exists and that $\nu_{\tau L}$ is nearly massless. We note that the reason that ${\bf \bar 5} + {\bf 10}$ of $SU_5$ is anomaly-free is that the complex spinor representation ${\bf 16}$ of $SO_{10}$ breaks-up into ${\bf 1} +{\bf  \bar 5} + {\bf 10}$ of $SU_5$, where the singlet can give rise to no $SU_5$ anomaly, and all representations of $SO_n$ $( n \ne 6)$ are anomaly-free.

The 16-dimensional spinor possesses a $C$ symmetry to start with, connecting $\tau^-_L \leftrightarrows \tau^+_L$, $b_L\leftrightarrows (\bar b)_L$, etc., and the 16th particle is just the missing $(\bar\nu_\tau)_L$. Symmetry violations giving fermion masses must correspond to representations contained in the symmetrized square of the fermion representation. We note that in $S0_{10}$ we have $({\bf 16})^2_s ={\bf 10} + {\bf 126}$ and that with respect to $SU_5$ we have ${\bf 10}\rightarrow {\bf 5}+{\bf \bar 5}$ and ${\bf 126}\rightarrow {\bf 1}+{\bf 45}+{\bf 10}+{\bf \overline{15}}+{\bf \bar 5}+{\bf 50}$. An operator transforming like the $SU_5$ singlet piece of ${\bf 126}$ would break the $SO_{10}$ symmetry down to $SU_5$ and would give a Majorana mass term of the form $(\bar\nu_\tau)^2_L+\nu^2_{\tau R}$ to the unobserved neutrino, one that had better be very large if the scheme is to work.

Meanwhile, the ${\bf 10}$ of $SO_{10}$ would give rise to equal Dirac masses for $b$ and $t$ (apart from renormalization) and also to equal Dirac masses for $t$ and $\nu_\tau$. The Dirac mass for the neutrino leads directly to a small effective mass $m{\nu_{\tau L}}
\sim  m_{\rm Dirac}^2/m_{\rm Majorana}$. If $m(\nu_{\tau L}) \approx 1 ~eV$, then neutrinos account for a modest fraction of the missing matter in the universe and give a moderate contribution to the gravitational closure of galaxies and clusters of galaxies. Putting $m_{\rm Dirac}\approx  m_t \approx
~ 30~ GeV$ at a guess, the corresponding value of $m_{\rm Majorana}$ would be $\approx 10^{12}~ Gev$. If $m_{\rm Majorana}$ is very much smaller than that, the cosmological effects become too large; If $m_{\rm Majorana}$ is much larger, that is harmless, but the cosmological effects become negligible.

We can examine $SO_{10}$ in a different way by using the decomposition $SO_{10}\supset SO_6\times SO_4$, where algebraically $SO_6$ is equivalent to $SU_4$ and $SO_4$ to $SU_2\times   SU_2$. We have, then, effectively $SO_{10}\supset SU_2\times SU_2\times SU_4$, where the first $SU_2$ is that of the weak interactions, the second one the corresponding $SU_2$ for left-handed antiparticles or right-handed particles and $SU_4$ is the generalization of $SU^c_3$ introduced in a different connection by Pati and Salam, in which leptons appear as having a fourth colour. The representations of $SO_{10}$ then decompose as follows: $${\bf 16}\rightarrow ({\bf 2},{\bf 1},{\bf 4})+({\bf 1},{\bf 2},{\bf \bar 4}), \qquad {\bf \overline{16}}\rightarrow ({\bf 2},{\bf 1},{\bf \bar 4})+({\bf 1},{\bf 2},{\bf 4}),\qquad  {\bf 10}\rightarrow ({\bf 2},{\bf 2},{\bf 1})+({\bf 1},{\bf 1},{\bf 6})$$ $${\bf 126}\rightarrow ({\bf 2},{\bf 2},{\bf 15})+({\bf 1},{\bf 1},{\bf 6})+({\bf 3},{\bf 1},{\bf 10})+({\bf 1},{\bf 3},{\bf \overline{10}}).$$ The representations ${\bf 1}$, ${\bf 4}$, ${\bf \bar 4}$, ${\bf 15}$, ${\bf 10}$ and ${\bf \overline{10}}$ of $SU_4$ each contain one colour singlet. We see that the Dirac mass term coming from ${\bf 10}$ is just of the form $({\bf 2},{\bf 2},{\bf 1})$, while the Majorana mass term for $(\bar\nu_e)_L$ or $\nu_R$ coming from ${\bf 126}$ is a component of $({\bf 1},{\bf 3},{\bf \overline{10}})$. We must not use $({\bf 3},{\bf 1},{\bf 10})$, which would introduce an unwanted triplet violation of $SU^{\rm weak}_2$ and would give a mass directly to the left-handed neutrino. A possible danger is that radiative corrections might give rise to a large or uncontrollable term of that kind anyway, in addition to the term $m^2_{\rm Dirac}/m_{\rm Majorana}$, since the left-handed neutrino Majorana mass is not prohibited by a selection rule.

Such a selection rule exists in the $SU_5$ scheme, where an ungauged quantity that distinguishes ${\bf \bar 5}$ from ${\bf 10}$ and a gauged generator of $SU_5$ are simultaneously violated, preserving a linear combination, which is the baryon number minus the lepton number. The conservation of this quantity prohibits neutrino mass altogether. Here an unwanted massless spinless Goldstone boson is fed to an unwanted massless spin $1$ gauge boson to give a massive spin $1$ boson. This trick, which we have studied in connection with conserving baryon number (perhaps an obsolete idea now) can be applied whenever there is a reducible representation of $G$ for the fermions (or even for spinless elementary particles if there are some).

A further generalization of $SU_5$ for one family might make use of the lowest complex spinor representation ${\bf 27}$ of $E_6$, which breaks down to ${\bf 16} + {\bf 10 }+ {\bf 1}$ of $SO_{10}$. Here one would have to marry the new $SO_{10}$-singlet neutrino to the unwanted $(\bar\nu)_L$ of ${\bf 16}$, allowing them to share a huge Dirac mass, and one would have to do it in such a way as to leave the $\nu_L$ of ${\bf 16}$ with a small mass or none at all. At the same time one would have to assign high masses to the members of the ${\bf 10}$ of $SO_{10}$ in the ${\bf 27}$ of $E_6$, in order to get them out of the way, leaving just the fifteen fermions of the $SU_5$ scheme.

What we have seen from the example of one family is that a complex spinor representation, while it involves us in delicate questions of neutrino mass, does permit the description of left-handed fermions by a single irreducible representation of $G$ and in such a way that the asymmetry between the $SU_2^{\rm weak}$ assignments of particle and antiparticle is rather natural, while the whole system possesses an initial symmetry $C$ between left-handed particles and left-handed antiparticles, a symmetry that interchanges $SU_2^{\rm weak}$ and another $SU_2$.

It is also clear that in .such a scheme the dimensions of the representations that violate the symmetry, for example in the generation of fermion masses, tend to be large and that the arbitrary character of the violation scheme employing elementary Higgs bosons is strongly emphasized. It seems to us that one must hope for a situation in which, somehow, spontaneous symmetry violation is achieved dynamically.

Although we do not, of course, exclude the existence of some spinless elementary fields, provided they are not the arbitrary ones of the elementary Higgs boson method, we may look as an example at a theory with just gauge bosons and elementary spin $1/2$ fields and imagine what hypothetical dynamical spontaneous symmetry breaking would be like.

We would like to point out first that if the (say) left-handed fermion representation in such a theory is reducible, then ungauged quantum numbers arise that commute with the gauge group. When these are violated spontaneously, that necessarily leads to unwanted massless Goldstone bosons unless the trick described above is used and global conservation laws result. If all the irreducible representations are inequivalent, then such globally conserved quantities are Abelian and tolerable, but if there are equivalences among representations, as in the case of several families transforming alike, then an ungauged non-Abelian family group arises and that would have to be matched with an isomorphic subgroup of $G$ with resulting global conservation of a third isomorphic non-Abelian group relating the families. That would not agree at all with observation, and we conclude therefore that having united each family in an irreducible representation of $SO_{10}$ or $E_6$ we had better consider all the fermions as belonging to a single irreducible representation of the gauge group $G$.

This can be done in two different ways. Either we go to a higher-dimensional representation of the same group that we used for one family or else we enlarge the group and assign the fermions to a relatively low-lying representation of the bigger group. In the case of complex spinor representations, we could try, as an example of the first approach, the ${\bf 1728}$ of $E_6$, contained in ${\bf 27}\times {\bf 78}$. As examples of the second approach, we can take the lowest-dimensional complex spinors of larger groups, and the only larger groups possessing such spinors are $SO_{14}$, (lowest dimensional spinor ${\bf 64}$), $SO_{18}$ (lowest-dimensional spinor ${\bf 256}$), $SO_{22}$ (lowest-dimensional spinor ${\bf 1024}$), etc. We have studied both possibilities but we shall describe here the case of the lowest spinors of $SO_{4n+2}$.

A great deal of thought has been devoted to the question of what dynamical spontaneous symmetry breaking would be like for a theory containing elementary fields for gauge bosons and fermions only. Weinberg, Dimopoulos and Susskind, and various other theorists have drawn some important conclusions, including the following, which we specialize to the case of an irreducible fermion representation.

Symmetry reduction occurs through "condensations", that is non-zero vacuum expected values of operators that break symmetries. If the symmetry group of the kinetic energy is $H$ and if $G_1\subset   G$ and $H_1 \subset H$ are the subgroups left invariant by these condensations, then the generators of $G_1$ correspond to exact conservation laws and massless gauge bosons, those of $G/G_1$ to massive gauge bosons, those of $H_1/G_1$ to modified Goldstone bosons that acquire mass a a result of the gauge coupling and those of $(H/G)/(H_1/G_1)$ to approximate conservation laws, broken by the gauge coupling . The flavour-non-singlet pseudoscalar mesons would be modified Goldstone bosons, and the $PCAC$ condensation $<\bar {q} q>_{vac} \sim \Lambda^3$ presumably occurring in $QCD$ would contribute only $\sim ~ e\Lambda$ to the masses of the weak intermediate bosons. If one or more additional factors of the exactly conserved strong colour group exist (we prefer to call them primed colour, etc.), then these could have higher renormalization-groupvinvariant masses $\Lambda'$, etc., and a primed colour group with $\Lambda'\sim  10^3 ~GeV$ could give a condensation of fermions possessing primed colour that would account for the weak intermediate boson masses. Some of the corresponding pseudoscalar primed mesons would serve as effective Higgs bosons to be eaten by these intermediate gauge bosons. The mixing between these primed mesons and ordinary pseudoscalar mesons would be rather small. There would be no real ultra-violet fermion masses, but only medium-frequency or infra-red masses of order $\Lambda$ for quarks, $\Lambda'$ for fermions possessing primed colour, and so forth, and then masses obtained by sharing these medium-frequency masses through radiative corrections -these last would simulate ultraviolet masses up to fairly high energies.

A great deal of the algebraic behaviour of such symmetry-breaking schemes should be simulated by generalized non-linear $\sigma$-models. If those are embedded in linearized $\sigma$-models, then one has some connection with the algebraic properties of explicit Higgs boson theories.

An important question is whether the many condensations required for symmetry breaking in a unified theory can be explained by the strong long-range interactions that appear in the same theory. This is a problem, for example, in connection with any condensation leading to Majorana masses for the unwanted neutrinos.

Now let us return to the notion that $G$ might be $SO_{4n+2}$ with the left-handed fermions placed in the $2^{2n}$-dimensional complex spinor representation. Let us consider the example of $SO_{18}$, which evidently contains $SO_8 \times SO_{10}$. We can decompose the ${\bf 256}$ of $SO_{18}$ as $({\bf 8}_{sp},{\bf \overline{16}})+({\bf 8}_{sp}',{\bf 16})$ of $SO_8\times SO_{10}$, where ${\bf 8}_{sp}$ and ${\bf 8}_{sp}'$ are the two real inequivalent spinors of $SO_8$. We can now write $SO_8\supset   Sp_4\times SU_2$, where the vectorial octet ${\bf 8}_v$ of $SO_8$ can be made to give $({\bf 4},{\bf 2})$ of $Sp_4\times SU_2$ and ${\bf 8}_{sp}'$ of $SO_8$ likewise, while ${\bf 8}_{sp}$ of $SO_8$ gives $({\bf 1},{\bf 2}) + ({\bf 2},{\bf 1})$ of $Sp_4\times SU_2$. The ${\bf 256}$ of $SO_{18}$ then becomes $({\bf 1} ,{\bf 3},{\bf 16}) + ({\bf 5},{\bf 1},{\bf 16}) + ({\bf 4},{\bf 2},{\bf \overline{16}})$ of $Sp_4\times SU_2 \times SO_{10}$. If we interpret $Sp_4$ as a supplementary factor of the exactly conserved colour group, which becomes $SU_3^c\times Sp_4^{c'}$, and $SU_2$ as a gauged family subgroup of $SO_{18}$, we see that the only fundamental left-handed fermions without primed colour are three families of 16-dimensional spinors of $SO_{10}$, and we glimpse a possible agreement with experiment.

We note that $Sp_4^{c'}$, if there were no fermions would have the same renormalization group behaviour as $SU_3^c$ in lowest order, and we would need a special explanation for its reaching the strong-coupling regime at a much higher mass than $SU_3^c$. The differing fermion corrections might make a difference; so might the possibility that as we come down in mass from the unification region $SO_8$ remains undivided over a considerable interval before splitting into $Sp_4\times  SU_2$.

In the same way, $SO_{14} \supset SO_4\times SO_{10}$ and $SO_4$ is actually algebraically equivalent to $Sp_2 \times SU_2$; the ${\bf 64}$ of $SO_{14}$ decomposes into $({\bf 1},{\bf 2},{\bf 16}) + ({\bf 2},{\bf 1},{\bf \overline{16}})$ of $Sp_2\times SU_2\times SO_{10}$ and we would have two families lacking primed colour. Similarly, $SO_{22}\supset   SO_{12}\times SO_{10}$ and $SO_{12}\supset Sp_6\times SU_2$; the ${\bf 1024}$ of $SO_{22}$ decomposes into $({\bf 1},{\bf 4},{\bf 16}) + ({\bf 14},{\bf 2},{\bf 16}) + ({\bf 6},{\bf 3},{\bf \overline{16}}) + ({\bf 14'},{\bf 1},{\bf \overline{16}})$ and we would have four families lacking primed colour.

As far as representations giving fermion mass are concerned, we have the following situation: $$\begin{aligned}
SO_{10}:&&   ({\bf 16})_s^2= ({\bf 10})_A^5({\rm self-dual})+ ({\bf 10})^1={\bf 126}+{\bf 10};\\
SO_{14}:&&   ({\bf 64})_s^2= ({\bf 14})_A^7({\rm self-dual})+ ({\bf 14})^3_A={\bf 1716}+{\bf 364};\\
SO_{18}:&&   ({\bf 256})_s^2= ({\bf 18})_A^9({\rm self-dual})+ ({\bf 18})^1={\bf 24310}+{\bf 8568}+{\bf 18};\\
SO_{22}:&&   ({\bf 1024})_s^2= ({\bf 22})_A^{11}({\rm self-dual})+ ({\bf 22})^7_A={\bf 352,716}+{\bf 170,544}+{\bf 1540};
\end{aligned}$$ and so forth. It looks in each case as if the Majorana mass term comes from the highest-dimensional representation and the Dirac masses of the familiar fermions from the next-highest-dimensional one, if such a scheme is to work. The Dirac masses then obey an important constraint, which equates a function of the charged lepton mass matrix with the same function of the $Q = -1/3$ quark mass matrix. Since in each case the matrix is dominated, according to experience, by the highest mass, these two masses must be roughly equal, and for three families that explains the relation $m_b=m_\tau$ after renormalization.

The question is, of course, left open as to why the mass matrix for three families is so close to $$\begin{pmatrix} {c&0&0\cr 0&0&0\cr 0&0&0}\end{pmatrix}$$ for each kind of particle. With the families described as a triplet of $SU_2$ rather than a triplet of $SU_3$, that is rather mysterious, since it corresponds to a miraculous compensation of a scalar and a quadrupole term under $SU_2$. Under $SU_3$, of course, they would combine to form a ${\bf 6}$ of $SU_3$ and the approximate matrix above would correspond to the intervention of the component of ${\bf 6}$ invariant under the maximal little group $SU_2$. Unfortunately we are not dealing here with a family $SU_3$.

The Dirac masses of neutrinos and of $Q = + 2/3$ quarks would obey the same relation as that for the charged leptons and $Q = -1/3$ quarks. The Majorana masses of the neutrinos are also subject to a constraint if they come from the highest-dimensional representation for the fermion mass. Of course, the mass matrix for neutrinos is not easy to detect and at best requires delicate experiments that we shall describe elsewhere.

In summary, the idea of assigning left-handed spin $1/2$ fermions to a complex spinor representation of a gauge group $SO_{4n+2}$ (or conceivably $E_6$) has a number of attractive features, although some difficulties as well. As alternatives for an irreducible representation, we have, of course, the possibility of a real or pseudoreal representation, giving a vector-like theory in which all the known fermions must be accompanied by heavy partners that have weak interactions of opposite handedness; or a complex representation of a unitary group, which when irreducible generally leads to anomalies and thus to divergences, and is also rather hard to reconcile with observation.

If we suppose that the familiar quarks and leptons are really to be assigned to a complex spinor representation of a group $SO_{4n+2}$ or $E_6$, can we reconcile that idea with the notion that there is some truth in extended supergravity where the spin $1/2$ fermions are placed in a third rank antisymmetric tensor representation of $SO_N$?

We have looked, in collaboration with Jon Rosner, for an analogue of supersymmetry that might lead to a theory with assignments like ${\bf 1}$ of $E_6$ for $J_z= 2$, ${\bf 27}$ of $E_6$ for $J_z= 3/2$, adjoint ${\bf 78}$ of $E_6$ for $J_z=1$, and for $J_z= 1/2$ some representation contained in ${\bf 27}\times {\bf 78}$, like ${\bf 1728}$ of $E_6$. We have searched for the same kind of scheme using $SO_{10}$, and we have-even tried non-associative systems in an effort to find something that would work. So far we have had no success.

It seems likely anyway that if supergravity or some similar future theory is correct, then there must be only an indirect relation between the elementary fields of the theory and the particles that appear to us today to be elementary. If the known fermions behave, for a given handedness, like a complex spinor representation of $SO_{4n+2}$ or $E_6$, then the relation is not even that of a composite model. All or most of the familiar particles would have to correspond to particle-like solutions of the fundamental equations, with a different algebraic behavior from that of the fundamental fields.

In this talk we have only sketched the subject of complex spinor representations and related topics. Elsewhere we present a proper account of our own work and adequate references to the work of others.

We have also taken a rather schizophrenic approach, shuttling back and forth between extended supergravity on the one hand and a particular kind of unified Yang-Mills theory on the other. The ideas underlying the two approaches have to be compared more carefully.

ERRATUM:

On top of page 8, the expressions '$H_1/G_1$' and '$(H/G)/(H_1/G_1)$' should be interchanged

[^1]: Permanent address: 452-48 California Institute of Technology, Pasadena, California 91 125 U.S.A. Work supported in part by a grant from the Alfred P. Sloan Foundation.

[^2]: Work supported in part by the U.S. Department of Energy under contract No. DEAC- 03-79ER0068


abstract: Nonlinear dynamics of surface gravity waves trapped by an opposing jet current is studied analytically and numerically. For wave fields narrowband in frequency but not necessarily with narrow angular distributions the developed asymptotic weakly nonlinear theory based on the modal approach of leads to the one-dimensional modified nonlinear Schrödinger equation of self-focusing type for a single mode. Its solutions such as envelope solitons and breathers are considered to be prototypes of rogue waves; these solutions, in contrast to waves in the absence of currents, are robust with respect to transverse perturbations, which suggests potentially higher probability of rogue waves. Robustness of the long-lived analytical solutions in form of the modulated trapped waves and solitary wave groups is verified by direct numerical simulations of potential Euler equations.
author: V.I. Shrira$^1$; A.V. Slunyaev$^{2,3}$
date: 2024-10-02
title: Nonlinear dynamics of trapped waves on jet currents and rogue waves

*Introduction.—* In the last two decades there was a surge in interest in the phenomenon of extreme or rogue waves in various areas of physics, e.g. . In the most studied area of water waves in the ocean the main thrust of the studies was onto the search of mechanisms of rogue waves in the absence of currents . The modulational (or Benjamin-Feir, BF) instability of narrowband wave fields has been identified as the mechanism leading to formation of anomalously high waves and significant increase of their probability. This route has been most intensively studied and is most effective for one-dimensional (1-d) wave propagation; the theoretical and experimental modeling is the simplest. However, in reality the 1-d patterns are transversally unstable and, hence, short-lived; there is a dramatic difference in the probability of rogue events due to the BF instability for the strictly 1-d and 2-d wave propagation: the likelihood of rogue events is much higher for the 1-d propagation, moreover it completely vanishes for wave fields with the angular spectra width exceeding a certain threshold . On the other hand, it is known that the rogue waves are much more frequent on currents; the Agulhas current gained notoriety in this respect .

To explain increased probability of rogue waves on currents the prevailing approach exploits separation of scales between the typical wavelength and current, which leads to the WKB or ray description with a special consideration of caustics, see the literature reviews in .

With the focus on nonlinear dynamics, various versions of nonlinear Schrödinger equations (NLSE) were derived and analyzed under general assumptions of slow current, weak nonlinearity and narrow-banded spectrum, see and references there. In these works the BF instability was found to be strengthened for waves on adverse intensifying current. The triggering of the BF instability of a narrow band fields due to the intensification of the current was considered in . However there is an essential feature not captured by the existing NLS-type models: waves propagating upstream can be trapped by the current, there are multiple caustics. Such trapped waves have been observed on the Gulf Stream and were found to have considerably higher steepness than free waves on current . A more general and profound difficulty is that there is no technique enabling one to describe wave resonant interactions on currents; the waves refract on currents and hence vary in space, while the resonant interactions and, in particular, the resonance conditions have to be described in the wave-vector space. A new approach suggested in allows to overcome these obstacles for currents with a symmetry, in particular, for the common parallel jet currents. Instead of operating with rays we deal with the modes propagating on jet currents for which the standard nonlinear theory applies.

Here, based on the modal approach and weakly nonlinear asymptotic expansions, we derive equations governing one-dimensional wave evolution along the current; the transverse structure of the field is being provided by the modes. For a one mode we derive 1-d NLSE without the constraint of a too narrow angular spectra. In contrast to unguided waves, here the NLSE solitary wave type solutions are robust. The robustness of such wave patterns suggests dramatic increase of the probability of rogue waves. The predictions of the asymptotic model are validated by direct numerical simulations of the Euler equations.

*Modal representation for waves on jet currents.—* We consider wave motions on the free surface of an ideal incompressible fluid of unit density, waves are propagating along the $Ox$ direction on a given vertically uniform unidirectional steady current $U = \{U(y),0,0\}$. The motions are governed by the standard Euler and continuity equations in the domain occupied by the fluid $z \leq \eta$, where $\eta(x,y,t)$ is the water surface elevation; the water depth is assumed infinite for convenience. These equations are complemented with the standard boundary conditions for gravity waves: dynamic and kinematic boundary conditions on the surface and decay of velocities as $z \rightarrow -\infty$. We focus on the evolution of waves trapped by the current; trapped modes are selected by the boundary conditions stipulating lateral decay of velocities at $y \rightarrow \pm \infty$.

In the linear setting the problem formulated above was thoroughly examined in . Making use of uniformity of the problem with respect to $x$ and $t$, the Fourier transform may be applied with respect to these variables; then the complete linear solution has the form of a superposition of traveling waves propagating collinear to the current with some structure in $(y,z)$ plane, $w = \Re{ \sum_n{\hat{w}_n(y,z) \exp{\left( i \omega_n t - i k x \right)}} }$, where $w(x,y,z,t)$ is the vertical component of fluid velocity, index $n$ numerates the lateral modes. The modes $\hat{w}_n(y,z)$ are specified by the two-dimensional boundary value problem (BVP) $$\begin{aligned}
\label{2DBVP}
\frac{\partial ^2 \hat{w}_n}{\partial z^2} + \frac{\partial ^2
\hat{w}_n}{\partial y^2}  + \left( \frac{\Omega_n^{\prime\prime}}{\Omega_n}
- 2 \frac{\Omega_n^{\prime 2}}{\Omega_n^2} - k^2 \right) \hat{w}_n = 0 \,, \\
%\quad
\hat{w}_n|_{y \rightarrow \pm \infty}\rightarrow 0 \,,
\qquad
\hat{w}_n|_{z \rightarrow - \infty}\rightarrow 0 \,, \nonumber
\end{aligned}$$ with $\Omega_n(y) = \omega_n - k U$. Each mode is characterized by its cyclic frequency, $\omega_n$, and longitudinal wavenumber, $k$; we choose $k>0$ with no loss of generality.

The 2-d BVP () may be either solved numerically or reduced to a one-dimensional BVP using an asymptotic separation of variables . Here we adopt the second route, which assumes the following representation $\hat{w}_n = B_n Y_n(y) Z_n(z, y)$. Thus, the mode is specified by two real functions: $Z_n(z,y)$ determines the vertical structure which depends on $y$ slowly, $Y_n \partial Z_n / \partial y \ll dY_n / dy Z_n$, and is equal to one on the surface; $Y_n(y)$ determines the mode transverse dependence. Constants $B_n$ may be complex.

Asymptotic 1-d reductions of BVP () were derived in for two regimes: of 'weak' currents (compared to the wave phase velocity), and of 'broad' currents (compared to the longitudinal wave length). For the dominant wind waves and swell we are primarily interested in, all oceanic currents are weak in this sense. The corresponding 1-d BVP is of the Sturm - Liouville type, $$\begin{aligned}
  \label{BVP}
\frac{d^2 Y_n}{d y^2}  + 4 k^2 \left[ \frac{\omega_n}{\omega_g} - \left( 1 + \frac{kU}{\omega_g} \right) \right] Y_n =0 \,, \quad Y_n|_{y \rightarrow \pm \infty}\rightarrow 0\,.
\end{aligned}$$ Here $\omega_g = \sqrt{k g}$ denotes the frequency of linear gravity waves unaffected by the current. The wave frequency $\omega_n$ is the to-be-defined eigenvalue of the problem. For single-humped currents trapped modes require $kU < 0$ to exist; waves must run against the current.

*Weakly nonlinear model for modulated waves on jet currents.—* Now we concentrate on the regime where nonlinear interactions only between the trapped modes are essential. Employing a standard asymptotic procedure based upon a small parameter characterizing smallness of wave steepness: $\epsilon \sim k \max{|\eta|}$ or, equivalently, smallness of fluid velocities, $\epsilon \sim k \max{|w|} / \omega_g$, it is straightforward to derive a variety of evolution equations governing wave fields nonlinear dynamics; the resulting equations for constants $B_n$ now being slow functions of time and coordinate $x$ are determined by the choice of initial configurations of the field. In view of our interest in rogue wave formation we focus on narrowband spectra and consider wave fields with the longitudinal spectrum confined to $\epsilon$ vicinity of the carrier wavenumber, $k=k_c$.

We stress that jet currents profoundly modify the picture of wave resonances comparing to the freely propagating deep-water gravity waves; crucially, three-wave resonances become possible, which results in dynamical equations of the three-wave-interaction type. However, it can be shown that in the limit of weak currents, in which we are primarily interested, the four-wave interactions dominate. A key feature of the four-wave regimes is that to leading order wave dynamics is potential. Detailed analysis of all the regimes will be reported elsewhere. Here, on skipping the derivation, we provide and briefly discuss the new version of the NLSE which we obtain for the wave fields belonging to a single trapped mode with spectra narrow in longitudinal wavenumbers $k$; it reads: $$\begin{aligned}
 \label{NLS}
     - i \left( \frac{\partial A}{\partial t} +  V \frac{\partial A}{\partial x} \right) + \frac{\omega_n}{8k_c^2} \frac{\partial^2 A}{\partial x^2} + \kappa_n  \frac{\omega_n k_c^2}{2}  A |A|^2 = 0, \\
 %    V = \frac{\omega_n}{2k_c} + U(y),
    V=\int_{-\infty}^{\infty}{\left( \frac{k_c g^2}{2 \Omega_n^3} + U \right) Y_n^2 dy} \,,
     \quad
     \kappa_n = \frac{  \int_{-\infty}^{\infty}{  Y_n^4 dy}}{  \int_{-\infty}^{\infty}{Y_n^2 dy}}. \nonumber
\end{aligned}$$ Equation () describes evolution of the surface elevation, $\eta = \Re{ \left[ A(x,t) Y_n(y) \exp{\left( i \omega_n t - i k_c x \right)} \right]}$, $A$ is linked to the vertical velocity component as $B = i \omega_n A$; the subscripts for $A$ and $B$ are omitted for brevity. Sets of similar coupled equations appear when more than one mode is initially excited. The detailed derivation of () and its coupled generalizations will be reported elsewhere.

Equation () differs from the classical NLSE in the still water by the account for the Doppler frequency shift and a reduced nonlinear coefficient due to factor $\kappa_n$. By virtue of the Cauchy–Schwarz inequality $\kappa_n <1$. The NLSE () is of focusing type, hence it supports the BF instability and being integrable it admits a wide class of well studied exact solutions of variable degree of complexity . The basic solutions are uniform waves, localized solitary wave groups (envelope solitons) and breathers. Being a one-dimensional evolution equation, () provides a dramatic simplification of description of complicated dynamics of 2-d nonlinear trapped wave patterns.

*Nonlinear dynamics of trapped waves in simulations of the primitive equations.—* To verify the obtained asymptotic description a few key solutions of () are tested below by means of strongly nonlinear numerical simulations of the primitive hydrodynamic equations. The High Order Spectral Method (HOSM) to solve the potential Euler equations is adopted for the situations when the four-wave interactions dominate. The approach similar to is employed, where the terms responsible for current $U(y)$ have been introduced. The computational domain is periodic in both, $x$ and $y$ coordinates. The current is chosen to be close to $sech^2y$, it is taken to be periodical in $y$ with widely separated humps; it is specified as $U = U_{max} cn^2{(2K \frac{y}{L_y},s^2)}$, where $K(s^2)$ is the complete elliptic integral of the first kind, and $U_{max}=-2$ m/s and $s=0.9$ are used. The current varies from zero to $U_{max}$ at $y=0$; it is shown with arrows in Figs. , .

*1. Single trapped mode.* In the first experiment we verify the ability of trapped waves belonging to a single mode to propagate with no noticeable radiation in the fully nonlinear system. The initial condition has the form of a uniform train of 10 Stokes waves with the wavenumber $k=0.1$ rad/m and steepness $kH/2=0.15$ (where $H$ is the trough-to-crest wave height), modulated in the transverse direction according to the fundamental ($n=0$) mode function, $Y_0(y)$. The function $Y_0$ is found by solving () numerically (see details in ).

The evolution of an initially uniform wave train of trapped waves is simulated for about 80 wave periods . It propagates steadily with no evidence of significant radiation or structural deformation. The presence of some small-amplitude ripples is natural since the initial condition is not exactly an one-mode solution. The examination of the instantaneous wave height record gives some clues of two processes which lead to a slow decrease of the trapped wave height: i) about $10\%$ of wave height is lost during the first $\approx 10$ wave periods (we attribute this to the imperfect initial conditions); ii) a longer-term slow trend resulting in the total loss of about $2\%$ of energy over the simulated $80$ wave periods is apparently caused by interaction with noise. In other respects the train of trapped waves exhibits robustness.

The gravity wave angular frequency for the chosen initial condition is $\omega_g = 0.99$ rad/s. The nonlinear correction to the Stokes wave of moderate steepness, $kH/2=0.15$, adds extra $\approx 1\%$ to the frequency. The solution of the eigenvalue problem () yields $\omega_0 = 0.80$ rad/s, while the frequency spectrum of the simulated surface displacement has the maximum at $\omega = 0.80 \pm 0.01$ rad/s.

*2. Modulational instability of a single trapped mode leading to a rogue wave pattern.* A $5\%$ modulation along $Ox$ was applied to the wave train used in the previous simulation to initiate the modulational instability. Also, the reference simulation was performed when the current was set equal to zero, and the train had no modulation in the transverse direction. Supporting the NLSE prediction, the initial modulations grow in both cases with the maximal waves eventually reaching the breaking limit.

The initially modulated train of trapped waves undergoes further localization of wave energy, and the emerging large wave breaks at some instant, which leads to blowing up of the numerical iterations in time. The picture of surface elevation at the moment close to the wave breaking is given in Fig. . Due to the factor $\kappa$ in (), the evolution governed by the NLSE for trapped waves on currents is slower compared with the 1-d free gravity waves. Indeed, the curves of instantaneous maximal wave heights versus time in the simulations discussed above may be fitted onto each other, when the time is scaled with factor $0.65$ in the case of trapped waves. This value differs from $\kappa \approx 0.71$ calculated for the chosen profile of the current; the discrepancy is most likely due to the inaccuracy in prescribing the single-mode initial conditions for the trapped wave simulation.

*3. Solitary groups of trapped waves.* Envelope solitons are the most fundamental solutions supported by the 1-d NLSE, they represent asymptotics of the initial problem with generic localized initial data. The 2-d generalizations of the NLSE for deep water waves also admit planar envelope soliton solutions, but they are no longer asymptotics of the initial problem and are known to be unstable with respect to long transverse perturbations and, hence, relatively short lived.

To specify the initial condition for the strongly nonlinear simulations we use the exact analytic solution to equation () in the form of an envelope soliton, $A(x,t=0)=A_s sech{(\sqrt{2} k_c^2 A_s x)}$, with the transverse shape prescribed by the modal function, $Y_n(y)$; two cases ($n=0$ and $n=4$) were considered. In the runs the carrier wave has the same longitudinal wave number, $k_c = 0.1$ rad/m, though the intensity is smaller, $k_c A_s \approx 0.12$.

Fig.  shows the result of simulation of a solitary wave which corresponds to the fifth mode. The surface snapshot corresponds to the moment when the solitary group has passed the computational domain twice, that is about $40$ wave lengths. The longitudinal section of the solitary group is shown by lines above the maximum of the current and above the $x$-axis (red solid line); the transversal section of the group is shown in front of the surface above the $y$-axis (red solid line). These sections are compared with the corresponding sections of the initial condition (thin black curves). The amplitude of the solitary wave group ends up somewhat reduced compared to the initial condition; the radiated wave patterns are discernable in the figure.

The survived intense solitary wave group in Fig.  is an indication of robustness of such groups of trapped waves; although the solitary group produces some radiation each time it interacts with other wave patterns which exist in the simulation domain . In the course of evolution the group is slowly loosing energy; the total drop of the maximum wave height over $120$ wave periods is about $20-25\%$. The solitary group of waves belonging to the fundamental mode preserves energy better, though some radiation is observed as well . The decrease of maximum wave height over the same time of simulation is about $%\approx
10-15\%$. The patterns causing radiation by the groups can be viewed as an artifact of the imperfect choice of the initial conditions, or an unavoidable element of the imperfect reality the solitary waves are likely to encounter in nature.

*Concluding remarks—* In the context of oceanic waves the robust solitary groups of trapped waves found on jet currents represent a unique case of intense patterns of deep water gravity waves localized in both dimensions, which can propagate over long distances preserving their structure and energy. We stress that one-dimensionalization of wave dynamics occurs without assumption of a narrow angular spectrum; in the example shown in Fig.  the width of the spectrum is of order of the longitudinal wavenumber, $\Delta k_y / k_c \sim 1$. It makes relevant for oceanic conditions a huge corpus of theoretical and laboratory studies concerned with one-dimensional wave dynamics. As long as the trapped waves are described by the integrable 1-d NLSE to the leading order, all powerfull mathematical techniques and analytic solutions obtained since 1970-es may be applied (see e.g. and references therein). Correspondingly, the same well studied dynamics resulting in extreme waves for the planar waves (including rogue waves in optical fibers, where the NLSE and its solutions may work perfectly well ) is also expected for the trapped waves. As a particular conclusion, a higher likelihood of rogue wave events in the field of trapped waves should be expected. Exploiting integrability of the NLSE, elements of deterministic forecasting of oceanic rogue waves may be suggested.

The situations when several/many trapped modes and/or passing modes are present and interact the trapped waves require a dedicated study. Here we note that the phenomenon of one-dimensionalization of wave dynamics and its implications we discussed are not confined to water waves on currents; similar effects are likely in all branches of physics wherever there are guided waves.

The work was supported by the EC 7th Framework Grant PIIF-GA-2009-909389. AS acknowledges the support from RFBR via grants 11-02-00483 and 12-05-33087.