
This paper introduces the trajectory divergence rate, a scalar field which
locally gives the instantaneous attraction or repulsion of adjacent
trajectories. This scalar field may be used to find highly attracting or
repelling invariant manifolds, such as slow manifolds, to rapidly approximating
hyperbolic Lagrangian coherent structures, or to provide the local stability of
invariant manifolds. This work presents the derivation of the trajectory
divergence rate and the related trajectory divergence ratio for 2dimensional
systems, investigates their properties, shows their application to several
example systems, and presents their extension to higher dimensions.

The PHENIX experiment at the Relativistic Heavy Ion Collider has measured the
differential cross section of $\phi$(1020)meson production at forward rapidity
in $p$$+$$p$ collisions at $\sqrt{s}=$510 GeV via the dimuon decay channel. The
partial cross section in the rapidity and $p_T$ ranges $1.2<y<2.2$ and
$2<p_T<7$ GeV/$c$ is $\sigma_\phi=[2.28 \pm 0.09\,{\rm (stat)} \pm 0.14\,{\rm
(syst)} \pm 0.27\, {\rm (norm)}] \times 10^{2}$~mb. The energy dependence of
$\sigma_\phi$ ($1.2<y<2.2, \; 2<p_T<5$ GeV/$c$) is studied using the PHENIX
measurements at $\sqrt{s}=$200 and 510 GeV and the LargeHadronCollider
measurements at $\sqrt{s}=$2.76 and 7 TeV. The experimental results are
compared to various event generator predictions ({\sc pythia6, pythia8, phojet,
ampt, epos3,} and {\sc eposlhc}).

We present a detailed measurement of charged twopion correlation functions
in 0%30% centrality $\sqrt{s_{_{NN}}}=200$ GeV Au$+$Au collisions by the
PHENIX experiment at the Relativistic Heavy Ion Collider. The data are well
described by BoseEinstein correlation functions stemming from L\'evystable
source distributions. Using a fine transverse momentum binning, we extract the
correlation strength parameter $\lambda$, the L\'evy index of stability
$\alpha$ and the L\'evy length scale parameter $R$ as a function of average
transverse mass of the pair $m_T$. We find that the positively and the
negatively charged pion pairs yield consistent results, and their correlation
functions are represented, within uncertainties, by the same L\'evystable
source functions. The $\lambda(m_T)$ measurements indicate a decrease of the
strength of the correlations at low $m_T$. The L\'evy length scale parameter
$R(m_T)$ decreases with increasing $m_T$, following a hydrodynamically
predicted type of scaling behavior. The values of the L\'evy index of stability
$\alpha$ are found to be significantly lower than the Gaussian case of
$\alpha=2$, but also significantly larger than the conjectured value that may
characterize the critical point of a secondorder quarkhadron phase
transition.

Production of $\pi^0$ and $\eta$ mesons has been measured at midrapidity in
Cu$+$Au collisions at $\sqrt{s_{_{NN}}}$=200 GeV. Measurements were performed
in $\pi^0(\eta)\rightarrow\gamma\gamma$ decay channel in the 1(2)20 GeV/$c$
transverse momentum range. A strong suppression is observed for $\pi^0$ and
$\eta$ meson production at high transverse momentum in central Cu$+$Au
collisions relative to the $p$$+$$p$ results scaled by the number of
nucleonnucleon collisions. In central collisions the suppression is similar to
Au$+$Au with comparable nuclear overlap. The $\eta/\pi^0$ ratio measured as a
function of transverse momentum is consistent with $m_T$scaling
parameterization down to $p_T=$2 GeV/$c$, its asymptotic value is constant and
consistent with Au$+$Au and $p$$+$$p$ and does not show any significant
dependence on collision centrality. Similar results were obtained in
hadronhadron, hadronnucleus, and nucleusnucleus collisions as well as in
$e^+e^$ collisions in a range of collision energies $\sqrt{s_{_{NN}}}=$31800
GeV. This suggests that the quarkgluonplasma medium produced in Cu$+$Cu
collisions either does not affect the jet fragmentation into light mesons or it
affects the $\pi^0$ and $\eta$ the same way.

We have measured direct photons for $p_T<5~$GeV/$c$ in minimum bias and
0\%40\% most central events at midrapidity for Cu$+$Cu collisions at
$\sqrt{s_{_{NN}}}=200$ GeV. The $e^{+}e^{}$ contribution from quasireal
direct virtual photons has been determined as an excess over the known hadronic
contributions in the $e^{+}e^{}$ mass distribution. A clear enhancement of
photons over the binary scaled $p$$+$$p$ fit is observed for $p_T<4$ GeV/$c$ in
Cu$+$Cu data. The $p_T$ spectra are consistent with the Au$+$Au data covering a
similar number of participants. The inverse slopes of the exponential fits to
the excess after subtraction of the $p$$+$$p$ baseline are
285$\pm$53(stat)$\pm$57(syst)~MeV/$c$ and 333$\pm$72(stat)$\pm$45(syst)~MeV/$c$
for minimum bias and 0\%40\% most central events, respectively. The rapidity
density, $dN/dy$, of photons demonstrates the same power law as a function of
$dN_{\rm ch}/d\eta$ observed in Au$+$Au at the same collision energy.

PHENIX has measured the azimuthal correlations of muon pairs from charm and
bottom semileptonic decays in $p$+$p$ collisions at $\sqrt{s}=200$ GeV, using
a novel analysis technique utilizing both unlike and likesign muon pairs to
separate charm, bottom and DrellYan contributions. The dimuon measurements
combined with the previous electronmuon and dielectron measurements span a
wide range in rapidity, and are well described by PYTHIA Tune A. Through a
Bayesian analysis based on PYTHIA Tune A, we show that leading order pair
creation is the dominant $(76\%\pm^{14}_{19}\%)$ contribution for $b\bar{b}$
production, whereas the data favor the scenario in which nexttoleadingorder
processes dominate $c\bar{c}$ production. The small contribution of
nexttoleadingorder processes in $b\bar{b}$ production at the collision
energies of the Relativistic Heavy Ion Collider contrasts with the case at
LargeHadronCollider energies, where nexttoleadingorder processes are
expected to dominate.

The PHENIX collaboration presents first measurements of lowmomentum
($0.4<p_T<3$ GeV/$c$) directphoton yields from Au$+$Au collisions at
$\sqrt{s_{_{NN}}}$=39 and 62.4 GeV. For both beam energies the directphoton
yields are substantially enhanced with respect to expectations from prompt
processes, similar to the yields observed in Au$+$Au collisions at
$\sqrt{s_{_{NN}}}$=200. Analyzing the photon yield as a function of the
experimental observable $dN_{\rm ch}/d\eta$ reveals that the lowmomentum
($>$1\,GeV/$c$) directphoton yield $dN_{\gamma}^{\rm dir}/d\eta$ is a smooth
function of $dN_{\rm ch}/d\eta$ and can be well described as proportional to
$(dN_{\rm ch}/d\eta)^\alpha$ with $\alpha{\sim}$1.25. This new scaling behavior
holds for a wide range of beam energies at the Relativistic Heavy Ion Collider
and Large Hadron Collider, for centrality selected samples, as well as for
different, $A$$+$$A$ collision systems. At a given beam energy the scaling also
holds for high $p_T$ ($>5$\,GeV/$c$) but when results from different collision
energies are compared, an additional $\sqrt{s_{_{NN}}}$dependent
multiplicative factor is needed to describe the integrateddirectphoton yield.

The experimental study of the collisions of heavy nuclei at relativistic
energies has established the properties of the quarkgluon plasma (QGP), a
state of hot, dense nuclear matter in which quarks and gluons are not bound
into hadrons. In this state, matter behaves as a nearly inviscid fluid that
efficiently translates initial spatial anisotropies into correlated momentum
anisotropies among the produced particles, producing a common velocity field
pattern known as collective flow. In recent years, comparable momentum
anisotropies have been measured in smallsystem protonproton ($p$$+$$p$) and
protonnucleus ($p$$+$$A$) collisions, despite expectations that the volume and
lifetime of the medium produced would be too small to form a QGP. Here, we
report on the observation of elliptic and triangular flow patterns of charged
particles produced in protongold ($p$$+$Au), deuterongold ($d$$+$Au), and
heliumgold ($^3$He$+$Au) collisions at a nucleonnucleon centerofmass energy
$\sqrt{s_{_{NN}}}$~=~200 GeV. The unique combination of three distinct initial
geometries and two flow patterns provides unprecedented model discrimination.
Hydrodynamical models, which include the formation of a shortlived QGP
droplet, provide a simultaneous description of these measurements.

PHENIX reports differential cross sections of $\mu\mu$ pairs from
semileptonic heavyflavor decays and the DrellYan production mechanism
measured in $p$$+$$p$ collisions at $\sqrt{s}=200$ GeV at forward and backward
rapidity ($1.2<\eta<2.2$). The $\mu\mu$ pairs from $c\bar{c}$, $b\bar{b}$,
and DrellYan are separated using a simultaneous fit to unlike and likesign
muon pair spectra in mass and $p_T$. The azimuthal opening angle correlation
between the muons from $c\bar{c}$ and $b\bar{b}$ decays and the pair$p_T$
distributions are compared to distributions generated using {\sc pythia} and
{\sc powheg} models, which both include nexttoleading order processes. The
measured distributions for pairs from $c\bar{c}$ are consistent with {\sc
pythia} calculations. The $c\bar{c}$ data presents narrower azimuthal
correlations and softer $p_T$ distributions compared to distributions generated
from {\sc powheg}. The $b\bar{b}$ data are well described by both models. The
extrapolated total cross section for bottom production is $3.75{\pm}0.24({\rm
stat}){\pm}^{0.35}_{0.50}({\rm syst}){\pm}0.45({\rm global})$[$\mu$b], which is
consistent with previous measurements at the Relativistic Heavy Ion Collider in
the same system at the same collision energy, and is approximately a factor of
two higher than the central value calculated with theoretical models. The
measured DrellYan cross section is in good agreement with
nexttoleadingorder quantumchromodynamics calculations.

Dihadron and isolated direct photonhadron angular correlations are measured
in $p$$+$$p$ collisions at $\sqrt{s}=200$ GeV. The correlations are sensitive
to nonperturbative initialstate and finalstate transverse momentum $k_T$ and
$j_T$ in the azimuthal nearly backtoback region $\Delta\phi\sim\pi$. In this
region, transversemomentumdependent evolution can be studied when several
different hard scales are measured. To have sensitivity to small transverse
momentum scales, nonperturbative momentum widths of $p_{\rm out}$, the
outofplane transverse momentum component perpendicular to the trigger
particle, are measured. These widths are used to investigate possible effects
from transversemomentumdependent factorization breaking. When accounting for
the longitudinal momentum fraction of the awayside hadron with respect to the
nearside trigger particle, the widths are found to increase with the hard
scale; this is qualitatively similar to the observed behavior in DrellYan and
semiinclusive deepinelastic scattering interactions. The momentum widths are
also studied as a function of centerofmass energy by comparing to previous
measurements at $\sqrt{s}=510$ GeV. The nonperturbative jet widths also appear
to increase with $\sqrt{s}$ at a similar $x_T$, which is qualitatively
consistent to similar measurements in DrellYan interactions. To quantify the
magnitude of any transversemomentumdependent factorization breaking effects,
calculations will need to be performed to compare to these measurements.

We report the transverse singlespin asymmetries of $J/\psi$ production at
forward and backward rapidity, $1.2<y<2.2$, as a function of $J/\psi$
transverse momentum ($p_T$) and Feynman$x$ ($x_F$). The data analyzed were
recorded by the PHENIX experiment at the Relativistic Heavy Ion Collider in
2015 from $p$$+$$p$, $p$$+$Al, and $p$$+$Au collisions with transversely
polarized proton beams at $\sqrt{s_{_{NN}}}=200$ GeV. At this collision energy,
singlespin asymmetries for heavyflavor particle production of $p$$+$$p$
collisions provide access to the spindependent gluon distribution and
highertwist correlation functions inside the nucleon, such as the gluon
QiuSterman and trigluon correlation functions. Proton+nucleus collisions offer
an excellent opportunity to study nuclear effects on the correlation functions.
The data indicate negative asymmetries at the twostandarddeviation level in
the $p$$+$Au data for $p_T<2$ GeV/$c$ at both forward and backward rapidity,
while in $p$$+$$p$ and $p$$+$Al collisions the asymmetries are consistent with
zero within the range of experimental uncertainties.

We present measurements of elliptic and triangular azimuthal anisotropy of
charged particles detected at forward rapidity $1<\eta<3$ in Au$+$Au
collisions at $\sqrt{s_{_{NN}}}=200$ GeV, as a function of centrality. The
multiparticle cumulant technique is used to obtain the elliptic flow
coefficients $v_2\{2\}$, $v_2\{4\}$, $v_2\{6\}$, and $v_2\{8\}$, and triangular
flow coefficients $v_3\{2\}$ and $v_3\{4\}$. Using the smallvariance limit, we
estimate the mean and variance of the eventbyevent $v_2$ distribution from
$v_2\{2\}$ and $v_2\{4\}$. In a complementary analysis, we also use a folding
procedure to study the distributions of $v_2$ and $v_3$ directly, extracting
both the mean and variance. Implications for initial geometrical fluctuations
and their translation into the final state momentum distributions are
discussed.

Hindman's Theorem (HT) states that for every coloring of $\mathbb N$ with
finitely many colors, there is an infinite set $H \subseteq \mathbb N$ such
that all nonempty sums of distinct elements of $H$ have the same color. The
investigation of restricted versions of HT from the computabilitytheoretic and
reversemathematical perspectives has been a productive line of research
recently. In particular, HT$^{\leqslant n}_k$ is the restriction of HT to sums
of at most $n$ many elements, with at most $k$ colors allowed, and HT$^{=n}_k$
is the restriction of HT to sums of \emph{exactly} $n$ many elements and $k$
colors. Even HT$^{\leqslant 2}_2$ appears to be a strong principle, and may
even imply HT itself over RCA$_0$. In contrast, HT$^{=2}_2$ is known to be
strictly weaker than HT over RCA$_0$, since HT$^{=2}_2$ follows immediately
from Ramsey's Theorem for $2$colorings of pairs. In fact, it was open for
several years whether HT$^{=2}_2$ is computably true.
We show that HT$^{=2}_2$ and similar results with addition replaced by
subtraction and other operations are not provable in RCA$_0$, or even WKL$_0$.
In fact, we show that there is a computable instance of HT$^{=2}_2$ such that
all solutions can compute a function that is diagonally noncomputable relative
to $\emptyset'$. It follows that there is a computable instance of HT$^{=2}_2$
with no $\Sigma^0_2$ solution, which is the best possible result with respect
to the arithmetical hierarchy. Furthermore, a careful analysis of the proof of
the result above about solutions DNC relative to $\emptyset'$ shows that
HT$^{=2}_2$ implies RRT$^{=2}_2$, the Rainbow Ramsey Theorem for $2$colorings
of pairs, over RCA$_0$. The most interesting aspect of our construction of
computable colorings as above is the use of an effective version of the
Lov\'asz Local Lemma due to Rumyantsev and Shen.

For commutative rings, we introduce the notion of a {\em universal grading},
which can be viewed as the "largest possible grading". While not every
commutative ring (or order) has a universal grading, we prove that every {\em
reduced order} has a universal grading, and this grading is by a {\em finite}
group. Examples of graded orders are provided by group rings of finite abelian
groups over rings of integers in number fields. We generalize known properties
of nilpotents, idempotents, and roots of unity in such group rings to the case
of graded orders; this has applications to cryptography. Lattices play an
important role in this paper; a novel aspect is that our proofs use that the
additive group of any reduced order can in a natural way be equipped with a
lattice structure.

We have measured the cross section and single spin asymmetries from forward
$W^{\pm}\rightarrow\mu^{\pm}\nu$ production in longitudinally polarized
$p$$+$$p$ collisions at $\sqrt{s}=510$ GeV using the PHENIX detector at the
Relativistic Heavy Ion Collider. The cross sections are consistent with
previous measurements at this collision energy, while the most forward and
backward longitudinal single spin asymmetries provide new insights into the sea
quark helicities in the proton. The charge of the W bosons provides a natural
flavor separation of the participating partons.

It is now well established that a Josephson junction made from conventional
spinsinglet superconductors containing ferromagnetic layers can carry
spintriplet supercurrent under certain conditions. The first experimental
signature of that fact is the propagation of such supercurrent over long
distances through strong ferromagnetic materials. Surprisingly, one of the most
salient predictions of the theory has yet to be verified experimentally 
namely that a Josephson junction containing three magnetic layers with coplanar
magnetizations should exhibit a groundstate phase shift of either zero or pi
depending on the relative orientations of those magnetizations. Here we
demonstrate this property using Josephson junctions containing three different
types of magnetic layers, chosen so that the magnetization of one layer can be
switched by 180 degrees without disturbing the other two. Phasesensitive
detection is accomplished using a superconducting quantum interference device,
or SQUID. Such a phasecontrollable junction could be used as the memory
element in a fullysuperconducting computer.

It has been proposed recently that a previously unobserved neutron decay
branch to a dark matter particle ($\chi$) could account for the discrepancy in
the neutron lifetime observed in experiments that use two different measurement
techniques. One of the possible final states discussed includes a single $\chi$
along with an $e^{+}e^{}$ pair. We use data from the UCNA (Ultracold Neutron
Asymmetry) experiment to set limits on this decay channel. Coincident
electronlike events are detected with $\sim 4\pi$ acceptance using a pair of
detectors that observe a volume of stored Ultracold Neutrons (UCNs). The summed
kinetic energy ($E_{e^{+}e^{}}$) from such events is used to set limits, as a
function of the $\chi$ mass, on the branching fraction for this decay channel.
For $\chi$ masses consistent with resolving the neutron lifetime discrepancy,
we exclude this as the dominant dark matter decay channel at $\gg~5\sigma$
level for $100~\text{keV} < E_{e^{+}e^{}} < 644~\text{keV}$. If the
$\chi+e^{+}e^{}$ final state is not the only one, we set limits on its
branching fraction of $< 10^{4}$ for the above $E_{e^{+}e^{}}$ range at $>
90\%$ confidence level.

We present measurements of azimuthal correlations of charged hadron pairs in
$\sqrt{s_{_{NN}}}=200$ GeV Au$+$Au collisions after subtracting an underlying
event using a model that includes higherorder azimuthal anisotropy $v_2$,
$v_3$, and $v_4$. After subtraction, the awayside ($\Delta\phi\sim\pi)$ of the
highest transversemomentum trigger ($p_T>4$ GeV/$c$) correlations is
suppressed compared to that of correlations measured in $p$$+$$p$ collisions.
At the lowest associated particle $p_T$, the awayside shape and yield are
modified. These observations are consistent with the scenario of radiativejet
energy loss. For the lowest$p_T$ trigger correlations, an awayside yield
exists and we explore the dependence of the shape of the awayside within the
context of an underlyingevent model. Correlations are also studied
differentially versus eventplane angle $\Psi_n$. The angular correlations show
an asymmetry when selecting the sign of the triggerparticle azimuthal angle
with respect to the $\Psi_2$ event plane. This asymmetry and the measured
suppression of the pair yield out of plane is consistent with a
pathlengthdependent energy loss. No $\Psi_3$ dependence can be resolved
within experimental uncertainties.

In this paper, we investigate an electrodynamics in which the physical modes
are coupled to a Lorentzviolating (LV) background by means of a
higherderivative term. We analyze the modes associated with the dispersion
relations (DRs) obtained from the poles of the propagator. More specifically,
we study Maxwell's electrodynamics modified by a LV operator of mass dimension
6, that is, by a term that possesses two additional derivatives coupled to a
CPTeven tensor $D_{\beta\alpha}$ that plays the role of the fixed background.
The new term possesses the % form $\partial_{\sigma}F^{\sigma\beta}
\partial_{\lambda} F^{\lambda\alpha}D_{\beta\alpha}$. We first evaluate the %
propagator and obtain the dispersion relations of the theory. By doing so, we
analyze some configurations of the fixed background % and search for sectors
where the energy is welldefined and causality is assured. In the sequel, we
perform the same kind of analysis for a more general dimension6 model. We
conclude that the modes of both Lagrange densities are plagued by physical
problems, including causality violation, and that signal propagation becomes
physically meaningful only in the highmomentum regime.

We study how doublestranded DNA translocates through graphene nanogaps.
Nanogaps are fabricated with a novel capillaryforce induced graphene nanogap
formation technique. DNA translocation signatures for nanogaps are
qualitatively different from those obtained with circular nanopores, owing to
the distinct shape of the gaps discussed here. Translocation time and
conductance values vary by $\sim 100$%, which we suggest are caused by local
gap width variations. We also observe exponentially relaxing current traces. We
suggest that slow relaxation of the graphene membrane following DNA
translocation may be responsible. We conclude that DNAgraphene interactions
are important, and need to be considered for graphenenanogap based devices.
This work further opens up new avenues for direct read of single molecule
activitities, and possibly sequencing.

A Gadget, more precisely a scalar Gadget, is defined as a mathematical
calculation acting over a domain of one or more adinkra graphs and whose range
is a real number. A 2010 work on the subject of automorphisms of adinkra
graphs, implied the existence of multiple numbers of Gadgets depending on the
number of colors under consideration. For four colors, this number is two. In
this work, we verify the existence of a second such Gadget and calculate (both
analytically and via explicit computerenabled algorithms) its 1,358,954,496
matrix elements over 36,864 minimal valise adinkras related to the Coxeter
Group BC4.

As a straightforward generalization and extension of our previous paper, J.
Phys. A50 (2017) 215201 we study aspects of the quantum and classical dynamics
of a $3$body system with equal masses, each body with $d$ degrees of freedom,
with interaction depending only on mutual (relative) distances. The study is
restricted to solutions in the space of relative motion which are functions of
mutual (relative) distances only. It is shown that the ground state (and some
other states) in the quantum case and the planar trajectories (which are in the
interaction plane) in the classical case are of this type. It corresponds to a
threedimensional quantum particle moving in a curved space with special
$d$dimensionindependent metric in a certain $d$dependent singular potential,
while at $d=1$ it elegantly degenerates to a twodimensional particle moving in
flat space. It admits a description in terms of pure geometrical
characteristics of the interaction triangle which is defined by the three
relative distances. The kinetic energy of the system is $d$independent, it has
a hidden $sl(4,R)$ Lie (Poisson) algebra structure, alternatively, the hidden
algebra $h^{(3)}$ typical for the $H_3$ Calogero model as in the $d=3$ case. We
find an exactlysolvable threebody $S^3$permutationally invariant,
generalized harmonic oscillatortype potential as well as a
quasiexactlysolvable threebody sextic polynomial type potential with
singular terms. For both models an extra first order integral exists. It is
shown that a straightforward generalization of the 3body (rational) Calogero
model to $d>1$ leads to two primitive quasiexactlysolvable problems. The
extension to the case of nonequal masses is straightforward and is briefly
discussed.

Distribution of quantum correlations among remote users is a key procedure
underlying many quantum information technologies. EinsteinPodolskyRosen
steering, which is one kind of such correlations stronger than entanglement,
has been identified as a resource for secure quantum networks. We show that
this resource can be established between two and even more distant parties by
transmission of a system being separable from all the parties. For the case
with two parties, we design a protocol allowing to distribute oneway Gaussian
steering between them via a separable carrier; the obtained steering can be
used subsequently for onesided deviceindependent (1sDI) quantum key
distribution. Further, we extend the protocol to three parties, a scenario
which exhibits richer steerability properties including onetomultimode
steering and collective steering, and which can be used for 1sDI quantum secret
sharing. All the proposed steering distribution protocols can be implemented
with squeezed states, beam splitters and displacements, and thus they can be
readily realized experimentally. Our findings reveal that not only entanglement
but even steering can be distributed via communication of a separable system.
Viewed from a different perspective, the present protocols also demonstrate
that one can switch multipartite states between different steerability classes
by operations on parts of the states.

Recently, all 1,358,954,496 values of the Gadget between the 36,864 adinkras
with four colors, four bosons, and four fermions have been computed in
1701.00304. In this paper, we further analyze these results in terms of $BC_3$,
the signed permutation group of three elements, and $BC_4$, the signed
permutation group of four elements. It is shown how all 36864 adinkras can be
generated via $BC_4$ boson $\times$ $BC_3$ color transformations of two
quaternion adinkras that satisfy the quaternion algebra. An adinkra inner
product has been used for some time, known as the gadget, that is used to
distinguish adinkras. We show how 96 equivalence classes of adinkras that are
based on the gadget emerge in terms of $BC_3$ and $BC_4$. We also comment on
the importance of the gadget as it relates to separating out dynamics in terms
of K\"ahlerlike potentials.

Josephson junctions containing two ferromagnetic layers are being considered
for use in cryogenic memory. Our group recently demonstrated that the
groundstate phase difference across such a junction with carefully chosen
layer thicknesses could be controllably toggled between zero and $\pi$ by
switching the relative magnetization directions of the two layers between the
antiparallel and parallel configurations. However, several technological issues
must be addressed before those junctions can be used in a largescale memory.
Many of these issues can be more easily studied in single junctions, rather
than in the Superconducting QUantum Interference Device (SQUID) used for the
phasesensitive measurements. In this work, we report a comprehensive study of
spinvalve junctions containing a Ni layer with a fixed thickness of 2.0 nm,
and a NiFe layer of thickness varying between 1.1 and 1.8 nm in steps of 0.1
nm. We extract the field shift of the Fraunhofer patterns and the critical
currents of the junctions in the parallel and antiparallel magnetic states, as
well as the switching fields of both magnetic layers. We also report a partial
study of similar junctions containing a slightly thinner Ni layer of 1.6 nm and
the same range of NiFe thicknesses. These results represent the first step
toward mapping out a ``phase diagram" for phasecontrollable spinvalve
Josephson junctions as a function of the two magnetic layer thicknesses.