Measurement of identified pi^0 and inclusive photon v_2 and implication to the direct photon production in sqrt(s_NN) = 200 GeV Au+Au collisions

Abstract: The azimuthal distribution of identified pi^0 and inclusive photons has been measured in sqrt{s_{NN}} = 200 GeV Au+Au collisions with the PHENIX experiment at the Relativistic Heavy Ion Collider (RHIC). The second harmonic parameter (v_2) was measured to describe the observed anisotropy of the azimuthal distribution. The measured inclusive photon v_2 is consistent with the value expected for the photons from hadron decay and is also consistent with the lack of direct photon signal over the measured p_T range 1-6 GeV/c. An attempt is made to extract v_2 of direct photons.

\section{Acknowledgements}

We thank the staff of the Collider-Accelerator and Physics
Departments at BNL for their vital contributions.  We acknowledge
support from the Department of Energy and NSF (U.S.A.),
MEXT and JSPS (Japan), CNPq and FAPESP (Brazil), NSFC (China),
CNRS-IN2P3 and CEA (France),
BMBF, DAAD, and AvH (Germany),
OTKA (Hungary), DAE and DST (India), ISF (Israel),
KRF and CHEP (Korea), RMIST, RAS, and RMAE (Russia),
VR and KAW (Sweden), U.S. CRDF for the FSU,
US-Hungarian NSF-OTKA-MTA, and US-Israel BSF.


Blue Straggler Formation in Clusters

Abstract: Blue stragglers are thought to be formed from the merger or coalescence of two stars, but the details of their formation in clusters has been difficult to disentangle. We discuss the two main formation mechanisms for blue stragglers (stellar collisions or mass transfer in a binary system). We then look at the additional complications caused by the stars living in the dynamically active environment of a star cluster. We review the recent observational and theoretical work which addresses the question "which mechanism dominates?" and conclude that the most likely answer is that both mechanisms are at work, although with different importance in different environments and at different times in the cluster lifetime. We finish with a short discussion of some avenues for future work

\section{search for formation mechanisms}

Blue stragglers are stars which are bluer and/or brighter than the
main sequence turnoff in clusters. They are found in all resolved
stellar populations, including globular and open
clusters, dwarf galaxies \cite{2007MNRAS.380.1127M}, and even
the field \cite{2010arXiv1002.3009D}. They are main sequence stars, lying between the zero-age main sequence and the giant branch in clusters. Typically they do not lie on a tight sequence, but instead fill much of that region of the colour-magnitude diagram. Where we can identify companions, we know that blue stragglers often occur in binary systems or even in triples
\cite{2009Natur.462.1032M, 2010arXiv1008.4347T}.  Typically there are
a few tens of blue stragglers in open clusters, and up to a few
hundred in globular clusters. While they only represent 1\% or less of
the stellar population in any cluster, they are easily identified and
provide a way to probe the history of the cluster and/or its binary
population, depending on our understanding of how this unexpected
population was formed.

Any formation mechanism for blue stragglers must find a way to either
increase the mass of a normal (low-mass) cluster main sequence star,
or delay the evolution of a higher-mass star such that it stays on the
main sequence longer than its peers. Currently, the two main formation
mechanisms which are thought to create blue stragglers are direct
stellar collisions, or mass transfer in a binary system. We discuss
each of them below, and then look at the importance of the cluster
environment.

\subsection{Stellar Collisions}

In the simplest form of this scenario, two previously unrelated stars
collide as they are moving through the cluster. A fairly dense
environment is needed, but calculations of collisional cross sections
in globular \cite{1976ApL....17...87H} and open
\cite{1989AJ.....98..217L} clusters have shown that collisions are
common enough in these environments to produce approximately the right
number of blue stragglers. In some clusters, it is expected that every
star has undergone at least one close encounter. There is also
evidence that another exotic population, cataclysmic variables, is a
collision product and they are over-abundant in clusters
\cite{2006ApJ...646L.143P}. Blue stragglers are also concentrated in
the cores of clusters \cite[e.g.][]{1992A&A...262...63L}, suggesting that
they are more massive than the average star, but also indicative that
they are formed in the densest environments.

Hydrodynamic simulations of stellar collisions
\cite[e.g.][]{1996ApJ...468..797L} show that for the relatively gentle
collisions which occur in globular clusters, the collision product is
a star-like object. Only a few percent of the total mass of the two
parent stars is lost during the collision, and the rest settles into a
spherical, hydrostatic object. The collision product is typically not
fully mixed, but instead the cores of the parent stars end up at the
core of the collision product \cite{1996ApJ...468..797L}. Subsequent
stellar evolution calculations of the collision product
\cite{2001ApJ...548..323S, 2008A&A...488.1017G} have shown that after
some initial thermal-timescale relaxation, the star evolves relatively
normally. One such example is shown in Figure
\ref{EvertCollision}. The red line is the evolutionary track of the
product of a collision between a 2.4 and 1.9 M$_{\odot}$  star, and the blue
line shows the evolution of a normal star of 4.3 M$_{\odot}$ . The initial
thermal relaxation of the collision product has not been shown. In
general, collision products are quite similar to normal stars of the
same mass. However, their main sequence lifetime is generally shorter
than that of the normal star. Since the parents of the collision
product have been evolving for some time before the collision, they
have converted some of their hydrogen to helium, and so the collision
product has a smaller hydrogen reservoir and a shorter main sequence
lifetime. Their position in the colour-magnitude diagram is therefore
usually above the zero-age main sequence, reflecting their advanced
stage of evolution when they are formed. However, the collision
products do typically spend considerable time in the appropriate
region of the colour-magnitude diagram so that we would expect to see
them as blue stragglers.

\begin{figure}
  \includegraphics[height=.4\textheight]{cmd_3400_08_24.eps}
  \caption{Evolution of the product of a collision between a 2.4 M$_{\odot}$
    and a 1.9 M$_{\odot}$  star (red), compared to the evolutionary track of
    a normal 4.3 M$_{\odot}$  star (blue). Models courtesy of E. Glebbeek,
    based on the calculations in Glebbeek \& Pols\cite{2008A&A...488.1017G}.}
  \label{EvertCollision}
\end{figure}

\subsection{Binary Mass Transfer}

The other plausible mechanism for blue straggler formation involves
mass transfer in a binary system. We know that binaries exist in all
stellar populations, and we also know that if the orbits are small
enough, mass will be transfered from one star to another. Binary mass
transfer has been invoked to explain a wide variety of unusual stars,
starting with the Algol system and encompassing objects at all
evolutionary stages. Blue stragglers are seen in all well-studied
stellar populations, including sparse clusters and the field. Stellar
collisions are not expected to contribute to blue stragglers in these
environments, and therefore binary mass transfer provides an obvious
and likely alternative formation mechanism.

If the accreting star is a main sequence star, then the product will
be a blue straggler. Not all of the donor star needs to become part of
the blue straggler -- only enough mass to significantly change the
accretor's position in the colour-magnitude diagram needs to be
transferred. For relatively high-mass accretors near the turnoff, this
could be only a few tenths of a solar mass. The other extreme is a
mass transfer event in which the entire donor star is accreted onto
the blue straggler. Some authors refer to the latter as ``binary
coalescence'' or ``binary merger''. For the purposes of this review,
we will consider any amount of mass transfer as belonging to this
mechanism.

The position of the blue straggler in the colour-magnitude diagram and
its main sequence lifetime depend heavily on the masses and
evolutionary states of the two stars involved. The initial binary
orbit also determines how much mass will be transfered, and when. The
parameter space of binary mass transfer blue stragglers has not yet
been fully explored. The most comprehensive studies to date follow
binaries which are appropriate for producing blue stragglers in old
open clusters like M67 \cite{2006A&A...455..247T,
  2008MNRAS.384.1263C}. Both case A and case B mass transfer have been
studied. Again, the evolutionary tracks for these objects spend a
significant amount of time in the blue straggler region of the
colour-magnitude diagram. Interestingly, many of the mass transfer systems spend
the majority of their time at or near the equal-mass binary line (0.75
magnitudes above the zero age main sequence), regardless of the actual
mass ratio of the system.

The present-day appearance of binary mass transfer products will
depend strongly on when the mass transfer began, whether the mass
transfer was conservative, and whether the mass transfer was
complete. Incomplete mass transfer will result in a blue
straggler in a binary system, and the nature of that companion star
can tell us a lot about the initial configuration. For example, a main
sequence companion means that case A mass transfer was the cause,
while a helium white dwarf points to a giant branch donor instead.

\subsection{Influence of the cluster environment}

In the discussion above, we have considered the effects of stellar
dynamics and of binary evolution separately. In clusters, however,
this simplification is invalid. We know that all clusters have a
population of binaries, and these binaries can interact dynamically
with each other and with single stars. These interactions will modify
our expectations for the blue straggler populations in clusters in a
number of ways.

Dynamical interactions will modify the binary population in
clusters. The standard wisdom can be summarized as ``hard binaries
harden'' \cite{1975MNRAS.173..729H}. In other words, during long-range
encounters with other objects in the clusters, binaries whose orbital
energy is larger than the average kinetic energy of stars in the
cluster will have their orbital energy increased. Soft binaries, on
the other hand, have their orbital energy decreased by encounters,
even to the point of disruption of the binary completely. What
implications will this have for those binaries which would have formed
blue stragglers through mass transfer? In order to produce a blue
straggler that we identify at the current time, mass transfer needs to
have happened in the somewhat recent past.  Also, the accretor needs
to be a main sequence star of the appropriate mass. Otherwise, the
blue straggler will have had time to leave the main sequence and
evolve off to the giant branch. If the distribution of binary orbital
energies (and hence separations) is being modified by interactions,
then the number of appropriate binaries which would have undergone
mass transfer at the appropriate time will have changed. Some systems
will be hardened to the point where mass transfer will happen earlier
and we will not see them as blue stragglers at the current
time. Other, wider systems, which would not have become blue
stragglers, may be hardened enough that those objects will contribute
to the blue straggler population. However, we must also consider the
evolutionary state of the donor when determining which objects will
become blue stragglers -- the initial orbital separation is not
enough information.

In addition to long-range, gentle encounters, binaries in clusters
undergo close encounters with single stars or other binaries which
result in much more violent outcomes. Binaries can be broken apart
entirely, which means that those objects cannot later produce blue
stragglers. These close encounters are typically resonant encounters
\cite{1983ApJ...268..319H}, meaning that the individual stars spend a
long time involved in a complicated dance around each other. Exchange
interactions are quite common, and the typical outcome is that the
heavier stars are exchanged into the resultant binary. Therefore, the
stellar evolution timescale of this new binary is shorter than the
original one. Even if the new binary does produce a blue straggler, it
will do so sooner and the blue straggler is likely to have evolved
away to the giant branch before the present day
\cite{2006A&A...459..489D}.

The major effect of the cluster environment is to blur the distinction
between the two formation mechanisms for blue stragglers. Direct
collisions are much more likely between two {\it binary} stars rather
than two single stars, since the cross section of a binary star goes
as its semi-major axis instead of the stellar radius. During resonant
interactions, two or more of the stars involved could run into each
other. Even though the velocities are larger than in the single star
case (because of the addition of the orbital velocities of the
binaries), the collisions are still sufficiently gentle that our
models of collision products still apply. Therefore, identifying an
individual blue straggler as a collision product or a binary product
may not be an easy, or even a valid, distinction.

\section{Which Mechanism Dominates?}

Standard practice in science, when presented with two likely
explanations for any phenomenon, is to figure out which explanation is
correct. In astronomy, observations are made in many wavelengths and
modes; theorists produce detailed models of individual objects and of
whole populations. Comparisons are made, and after a while, one
explanation emerges triumphant over the other. However, in the
question of blue straggler formation, we are not yet at the final
stage of this story. In fact, it seems that the question we need to
answer is not ``which mechanism dominates'' but rather ``does one
mechanism dominate?'' And, in fact, we many need to look for different
dominant mechanisms in different environments in the universe and
perhaps even in different regions in the same star cluster.

As we have discussed above, both collisions and mass transfer can
reproduce observed properties of blue stragglers. These mechanisms
create main sequence stars, more massive than the turnoff, with
lifetimes that are long enough that we expect to see significant
numbers of these objects in present-day clusters. Therefore, in order
to distinguish between the different formation mechanisms, we must
look at the entire population of blue stragglers in clusters. Can either
mechanism correctly predict the distribution of blue stragglers in the
cluster, and in the colour-magnitude diagram? Which mechanism can
correctly predict the number of blue stragglers found in each cluster,
or predict the sizes of populations as a function of global cluster
properties? Thanks to both increasingly detailed observations of blue
stragglers, and to more and more models, we are starting to be able to
answer these questions.

An early comparison of collisional blue straggler populations to
observations of cluster populations was done for 6 globular clusters
\cite{2003ApJ...588..464F} with a range of central densities. The
authors compared the luminosity functions of blue stragglers with
collisional models, as well as their distribution in
temperature. Because these were HST observations, only core blue
stragglers were included in the comparison. The collisional models
could fit the data reasonably well for most clusters, but not for the
lowest-density cluster in the sample (NGC 288). This was expected
since the collision rate in that cluster is quite low. Most of the
clusters have a population of blue stragglers which are redder than
the collisional models can explain, as it is very difficult to get a
star to remain in the Hertzsprung gap for any substantial length of
time.

A similar study of the blue straggler population throughout the nearby
cluster 47 Tucanae \cite{2006ApJ...650..195M} made use of HST and
ground-based data, thereby covering the entire extent of the blue
straggler population in this cluster. In this case, the blue straggler
population was assumed to be collisional, and the models allowed for
two free parameters: the mass function of the main sequence stars
which went into the collisions, and the duration over which the
collisions could occur. The collisional models fit the observations
very well. Interestingly, the blue straggler models implied that the
mass function in the core of the cluster is very heavily weighted
towards high mass stars, and that it becomes shallower towards the
cluster edge. This result is consistent with the direct measurement of
the main sequence mass function as a function of radius in 47
Tuc. Therefore, using blue stragglers as tracers of the dynamical
state of the cluster seems plausible, which is really only possible if
collisions contribute significantly to the formation of these objects.

However, around the same time, some completely contradictory evidence
was published. By comparing the fraction of blue stragglers in
clusters to the current collision rate, Piotto et
al. \cite{2004ApJ...604L.109P} showed that there was not the expected
correlation. In fact, there seemed to be an {\em
  anti}-correlation. That paper concluded that collisions could not be
the dominant mechanism for producing blue stragglers. A consistent
selection of core blue stragglers, and comparison to many cluster
properties (e.g. velocity dispersion, central density, etc), concurred
with that result \cite{2007ApJ...661..210L,
  2008ApJ...678..564L}. Finally, a robust correlation with a cluster
property was found \cite{2009Natur.457..288K}: The number of core blue
stragglers is proportional to the cluster's core mass. The authors
conclusively proved that the number of blue stragglers could not be
proportional to the collision rate. In addition, they showed that the
blue straggler number was not linearly proportional to the core mass,
but instead goes as $\sim M_{c}^{0.5}$. If the binary fraction in
globular clusters also depends on the (core) mass as $\sim
M_{c}^{-0.5}$, then it seems reasonable that the number of blue
stragglers in globular cluster cores is proportional to the number of
binaries in the core, $f_b M_c$. In fact, preliminary evidence does
show approximately this dependence of binary fraction on cluster mass
\cite{2007MNRAS.380..781S, 2008MmSAI..79..623M}. In addition, a
correlation between blue straggler populations and binary fraction was
found for some sparse globular clusters
\cite{2008A&A...481..701S}. Therefore, binary stars must be involved
in the formation of blue stragglers, even in the cores of clusters.

Additional evidence supporting a binary evolution for blue stragglers
comes from detailed observations of surface chemical abundances in
these stars.  Ferraro {\it et al.} \cite{2006ApJ...647L..53F} did a spectroscopic survey
of blue stragglers in 47 Tuc, and found some evidence for mass
transfer -- a subset of the blue stragglers were depleted in both
carbon and oxygen compared to normal stars in the cluster. One way to
do this is to have CNO-processed material from the inner regions of a
star be dumped onto the surface of a companion, where it is now
visible. Unfortunately, a similar study of blue stragglers in M4
\cite{2010ApJ...719L.121L} found no evidence of such a population so
this signature seems not to be ubiquitous.

So, collisions can explain the observations of blue stragglers and can
predict cluster properties....except that blue stragglers cannot be
formed through collisions and must have a binary origin. The picture
is not becoming clearer. And, indeed, we have evidence that {\em both}
mechanisms are going on in the same cluster. First, we have the
long-known ``bimodal distribution'' of blue stragglers in globular
clusters. HST plus ground-based observations of a number of clusters,
starting with M3 \cite{1997A&A...324..915F}, have shown that the blue
straggler population is largest in the centre of the cluster. Then
there is a significant drop of blue straggler numbers with increasing
cluster radius, and after some point the numbers go back up again,
although not usually to values as high as the central value. The best
explanation for this distribution was presented by Mapelli et
al. \cite{2006MNRAS.373..361M}. A blue straggler population which
consists of a collisional population and a binary population was
allowed to react dynamically with the underlying cluster stars for the
typical lifetime of a blue straggler. The collisional population was
initially found only in the core of the cluster, where the binary
population was initially spread evenly throughout the cluster. These
binary blue stragglers started to fall into the core, because of mass
segregation. But because of the finite lifetime of blue stragglers,
not all of them have a chance to reach the core in the lifetime of the
cluster. The radius of the minimum of the blue straggler population is
called the ``zone of avoidance'', and is where the dynamical friction
timescale is short enough that the blue stragglers can have reached
the core. Outside of that radius, the blue stragglers remain
more-or-less at their initial radii. As shown in Mapelli {\it et al.}
\cite{2006MNRAS.373..361M} and later papers
\cite{2007ApJ...663.1040L,2008ApJ...679..712B}, the predicted zone of
avoidance matches the minimum blue straggler position within a factor
of 2 for most clusters.

The second strong piece of evidence for two formation mechanisms at work in the
same cluster comes from HST observations of M30
\cite{2009Natur.462.1028F}. In the core of this cluster, two distinct
sequences of blue stragglers can be seen (see figure
\ref{FerraroNature}). The bluer sequence lies slightly above the ZAMS,
with the brightest of these blue stragglers further from the ZAMS than
the fainter stars on this sequence. It is extremely well-modeled with
an isochrone of collision products, approximately 1 Gyr after the
collisions have occurred. The red sequence consists of stars which are
at least 0.75 magnitudes brighter than the ZAMS, and some are
substantially brighter than that and lie in the Hertzsprung gap. These
stars are consistent with the mass transfer models from
\cite{2008MNRAS.384.1263C}. The redder blue stragglers are more
centrally concentrated than the blue ones, which has interesting
implications for the binary properties of the two progenitor
populations. In another Gyr or so, normal stellar evolution will blur
the distinction between the two sequences. The fact that we see two
sequences at the present time suggests that the formation of this
population of blue stragglers was both recent and short-lived. It is
interesting to consider that this may be an indication of core
collapse or some other clearly-defined dynamical event.

\begin{figure}
  \includegraphics[height=.5\textheight]{FerraroNature2010.eps}
  \caption{Colour-magnitude diagram of M30, from HST WF/PC2
    photometry. The two sequences of blue stragglers are shown as red
    and blue dots. The solid red lines are a 0.5 Gyr isochrone
    (equivalent to the ZAMS) and a 13 Gyr isochrone, which is
    appropriate for this cluster. The black lines show collisional
    isochrones, 1 and 2 Gyr after the collision events occurred. The
    dashed red line is 0.75 magnitudes brighter than the ZAMS. Figure
    taken from Ferraro {\it et al.} \cite{2009Natur.462.1028F}}
  \label{FerraroNature}
\end{figure}

Finally, evidence that we must include both (or hybrid) formation
mechanisms comes from ``kitchen-sink'' dynamical models of globular
cluster evolution. This result was first presented in the context of
N-body models of the old open cluster M67 \cite{2005MNRAS.363..293H},
where the blue stragglers were found to have very complicated
formation histories involving both dynamics and binary
evolution. Recently, we have been investigating the blue straggler
populations of approximately 100 globular cluster models using the
Northwestern Monte Carlo code (Sills, Chatterjee \& Glebbeek, in
preparation). Our preliminary results are shown in figure
\ref{chatterjee}. Here we show the observed trend of core blue
straggler number with core mass (open circles), and the results of our
simulations overplotted (squares). These models, which reproduce the
observations very nicely, require a cluster binary fraction of
5-10\%. Cluster with no initial binaries simply do not produce enough
blue stragglers, especially in the high-mass clusters. We can also
look at the formation histories of all our modeled blue stragglers,
and we find that very few of them are purely collision products or
purely binary mass transfer products. Almost all of them have both
dynamical encounters and binary evolution in their history. We can also
distinguish between blue stragglers formed in a binary encounter from
those formed from direct collision, and we find that the binary
encounters significantly outnumber both direct collisions and binary
mass transfer.

\begin{figure}
  \includegraphics[height=.5\textheight]{Nmcore.eps}
  \caption{Number of blue stragglers in the core of globular clusters,
    vs. core mass. Observations are shown as open circles, model
    calculations are dark squares. The models require a binary
    fraction of 5-10\%, and contributions from both collisions and
    binary evolution, to match the observations.}
  \label{chatterjee}
\end{figure}

\section{Future Work and Open Questions}

We are converging on an explanation for the formation of blue
stragglers in clusters. Unfortunately it seems to be a complicated
explanation, involving multiple formation mechanisms which act with
different weights in different environments and at different times in
the cluster evolution. However, we are certainly making progress with
a combination of detailed observations and more sophisticated models
of both individual blue stragglers and cluster populations. As always,
there are still a number of additional complications and
considerations that we still need to consider. The following is a
combination of a ``wish list'' and a nod to additional physics we may
still need to include in our models.

For the most part, we have completely ignored the effects of any
system more complicated than binary stars.  However, there has been
some promising research recently pointing to the importance of triple
systems, particularly for blue straggler populations in old open
clusters \cite{2009ApJ...697.1048P, 2010arXiv1009.0461L}.  Modeling the
dynamical effects of triples is extremely computational intensive in
N-body codes, and they are completely neglected in currently Monte
Carlo codes, but we may have become cleverer if we really wish to
investigate all possible formation mechanisms of blue stragglers and
other stellar exotica in clusters in the context of full-blown cluster
models.

In the mid-90's, there was a concerted effort to understand
single-single and single-binary encounters. Calculations of encounter
outcomes and cross sections were automated and parameterized, and so
we have a fairly good understanding of these kinds of
encounters. However, because the parameter space of binary-binary
encounters was so large, and because it is much harder to determine
when an encounter is ``finished'', we know much less about the
outcomes of higher order encounters. And, of course, very little work
has been done on triples and higher order multiples, as discussed
above. Special-purpose codes such as FEWBODY
\cite{2004MNRAS.352....1F} have been built to follow encounters of a
small number of stars, and could be harnessed for this purpose. Brute
force, however, is unlikely to cover the interesting parameter space
in a useful amount of time, and so the problem of being selective in
initial conditions is likely to remain.

I have not mentioned the rotation rates of blue stragglers in
clusters, and the information that those measurements can give
us. There are very few measured rotation rates of blue stragglers,
particularly in globular clusters. The situation is rather better in
open clusters, where the fields are sparser and spectroscopic studies
have been going on for longer (e.g. for characterization of the binary
populations in those clusters\cite{2009Natur.462.1032M}). The
indications are that blue stragglers do not rotate extremely rapidly,
but that they rotate more rapidly than normal stars of the same
temperature (typically 10-40 km/s for blue stragglers, compared to
2-10 km/s for normal stars). Have they been spun up by collisions or
mass transfer? Have they had a chance to spin down from a very
significant rotation rate? Can the rotation rates tell us anything
about their time of formation? Does the rotation mix the stars
significantly, so that their evolutionary tracks are different from
non-rotating stars?  Models of the angular momentum evolution of both
collision products e.g. \cite{2005MNRAS.358..716S} and binary mass
transfer products are required before we can use these pieces of
observational evidence to understand anything more about blue
straggler formation in clusters.

We have detailed evolutionary models of some binary mass transfer
products. However, the (very large) parameter space has not yet been
sufficiently explored. We need to understand enough about the
evolution of mass transfer products that we can approximate their
evolution in binary evolution codes like BSE
\cite{2002MNRAS.329..897H}. This has been done for collision products,
to some extent -- the Make Me A Star code (MMAS)
\cite{2002ApJ...568..939L} will produce a collision product from two
parents, without the need for the detailed hydrodynamic
calculation. The subsequent evolution of this collision product can be
calculated using a stellar evolution code, or the lifetime and CMD
location can be estimated based on the core hydrogen fraction of the
product. We are not yet at this stage for binary mass transfer models.

And finally, it would be of great interest to a large population of
researchers to know the binary properties in globular clusters. Does
the binary fraction depend on the mass of the cluster, as suggested by
preliminary work? Does it depend on the location in the cluster? The
answer to this question seems to be ``yes'', based on some scant
observations and a number of dynamical models. Does it depend on which
stellar population one observes? Of particular interest to
understanding the formation of blue stragglers, we're looking for two
answers to this question. One, what is the binary fraction of main
sequence-main sequence binaries, which could form blue stragglers
later? Two, what is the binary fraction of the blue stragglers
themselves? And, of course, any information we can glean about the
mass ratio, semi-major axis, and eccentricity distributions of any of
these binaries will be crucial for many of us. Will there be a
multi-object spectrograph on a telescope with sufficient resolution to
disentangle the crowded conditions of globular clusters any time soon?
I know of no such instrument in the planning stages, so we may be
using approximations and indirect methods for a long time.

\bibliographystyle{aipprocl}

Solar and Stellar Astrophysics (astro-ph.SR); Astrophysics of Galaxies (astro-ph.GA)

Blue Straggler Formation in Clusters


On the Moat-Penumbra Relation

Abstract: Proper motions in a sunspot group with a delta-configuration and close to the solar disc center have been studied by employing local correlation tracking techniques. The analysis is based on more than one hour time series of G-band images. Radial outflows with a mean speed of 0.67 km s^{-1} have been detected around the spots, the well-known sunspots moats. However, these outflows are not found in those umbral core sides without penumbra. Moreover, moat flows are only found in those sides of penumbrae located in the direction marked by the penumbral filaments. Penumbral sides perpendicular to them show no moat flow. These results strongly suggest a relation between the moat flow and the well-known, filament aligned, Evershed flow. The standard picture of a moat flow originated from a blocking of the upward propagation of heat is commented in some detail.

\section{Introduction}

Sunspots, as complex magnetic structures embedded in a convective plasma,
show many active and changing features on multiple scales. Convective flows and
large-scale plasma circulation plays an important role in the dynamics and
evolution of solar active regions \citep[see, e.g.,][]{schrijver2000}
The granular convective pattern surrounding
sunspots is perturbed by the presence of magnetic elements known as Moving
Magnetic Features (MMFs), magnetic elements that move radially
outwards through an annular cell called the "moat" \citep[][for a recent view]{sheeley1972,harvey1973,hagenaar2005}
This moat flow also carries away
granules as inferred using local correlation techniques \citep[LCT,][]{bonet2005}.
It has been suggested that the moat could be a supergranule whose
center is occupied by a sunspot with a typical cell-scale of up to $10^4$ km
\citep{meyer1974}. The sunspot would act as a blocking agent for the
upward propagation of heat from below. The excess temperature and
pressure generated in this way has been proposed as the origin of the
moat flow \citep[see also][]{nye1988}. An interesting property of this supergranular cell
is that, in deeper layers, it would have a component directed towards the spot
that could help stabilizing the magnetic structure.
While this viewpoint has been prevalent in sunspots physics, the real nature of the
moat flow is not well understood. Already \cite{vrabec1974} pointed out that in an irregular
sunspot, the moat was observed only in the sector with a well developed penumbra.
The sunspot sector that had no penumbra, was an area where pores were being advected
(inflowing) {\it towards} the main spot and, thus, displaying no moat-like flow.

\begin{figure}
\centering
\begin{tabular}{c}
\includegraphics[width=1.\linewidth]{./fig1.eps}
\end{tabular}
\vspace{-1mm}
\caption{G-band image of the field studied ($57\farcs8 \times 34\farcs4$),
reconstructed with the MOMFBD technique (see the text). The arrow in the upper
left corner points to the solar disk center.
\label{F:1}}
\end{figure}
This old result already suggests a link between the presence of a penumbra and the
moat flow, but no studies exist confirming such a relationship.
It is interesting to point out that pores
are not seen to be surrounded by a moat (whereas mature sunspots
always are). \cite{sobotka1999} studied (using LCT) six solar pores (including one with
6 arcsec diameter) and found the surrounding motions to be dominated by mesogranular
flows, as found elsewhere in the photosphere. No trace of a moat flow was observed.
It is clear that any possible heat blockage
by the pores magnetic body was unable to generate a moat flow.

The present paper studies the relation between moat flows and the
existence of penumbrae.
It is based on an excellent $\sim$79 minutes time series of a
sunspot group with $\delta$-configuration. The $\delta$ configuration
is ideal for this study because it displays complex umbral regions
harboring fully developed penumbra on one side and none on the other.
The long duration of these series,
besides the high and stable quality throughout the entire period, substantially
improved after image reconstruction (see Fig.~\ref{F:1}), make of this material an
excellent data set to study the morphology and dynamical behaviour of sunspots
and their surroundings.

\begin{figure*}
\centering
\begin{tabular}{ll}
\includegraphics[width=.8\linewidth]{./fig2.eps}
\end{tabular}
\vspace{-6mm}
\caption{Map of horizontal velocities in the entire FOV (71 min average). The
black contours outline the borders of umbrae, penumbrae and pores. The white
lines delineate the LOS neutral lines infered from a magnetogram. The
coordinates are expressed in Mm. The length of the black bar at coordinates
(0,0) corresponds to $0.4$ km s$^{-1}$. The background represents the average
image of the G-band series.
\label{F:2}}
\end{figure*}
\vspace{3mm}

\section{Observations and data processing}

The NOAA Active Region 10786 was observed on July 9, 2005 at the Swedish
1-meter Solar Telescope \citep[SST,][]{scharmer2003a} equipped with Adaptive
Optics (AO, Scharmer et al. 2003b), under collaboration in the International
Time Program. This complex region, corresponding to a $\delta$-configuration,
was placed at heliocentric position $\mu = 0.9$. A dichroic beamsplitter separates
setups in the red and the blue. In the blue beam, two simultaneous time
sequences of high-resolution images were taken in wavelength bands at the
G-band ($430.5 \pm 0.54$ nm) and nearby continuum ($436.3 \pm 0.57$ nm,
G-cont). The sequences span for more than one hour, from 7:47 UT to 9:06 UT,
following the evolution of the sunspots. The images were acquired using Kodak
Megaplus 1.6 CCD cameras with a 10-bits dynamical range and 1536 $\times$ 1024
pixels. The pixel size was $0\farcs041$ square. Real-time corrections with the
AO and further post-processing techniques rendered up images near the
diffraction limit. Standard procedures for flatfielding and dark-current
subtraction were applied before restoration.

Post-processing for image restoration was performed by employing the
Multi-Object Multi-Frame Blind Deconvolution (MOMFBD) method \citep{noort2005}. The observational strategy to apply this
technique consisted in taking G-band images and simultaneous G-cont phase
diversity image-pairs \citep{gonsalves1982,paxman1992,lofdahl1994} by using an optical configuration with three channels. Sets of
about 18 images per channel (i.e. 3 $\times$ 18 images) were combined to
produce a pair of simultaneous G-band and G-cont restored images. Following
this procedure, two time series of reconstructed images were produced. Figure \ref{F:1} shows one of these images in the G-band.

The restored images were de-rotated to compensate the diurnal field rotation
and rigidly aligned using cross-correlation. Additional image processing
consisted of destretching and p-mode filtering (threshold phase velocity $4$
km s$^{-1}$). The final product was two movies (G-band and G-cont) of 428
frames each, spanning over 71 min with a cadence of 10.0517 s, and covering a
field-of-view (FOV) of $57\farcs8 \times 34\farcs4$ (see movies at the web
site: \emph{http://www.iac.es/galeria/svargas/movies.html}).

In the red beam a tunable birefringent filter \citep[SOUP,][]{title1986}
working at FeI 6302\AA~in combination with liquid crystal retarders produced
longitudinal magnetograms. A beamsplitter in front of the SOUP filter splitted off
10\% of the light to obtain phase diversity image-pairs that combined with the
narrow-band images of the SOUP, allowed reconstructions MOMFBD to finally
provide magnetograms with $\sim 0\farcs25$  of resolution (see Fig.~\ref{F:4}).

\begin{figure*}
\centering
\includegraphics[width=.8\linewidth]{./fig3.eps}
\vspace{-6mm}
\caption{Map of the horizontal velocities with magnitude  $> 0.4$ km s$^{-1}$
(71 min average). The contrast within penumbrae has been enhanced (by removal
of a spatial running mean of the original image. Strong radial outflows (moats)
are evident surrounding all filamentary penumbrae. They are confined by the black contours.
These moats are found to be closely associated with the existence of penumbra.
Notice that moat region surrounding coordinates (1,5) corresponds to a
penumbra located to the left of the image and not visible in our FOV. The
coordinates are expressed in Mm. The length of the black bar at coordinates
(0,0) corresponds to $0.4$ km s$^{-1}$
\label{F:3}}
\end{figure*}

\begin{figure}
\hspace{-8mm}
\includegraphics[width=1.2\linewidth]{./fig4.eps}
\vspace{-8mm}
\caption{Corresponding magnetogram of the region under study obtained with the
SOUP filter. The $\delta$-spot neutral lines in Fig.~\ref{F:2} were obtained from this
frame. The coordinates are expressed in Mm.
\label{F:4}}
\end{figure}

\section{Data analysis and results}
The G-band series has been used to study the proper motions of the structures
in the FOV by the local correlation tracking algorithm of \cite{november1988}, as implemented by \cite{molowny1994}. We have chosen a Gaussian
tracking window of FWHM $0\farcs78$ (half of the typical granular size)
suitable for tracking mainly granules. With this procedure, we obtain maps of
horizontal displacements or proper motions per time step (horizontal
velocities), which we average in time. Averages of horizontal velocities have
been performed over 5 min and 71 min intervals.

Figure \ref{F:2} shows the resulting flow map averaged over the whole series (71 min).
Moat flows are seen in the spot on the lower left of the FOV but more prominently
in the spot on the right side (upper and lower penumbral regions).
Several centers of diverging horizontal motions are present in the entire FOV.
Some surround sunspots (lower left and right sides of the figure) but the most
conspicuous ones are seen in the lesser magnetized area
(lower part of the FOV between the spots),
displaying greater velocities and more symmetrical shapes.
These velocity structures, related to recurrent expansion and
splitting of granules, are commonly associated with mesogranulation \citep{roudier2003,roudier2004,bonet2005}.

The white lines in Fig.~\ref{F:2} delineate the line-of-sight (LOS) neutral lines
inferred from a SOUP magnetogram (see Fig.~\ref{F:4}). The regions crossed by this line
correspond naturally to horizontal fields, similar to those found in penumbrae.
However, we caution that these horizontal fields can have different flows as
those that occur in normal penumbra (such as shear flows, as found by, \cite{deng2006} or the supersonic nozzles observed by \cite{pillet1994}.
Close inspection also shows that areas near neutral lines have a filamentary
appearance that is much less distinct than observed in normal penumbrae
extending radially from umbral cores.

The map in Fig.~\ref{F:3} shows only those velocities in the granulation field with
magnitude above $0.4$ km s$^{-1}$. Strong radial outflows (moats) are evident
surrounding the sunspots and $0.4$ km s$^{-1}$ appears to be a characteristic in
the velocity value to define the moats. The black contour in the figure outline
the area of the moats. It has been drawn by hand following visual criteria of
proximity to the sunspot and avoiding a few strong exploding events at the
lower part of the FOV that contaminate the large-scale moat flow.  Empty areas
in the moat region in the velocity maps are a result of the chosen value of the
velocity threshold and correspond to the interaction with exploding events
that occur inside the moat.

From inspection of Fig.~\ref{F:3}, it becomes evident that the moats are closely
associated with the presence of a penumbra. The examples in the lower left, lower
right and upper right corners are quite revealing. The right pointed velocity
vectors on the lower left corner correspond to the moat flow of a penumbra that
is not visible in the FOV of this figure. Moats are missing in the sunspot
sides with no penumbral structure. This is readily seen in the right side of the
umbral core at coordinates (9,10), but also in the left side of the one
at (25,18). None of the pores in the middle of the FOV is associated with any moat-like
flow. More importantly, moats are also absent in granulation regions
that are located side by side with a penumbra, but that are in the
direction perpendicular to that marked by the penumbral filaments.
A clear example is the region around coordinates (9,6-7) in Fig.~\ref{F:3}. As a 'normal'
(meaning here far from the neutral line) penumbra, it is a candidate to develop
a moat flow, but none is seen in this region. The moat flow only appears along the
direction delineated by the penumbral filaments. Moats, thus, appear to
be a natural extension of the flow along the direction of penumbral filaments
and do not exist as outflows in the transverse direction.

The penumbral region in Fig.~\ref{F:3}, centered at coordinates (5,20), does not show clear
evidences of a moat flow (although localized outflows are seen above and below it).
This region corresponds to a strongly sheared neutral line commonly
seen in $\delta$ spots \citep{deng2006}. Figure \ref{F:4} shows a SOUP magnetogram obtained
simultaneously with the G-band series. The spot centered on the left side of the FOV
(white polarity) does not show a radially outward directed penumbra, but penumbral
filaments that run parallel to the neutral line. This is a configuration
commonly referred as sheared. This configuration is also inferred from the
filaments observed in the G-band frames. We suspect that both conditions, the
presence of a neutral line and the absence of radially oriented filaments, are
responsible for the absence of a clearly developed moat flow.
The nature of neutral lines with sheared configurations in $\delta$ spots
is not well understood. But we stress that if they correspond to relatively
shallow structures, they would also block the upward propagation of heat and
should be able to generate a moat flow. Such a moat
flow is not observed in this sunspot.

\section{Conclusions}
Time series of G-band and G-cont images of a sunspot group with
$\delta$-configuration (NOAA active region 10786), spanning over 71 min, were
observed and corrected for atmospheric and instrumental degradation. Proper
motions in the entire FOV have been measured in the G-band series by means of
local correlation tracking techniques. In this paper, we concentrate in the
results obtained in the identified moat flows (or lack thereof) of several
structures within the FOV. The main conclusions from our analysis
can be summarized as follows:

1) We have detected strong (mean speeds of $0.67$ km s${-1}$) outflows streaming
from penumbrae radially oriented from an umbral core, the so-called sunspot moats.
Umbral cores sides with no penumbra {\em do not display moat flows}.

2) Furthermore, the moats are also absent on penumbral sides {\em perpendicular} to
the direction defined by the penumbral filaments. They are not found in
directions transverse to them. A special case is sheared penumbral configurations
tangential to the umbral core. No moat flow is found there either.

These evidence is clearly suggestive of a link between the moat flow
and flows aligned with the penumbral filaments. The possible connection with the
Evershed flow is inescapable. Although it can be argued that a physical
relationship between these two flows is not firmly established in this work,
we believe the evidence is clear-cut and that a statistical analysis of
a larger number of regions, as done here, can establish this connection
more solidly. In particular, work is being done on analysing the moat flow
configuration of round unipolar sunspots lacking penumbra in one of their sectors
to find more evidence supporting points 1) and 2) above.

This result should be put in the context of the recent findings by \cite{sainz2005}; see also \cite{ravindra2006}. These authors find that the
penumbral filaments extend beyond the photometric sunspot boundary and cross
the region dominated by the moat flow. This region is also the location where the MMF
activity is normally found. Indeed some of these MMFs are seen by these authors
to start inside the penumbra. All these results point to an origin linking
the moat flow and MMF activity with the well-known Evershed flow. \cite{cabrera2006} also suggest that the Evershed clouds inside the penumbrae propagate to the moat becoming MMFs once they leave the sunspot.
It is indeed somewhat paradoxical that while the final fate of the Evershed
flow has remained unknown for decades, independent explanations and physical
scenarios have been proposed to generate the moat flow (that starts exactly
where the Evershed flow is seen to vanish).

Similarly, the results obtained from local helioseismology near sunspots, in particular those
related to the f-mode \citep{gizon2000}, are also quite
relevant to the results presented here. The interpretation of these results
in terms of convective cells surrounding sunspots \citep[e.g.,][]{bovelet2003}
may be different if the Evershed flow turns out to be the major process
that injects mass into the moats surrounding sunspots.

We believe the hypothesis put
forward in this paper, the moat flow is the continuation outside the spots
of the Evershed flow, deserves adequate attention both from an observational and a theoretical point of view.

\acknowledgments{
The authors are grateful to Mats L\"ofdahl and Michiel Van Noort, for their
inputs about the restoration process. The Swedish 1-m Solar Telescope is
operated on the island of La Palma by the Institute of Solar Physics of the
Royal Swedish Academy of Sciences in the Spanish Observatorio del Roque de los
Muchachos of the Instituto de Astrof\'isica de Canarias. Partial support by the
Spanish Ministerio de Educaci\'on y Ciencia through project ESP2003-07735-C04 and
financial support by the European Commission throught the SOLAIRE Network (MTRN-CT-2006-035484)
are gratefuly \mbox{acknowledged}.
}

On the Moat-Penumbra Relation


High-frequency diode effect in superconducting Nb$_3$Sn micro-bridges

Abstract: The superconducting diode effect has been recently reported in a variety of systems and different symmetry breaking mechanisms have been examined. However, the frequency range of these potentially important devices still remains obscure. We investigated superconducting micro-bridges of Nb$_{3}$Sn in out-of-plane magnetic fields; optimum magnetic fields of $\sim$10 mT generate $\sim $10% diode efficiency, while higher fields of $\sim$15-20 mT quench the effect. The diode changes its polarity with magnetic field reversal. We documented superconductive diode rectification at frequencies up to 100 kHz, the highest reported as of today. Interestingly, the bridge resistance during diode operation reaches a value that is a factor of two smaller than in its normal state, which is compatible with the vortex-caused mechanism of resistivity. This is confirmed by finite element modeling based on time-dependent Ginzburg-Landau equations. To explain experimental findings, no assumption of lattice thermal inequilibrium was required. Dissimilar edges of the superconductor strip can be responsible for the inversion symmetry breaking by vortex penetration barrier; visual evidence of this opportunity was revealed by scanning electron microscopy. Estimates are in favor of much higher (GHz) range of frequencies for this type of diode.

\section{Introduction}

For certain important problems of fundamental physics, for example,
exploration of quasi-local action of curl-less vector potential \cite{Gulian18}%
, it is necessary to use superconducting micro-bridges for which the
magnitudes of critical currents of superconductor bridge are equal for both
polarities. However, because of simultaneous breaking of time reversal
symmetry (TRS) and inversion symmetry (IS), the phenomenon called the
superconducting diode effect (SDE) arises which violates this
equality \cite{Tokura18,Wakatsuki18,Hoshino18,Ando20,Ideue20,BaumgartnerNN,Wu22,Strambini22}.
The latest research focus in superconductivity has been this effect
with a potential of being applied in significantly
less-dissipative superconducting electronics. In theoretical models, the TRS
is broken by externally applied magnetic field or internal inclusions of
magnetic micro-clusters, while the IS is broken by the out-of-plane Rashba
spin--orbit coupling \cite{Yuan22,Ilic22,Daido22,Karabassov22,He22}, valley-Zeeman interaction \cite{Bauriedl22}, etc.,
which results in the emergence of a chiral superconducting order
\cite{Yuan22,Ilic22,Daido22,He22}. Experimentally, systems based on van der Waals material $\mathrm{%
MoS}_{2}$ with noncentrosymmetric crystal potential \cite{Wakatsuki17}, synthetic
super lattice of $\mathrm{Nb/V/Ta}$ \cite{Ando20}, and some others \cite{Shin21} have
been reported as well as planar Josephson junction arrays of $\mathrm{Al}$
on $\mathrm{InAs}$ \cite{BaumgartnerNN,BaumgartnerIOP}. Yet some other systems reveal nonreciprocal
behavior in field-free environments, such as $\mathrm{NbSe}_{2}$-based
Josephson junction \cite{Wu22} and tri-layer graphene \cite{Lin22}. In view of
variety of experimental observations and theories on SDE in different
configuration of systems, it is quite reasonable to assume that more than
one mechanism can be responsible for the ubiquity of nonreciprocity observations
in superconductor thin films \cite{Hou22}. In a recent article \cite{Suri22}, the
SDE was found in $\mathrm{NbN}$ micro-bridges in an out-of-plane magnetic
field. This observation was attributed to the critical current being
determined by the vortex flow, confirming that the SDE is caused by unequal
vortex barriers on the two edges of the bridge \cite{Hope21,Vodolazov05}.

This mechanism is easy to understand. Consider a type-II superconductor film
strip. A magnetic field above $H_{c1}$ creates vortices; initially, they
nucleate at the strip edges \cite{Shmidt1,Shmidt2}. Morphology of these edges affects
the surface barrier, which prevents the vortices to enter into the strip
\cite{Bean64,FossheimBook}. However, if the current through the strip is
strong enough, the Lorentz force which it exerts onto the vortex overcomes
the surface barrier; vortices start moving across the strip, thus
dissipating energy and creating a resistive state. In practice, vortex
barriers are unequal because of non-identical nature of strip edges.
This circumstance has been used by Vodolazov and Peeters \cite{Vodolazov05}
at predicting superconducting diode effect in 2005 and
experimentally confirmed recently \cite{Suri22}.

The spike of activities in the area of superconducting diodes is paving a
road towards future practical application of these novel devices in
superconducting electronics. However, a very important topic is still open
and remains unaddressed: that is the frequency range of SDE. No reports on
the frequency range above few hertz can be found in the literature which
spreads doubts in practical applicability of the devices. To address this
topic, we explore SDE up to $100$ \textrm{kHz} frequencies, and express
arguments that the frequency range can be even orders of magnitude higher.
This is the main result of our report.

\section{Experimental details}

\begin{figure*}
    \centering
    \includegraphics[width=\linewidth]{Fig1.png}
    \caption{(\textbf{a}) Surface morphology of the films with
inclusion of small areas (nanomountains) with higher concentration of $Sn$
(zoomed tenfold in the inset). (\textbf{b}) Resistive and magnetic
transitions into superconducting state. (\textbf{c}) A bridge after ion
milling of 3D-printed pattern (initial stage of preparation). (\textbf{d})
Last stage of bridge preparation: the structure in panel (\textbf{c}) was
covered by positive photoresist, two circles were projected via
epi-fluorescent microscope and ion milled to reduce the working part of
bridges to micron-scale (2 \textmu m result is shown).
Inset illustrates micro-roughness of the edge line. (e) Effects of ion
milling and post annealing on transition temperature of bridges. (\textbf{f}%
) Typical voltage-current dependence of bridges above ($T=20$ $\mathrm{K}$)
and well-below ($T=1.7$ \textrm{K}) the superconducting transition.}
    \label{fig1}
\end{figure*}

Our study used conventional $\mathrm{Nb}_{3}\mathrm{Sn}$ superconducting
thin film bridges. The $\mathrm{Nb}_{3}\mathrm{Sn}$ films were prepared in a
DC/RF magnetron sputtering system (manufactured by AJA International, Inc.)
with a base pressure of $1 \times 10^{-8}$ \textrm{Torr}. The $\mathrm{Nb}$ target
(Kurt Lesker, $99.95\%$) was placed inside a DC gun while the $\mathrm{Sn}$
target (Kurt Lesker, $99.999\%$) was sputtered using an RF source to avoid
melting. The sapphire substrate (AdValue Technology, thickness $650$ \textmu
\textrm{m}, $C$-cut) was cleaned thoroughly with isopropyl alcohol before it
was mounted on the holder. In our chamber's configuration, the substrate
holder is at the center of the chamber facing upwards, while the sputtering
guns (up to $5$) are located at the top. The substrate is rotated in-plane
throughout the whole deposition process to ensure homogeneous deposition
layer over the whole surface. Our pre-deposition \textit{in-situ} cleaning
of the substrate involves heating it up to $900$ \textrm{C} for $10$ \textrm{%
min} then followed by a gentle bombardment of $\mathrm{Ar}^{+}$ ions at $400$
\textrm{C} for $5$ \textrm{min}. The substrate was oriented to face the ion
gun squarely. For the $\mathrm{Nb}_{3}\mathrm{Sn}$ film deposition, we used
a ``stack+anneal" process similar to the
one in \cite{SundahlThesis}. Taking advantage of $\mathrm{Nb-Sn}$ phase diagram
\cite{Charlesworth70}, an alternating $\mathrm{Nb}$ and $\mathrm{Sn}$ stacks with
an excess $\mathrm{Sn}$ content ($\mathrm{Nb:Sn}$ ratio $\sim 2$) ``phase-locks"
into the $\mathrm{Nb}_{3}\mathrm{Sn}$ composition when the excess $Sn$ evaporates during
post-deposition annealing. For the first/base layer, $\mathrm{Nb}$ was
sputtered in DC mode with $500$ \textrm{W} power and $3$ \textrm{mTorr}
pressure at $600$ \textrm{C} for $20$ \textrm{min}. Then Sn
layer was deposited on top of $\mathrm{Nb}$ in RF mode with $250$ \textrm{W}
power and $30$ \textrm{mTorr} pressure at $100$ \textrm{C} for $15$ \textrm{%
min}. $\mathrm{Nb}$ was deposited one more time as the capping (third) layer
using the same parameters as the base layer for $3$ \textrm{min}. Finally,
the substrate is \textit{in-situ} heated to $950$ \textrm{C} for \textrm{30}
minutes to let the excess $\mathrm{Sn}$ evaporate before cooled down to room
temperature. All the heating/cooling protocols always used a $30$ \textrm{%
C/min} ramp rate. Our films are $\sim 100$ \textrm{nm} thick
and resistivity measurements in PPMS (Quantum Design) confirmed its
stoichiometry and homogeneity from their consistent $T_{c}=17.6$ \textrm{K}
and RRR($300$\textrm{K}/$20$\textrm{K}) $\sim 4$ between
samples across several batches; the stoichiometry was also checked by the
SEM EDX (Hitachi SU3500, Oxford Instruments X-MAX-20).

Figure 1 demonstrates
certain physical properties of our films (panels \textbf{a} and \textbf{b})
and bridges (panels \textbf{c}-\textbf{f}).

Patterning of the films into bridges has been performed using combination of
$3D-$printing (Elegoo Mars-3 printer), photolithography and by ion milling with
macroscopic contact pads on $1\times 1$ \textrm{cm}$^{\mathrm{2}}$ sapphire
substrates. Reactive ion milling (Bal-Tec RES-101, $\mathrm{CF}_{4}$
etchant) was used to develop metallic film pattern, Fig.~\ref{fig1}(\textbf{c}).
After removal of the resin, the structure was covered by negative
photoresist, and the 100 \textmu m-scale bridge was narrowed down
to 12 \textmu m with the projective photolithography (using mask
projection via LUMAM epi-fluorescent microscope, $10\times $ objective).
After ion milling, the last stage of patterning was undertaken using
positive photoresist and the same projective technique with the $40\times $
objective and a different mask in the form of two holes. The final ion
milling delivered bridges about 10 \textmu m long and down
to 2 \textmu m wide, Fig.~\ref{fig1}(\textbf{d}).

\begin{figure*}[t]
    \centering
    \includegraphics[width=0.8\linewidth]{Fig2.png}
    \caption{(\textbf{a}) Resistive state at various values of bias
temperature for the applied current of opposite polarities. (\textbf{b}) SDE
with sinusoidal current amplitude $\approx 10$ \textrm{mA} and frequency $0.1
$ \textrm{Hz} in 2 \textmu m bridge. (\textbf{c}) Same as
in panel (\textbf{b}) with 5 \textmu m-wide bridge at frequency $10$ \textrm{kHz}. (\textbf{d}) SDE at
frequency $100$ \textrm{kHz} (average of $200$ acquisitions) at reversed
polarity of magnetic field. In all cases, a magnetic field orthogonal to the surface of the bridge was applied and optimized in the range 50-100 Oe.}
    \label{fig2}
\end{figure*}

Ion milling affects the physical properties of bridges thus reducing the
critical temperatures down to the $2-12$ K range. This $T_{c}$%
-reduction was also noticed in the literature \cite{Suri22}. We were able to restore
the $T_{c} \sim 17$ \textrm{K} values in some of the bridges by
high-vacuum post-annealing ($\sim 10^{-6}$ \textrm{Torr}) at
$900$ \textrm{C} during $30$ \textrm{min}, as shown in Fig.~\ref{fig1}(\textbf{e}).

\section{Results}

The tests demonstrated $V(I_{+})\neq V(I_{-})$, Fig.~1(f), thus indicating
SDE. To study the SDE, temperature dependence of $2$ $\mathrm{\mu m}$
bridge was measured, Fig.~\ref{fig2}(\textbf{a}).

This measurement confirms the fact mentioned in literature that at lower
temperatures the difference between the threshold values of $I_{+}^{res}$\
and $I_{-}^{res}$\ increases with the decreasing bias temperature \cite{Ando20,Hou22} (the
value of $I^{res}$ corresponds to the current at which the resistive state
sets up). The bridge then was biased at $2$ \textrm{K} temperature (PPMS
DynaCool cryostat) and the expected diode effect was confirmed, Fig.~\ref{fig2}(%
\textbf{b}), with externally applied current source and nanovoltmeter
(Keithley 6221 and 2182a respectively) at frequency $0.1$ \textrm{Hz}.

Since superconducting diodes are considered as important elements for
electronics, it was meaningful to register the effect at possibly higher
frequencies. The Keithley current source generates AC-currents up to $100$
\textrm{kHz}. Detection of voltage output was performed by oscilloscope
(Tektronix TDS 644A). The results are shown at Fig.~\ref{fig2}(\textbf{c}) and (\textbf{d})
for $10$ \textrm{kHz} and $100$ \textrm{kHz} frequencies, respectively. Higher frequency
measurements require specially designed circuitry in the cryostat since the
noise becomes an issue, as is visible from the comparison of Figs. 2(\textbf{%
d}) and (\textbf{b}). For reducing noise, we plotted $100$ \textrm{kHz}
data averaged over $200$ acquisitions.

\section{Discussion}

\begin{figure}
    \centering
    \includegraphics[width=\linewidth]{Fig3.png}
    \caption{$V(I)$ dependence of 2 \textmu m bridge. Resistive state jumps are noticeable
before the transition to the normal resistive state.}
    \label{fig3}
\end{figure}

To estimate the frequency range of this type of diodes, one needs to
understand what is the underlying mechanism of the SDE. Comparison of panels
(\textbf{a}) and (\textbf{b}) in Fig.~2 shows the diode polarity change
because of switching the external magnetic field from $100$ \textrm{Oe} to $%
-100$ \textrm{Oe}. The field is orthogonal to the bridge surface. The field
amplitude is optimal for SDE in our bridges. Higher values ($\sim 150$
\textrm{Oe}) quench the effect. The resistive state which yields the
non-zero voltage in Fig.~2(\textbf{a}) is caused by the vortex motion: the corresponding
value of resistance is about the half of its normal value. Indeed, as
follows from Fig.~3, the resistance ($V/I$) reaches the value of $25-40$
\textrm{Ohm}, while in normal state, just above the superconducting
transition, the resistance is $\sim 60$ \textrm{Ohm}. Figure 3 depicts a set
of $V(I)$ measurements at various temperatures. Multiple resistivity jumps
are noticeable at lower temperature curves. This kind of the $V(I)$ curves
behavior is in close correspondence with the results for Abrikosov and
kinematic vortices and phase-slip lines revealed by modeling in thin-film
bridges \cite{Suri22,Berdiyorov09,Lotero14,Vodolazov07}. As was shown
in via modeling in these works, different
values of the applied current correspond to different vortex lattice
patterns. We reproduced the computations by Vodolazov and Peeters \cite{Vodolazov07} using
time-dependent Ginzburg-Landau (TDGL) equations and COMSOL Multiphysics
package in a slightly different model: our bridge has geometrically broken
IS: its top edge has three shallow dents (with depth less than the vortex
lateral size), while the bottom edge is smooth. Such a model allows us not
only explain the $V(I)$ characteristics in Fig.~3 similarly to
\cite{Vodolazov07}, but also certain diode features.

We will accentuate on the latter ones here. Figure 4 demonstrates the vortex
lattice pattern in the given external magnetic field $H<H_{c2}$ at both
polarities of applied current.

\begin{figure}
    \centering
    \includegraphics[width=\linewidth]{Fig4.png}
    \caption{Vortex dynamics in the IS and TRS broken bridge modeled via
TDGL-COMSOL approach. Panels demonstrate vortices in the central, narrow part
of the bridge. Current flows horizontally from the left to the right in
panels (\textbf{a}) and (\textbf{c}), and in the reverse direction in two
other cases. Color coding corresponds to the Cooper pair density: from the
blue for normal metal state to the brown in the unperturbed pair-density
value. Magnetic field and transport current are switched on simultaneously.
Panels (\textbf{a}) and (\textbf{b})  correspond to the same moments of
elapsed time. Panels (c) and (d) correspond to twice higher values of the
transport current and four times larger elapsed times.}
    \label{fig4}
\end{figure}

As follows from this figure, the evolution is essentially different due to
the presence of (three) small dents in the top edge of the bridges. It
breaks the IC symmetry, and the evolution of the system differs for opposite
polarities of the current. At large enough values of the current, the
individual vortices are consolidating into lines crossing the width of the
bridge. They are classified as Phase Slip Lines (PSL). Again, the presence
of minor dents creates asymmetry between the opposite current polarities: to
generate the full width crossing by PSL the positive polarity current
required $50\%$ longer time than the negative current.

In case of our experiments, the IS breaking, which is the second required
condition for SDE, is being provided, most likely, by the breaking of
left-right edge symmetry in lithographically processed bridges. Closer look
to Fig.~1 (\textbf{d}, Inset) illustrates the fractured structure of the
edge of the current flow channel.

Associating the SDE with the vortex motion, we can make certain estimates
for the\ performance limits of SDE in the frequency domain. The $100$
\textrm{kHz} results have been obtained in 5 \textmu m-wide bridge.
One can relate this frequency with the reciprocal characteristic vortex
crossing time of the bridge. If the geometry is scaled down to $5-50$
\textrm{nm}-wide bridges, it can deliver $2-3$ orders of magnitude higher
frequencies because of geometric factors. Can even higher frequency range be
reached? The answer looks positive because of the following reasons: first,
the $100$ \textrm{kHz} may not be the true limit for the $5$ \textmu \textrm{m}
bridge. Second, the speed of kinematic vortices is $3$ orders of magnitude
higher than that of the Abrikosov vortices [33]. If the $I_{+/-}^{res}$
asymmetry spreads into the kinematic vortex state (second $V(I)\ $jump area
in Fig.~3), then the operational frequencies of superconducting diodes can
reach $\mathrm{GHz}$ and even higher levels. Obviously, more research is
required for reaching that stage.

\section{Summary}

Ordinary thin superconducting bridges (in our case, Nb$_{3}$Sn) can
demonstrate SDE because of simultaneous breakage of IS and TRS. Natural
reasons of breaking the IS is the geometric asymmetry between the bridge
edges, and the TRS can be broken by the remnant magnetic field in the
cryostat or by the current via the bridge itself. Finite element modeling of
physical phenomena in the bridge based on the TDGL equations reveal complex
microscopic nature of the resistive states in it. More experiments guided by
the results of this modeling can yield higher-performance two terminal and
multi-terminal devices useful for practice.

\begin{acknowledgments}

The work of Chapman U. research team is supported by the ONR grants N00014-21-1-2879 and
N00014-20-1-2442.

\end{acknowledgments}

\providecommand{\noopsort}[1]{}\providecommand{\singleletter}[1]{#1}%