Experiments and simulations demonstrating the rapid ultrasonic rewarming of frozen beef cryovials

Abstract: The development of methods to safely rewarm large volume cryopreserved biological samples remains a barrier to the widespread adoption of cryopreservation. Here, experiments and simulations were performed to demonstrate that ultrasound can increase rewarming rates relative to thermal conduction alone. An ultrasonic rewarming setup based on a custom 444 kHz tubular piezoelectric transducer was designed, characterized, and tested with 2 mL cryovials filled with frozen ground beef. Rewarming rates were characterized in the -20$^{\circ}$C to 5$^{\circ}$C range. Thermal conduction-based rewarming was compared to thermal conduction plus ultrasonic rewarming, demonstrating a ten-fold increase in rewarming rate when ultrasound was applied. The maximum recorded rewarming rate with ultrasound was 57$^{\circ}$C per minute, approximately 2.5 times faster than with thermal conduction alone. Coupled acoustic and thermal simulations were developed and showed good agreement with the heating rates demonstrated experimentally and were also used to demonstrate spatial heating distributions with small ($<$3$^{\circ}$C) temperature differentials throughout the sample when the sample was below 0$^{\circ}$C. The experiments and simulations performed in this work demonstrate the potential for ultrasound as a rewarming method for cryopreserved tissues, as faster rewarming rates may improve the viability of cryopreserved tissues and reduce the time needed for cells to regain normal function.

\section{Introduction}

There is an unmet need for effective cryopreservation methods that can be implemented in cellular therapies and for transplantable tissue sections\cite{giwa2017promise,urbani2018multi}, in part due to challenges in rewarming these materials without causing damage.
Cryopreservation of biological materials may be implemented whether via vitrification or by slow freezing.
Vitrification is the cooling of a medium into a disordered and non-crystalline state in the added presence of cryoprotective agents.
Vitrification reduces cellular damage by reducing the likelihood of ice crystal nucleation and growth\cite{fahy1987biological,fahy2015principles}.
Ice crystal nucleation is more probable at low sub-zero temperatures while ice crystal growth is more probable at higher sub-zero temperatures.
This presents a problem for rewarming a vitrified medium, as any nucleated ice formed during cooling will grow as the medium is rewarmed\cite{fahy1987biological}.
In consequence, rewarming rates that preserve vitrified cell viability generally exceed cooling rates by at least an order of magnitude\cite{fahy2015principles}.
An alternative to vitrification is slow freezing\cite{kilbride2014scale,urbani2017long} in which the medium is frozen slowly at a controlled and optimized rate and in the presence of ice nucleating and cryoprotectant agents\cite{massie2014gmp,kilbride2014scale}.
The medium can then be revived with rewarming rates as low as 1$^{\circ}$C per minute\cite{kilbride2017cryopreservation} depending on the rate of freezing.
However, faster rewarming may improve important outcomes such as cell viability, the restoration of cellular lactate production, and the reduction of bulk thermomechanical stresses if the volume is rewarmed homogeneously \cite{kilbride2017cryopreservation}.
Furthermore, faster rewarming may reduce the time needed post-rewarming to regain normal cellular function \cite{kilbride2017cryopreservation}.

The most common rewarming method is immersion in a warm water bath, which can effectively generate fast rewarming rates (up to 200$^{\circ}$C/min for volumes of the order of 100\,$\mu$L\cite{pegg1984effect}.
However, large ($\geq$3\,mL) volume rewarming cannot be performed effectively with a warm water bath due to the inherently low thermal conduction within biological media, along with the large thermal mass and thermal lag of a large cryopreserved volume \cite{taylor2019new,kilbride2016spatial}.
Furthermore, thermal gradients generate mechanical stresses in the medium that can cause gross damage to a cryopreserved organ or tissue\cite{fahy1990physical,rubinsky1980thermal}.
Alternative rewarming methods that don't rely on thermal conduction have been investigated.
Electromagnetic (dielectric) methods can generate fast rewarming of the order of 100$^{\circ}$\,C/min in larger sub-radiofrequency wavelength spherical and ellipsoidal volumes\cite{robinson2002electromagnetic}.
However, there is a positive relationship between most dielectric constants and temperature, making it difficult to control the heat distribution resulting in runaway heating effects\cite{lu2000combined,wusteman2004vitrification}.
An alternative electromagnetic heating method uses alternating magnetic fields to excite a ferrofluid\cite{rosensweig2002heating}.
Magnetic nanoparticles can be perfused into a volume without significantly changing the medium cooling rate\cite{etheridge2014rf}.
The medium can then be heated at a rate at least equal to dielectric heating, but with less potential for thermal runaway or the generation of large thermal gradients\cite{etheridge2014rf,manuchehrabadi2017improved,wang2015theoretical,solanki2017thermo}.
This method has been applied to vitrified rat and rabbit kidneys, and rat hearts to generate mean heating rates above 50$^{\circ}$C/min without significant thermomechanical damage or ice crystal formation \cite{sharma2021vitrification,gao2022vitrification}.
However, scaling inductive heating to larger human-scaled volumes will result in an increase in heat deposition away from the center of the heating coil due to the formation of eddy currents \cite{sharma2021vitrification}.
This may present a challenge for obtaining a uniform heating distribution.
An additional challenge with inductive heating is obtaining a homogeneous distribution of magnetic nanoparticles, particularly in organs that are less homogeneously vascularized than a kidney.
Another new method for rewarming cryopreserved materials is laser gold nanowarming\cite{khosla2018characterization}, which has been used to rewarm vitrified millimetre-sized droplets containing biological media and can maintain high cellular viability\cite{khosla2018characterization,zhan2021cryopreservation,zhan2022pancreatic}.
However, scaling this method past 2D scaffolds\cite{zhan2021cryopreservation} to 3D volumes remains a challenge and there remains a need for a rewarming method with a good penetration depth which minimises hot spots.

Here, ultrasound is explored as a potential rewarming technique.
Ultrasound has long been used to heat normothermic tissues, with many hyperthermia and high intensity focused ultrasound applications either in clinical trials or in use as approved medical procedures\cite{kennedy2003high,izadifar2020introduction}.
Ultrasound in the sub-MHz to low-MHz range has a high penetration depth in soft tissues at normal physiological temperatures, making it suited to heating the center of large volumes.
Ultrasound can be focused electronically using transducer arrays or geometrically\cite{o1949theory} by shaping the source or using lenses or acoustic holograms\cite{lalonde1993field}.
Electronic focusing and steering of array transducers can additionally be employed to shape or move the focal region over a large volume\cite{do1981annular,hindley2004mri,kohler2009volumetric}.
Ultrasound may be better suited to rewarming than electromagnetic methods because ultrasound attenuation (which is proportional to the volume rate of heat deposition) decreases as a biological tissue rewarms, reducing the likelihood of thermal runaway\cite{shore1986attenuation}.
Ultrasonic attenuation results from scattering and absorption, the latter of which generates heating\cite{cobbold2006foundations}.
Heating rates in normothermic biological tissues are well characterized and implemented in high intensity focused ultrasound treatments for many conditions including uterine fibroids, pancreatic cancer, and prostate cancer \cite{izadifar2020introduction}, although the contributions of absorption and scattering to attenuation in these tissues is not fully characterised.
It is not known for frozen tissues what proportion of ultrasound is absorbed rather than scattered and this will require investigation.
An experimental demonstration of the rewarming of a frozen tissue sample is needed to demonstrate the viability of ultrasound as a method to rapidly rewarming cryopreserved tissues.

Sound has previously been used to thaw frozen meat and fish, although these experiments were often done in water baths and at sub-100\,kHz frequencies\cite{rosenberg1974method,kissam1982water,guo2021ultrasound}.
In this early work, most of the energy supplied to thaw the materials was supplied by the water bath, with the ultrasound waves assisting in the heat transfer into the materials, without significant energy deposition from the ultrasound itself.
Similar ultrasound-assisted systems are used as freezers, and it has been shown that intermittent ultrasound can increase freezing rates\cite{li2002effect} and increase ice nucleation during freezing\cite{zhang2018using}.
It has also been shown that 0.4 - 1\,MHz ultrasound can directly thaw frozen meat and fish, using an experimental design where the frozen medium was insulated and not subject to thermal conduction from the surrounding medium\cite{miles1999high}.
However, the thawing times from -10$^{\circ}$C were between tens of minutes and hours depending on the distance from a piston source, and with a close-to linear increasing relationship between thawing times and the distance from the source \cite{miles1999high}.
A common trait of these prior works is that the ultrasonic sources were unfocused, resulting in increased energy deposition in the near field, and focused sources were not investigated.
A simulation study used multiple sources focused to an ovary to simulate fast (100$^{\circ}$\,C/min) rewarming from low sub-zero temperatures to high sub-zero temperatures\cite{olmo2020use}.
However, simulations did not model rewarming above -20$^{\circ}$C, where the acoustic and thermal properties of biological media have a strong temperature dependence, and experiments are required to confirm the simulated rewarming rates.

In this work, focused ultrasound at 444\,kHz was used to rapidly rewarm biological samples from -20$^{\circ}$C to above 5$^{\circ}$C, and ultrasonic rewarming rates were measured for comparison against other rewarming methods.
The ultrasonic rewarming setup was developed specifically to rewarm biological samples in 2\,mL cryovials.
Cryovials are commonly used for the cryopreservation of samples (e.g., sperm\cite{li2014ivf}, embryonic stem cells\cite{nagy2006freezing}, and mouse ovaries\cite{candy1997effect}).
The cryovial rewarming experiments were supported by additional experiments that reduce the rewarming effect of thermal conduction from the rewarming setup, demonstrating rewarming primarily caused by ultrasound absorption.

\section{Methods}

\subsection{Design Criteria}

A custom ultrasonic rewarming setup was designed and built.
The design process for the ultrasonic rewarming setup began with three criteria that were intended to ensure that the rewarming setup achieved the rewarming rates necessary for rewarming a slow-cooled material.
The rewarming setup should:\\
1. Repeatably deliver ultrasound to a frozen, filled cryovial without excessive reflective losses.\\
2. Generate a radially centered heat distribution within the cryovial, with limited thermal conduction-based warming.\\
3. Maximize the rewarming rate while maintaining a rewarming volume of sufficient size.\\
The first criterion requires that the ultrasound source be acoustically coupled to the cryovial via a medium that minimizes the impedance mismatch with the frozen biological material stored in the cryovial.
Olmo $et$ $al.$ used ice as the coupling medium in their simulation study\cite{olmo2020use}, but ice may cause expansion damage during freezing and may decouple during rewarming.
Natural nylon 6.6 (Direct Plastics UK Ltd.) was chosen instead as the coupling medium, as it has a similar acoustic impedance to frozen beef (See Table \ref{tbl:MediumProps}), reducing reflective losses, and a relatively low acoustic attenuation compared to other plastics.

A prospective simulation study was performed to obtain a source geometry satisfying the second criterion, and a source frequency satisfying the third criterion.
Three source geometries were investigated: planar sources (cryovial placed at the location of the last axial maximum), focused arc sources (cryovial placed at the geometric focus), and tubular sources (cryovial placed at the center).
For simplicity, only single-source setups were considered.
Figure \ref{fig:DesignSims} shows examples of the pressure magnitude fields of the three types of sources (equal source powers, 444\,kHz) embedded in nylon 6.6, continuously sonicating a cryovial containing frozen beef.

\begin{figure}
    \centering
    \includegraphics[width = \columnwidth]{Figure1.eps}
    \caption{Pressure fields in 12.5\,mm diameter cryovials placed at the a) last axial maximum of a 444\,kHz planar source, b) geometric focus of a 444\,kHz focused source, and c) center of a 444\,kHz tubular source. Pressure fields normalized by c).}
    \label{fig:DesignSims}
\end{figure}

Figure \ref{fig:DesignSims} shows that the tubular source generates the highest pressure, and the pressure distribution is radially centered.
The tubular source is advantaged by its matched geometry to the cryovial.
The frequency displayed in Fig. \ref{fig:DesignSims} shows a central focus and an outer ring encapsulated within the cryovial.
Halving the frequency would result in a central focus covering the entire cryovial, but would likely more-than halve the rewarming rate due to the frequency-dependence of ultrasonic attenuation and corresponding heat deposition.
In this work the decision was made to use a higher frequency source (400+ kHz as shown in Fig. \ref{fig:DesignSims}) to maximize the rewarming rate at the center of the cryovial.

\subsection{Device Construction}

A tubular modified lead zirconate-lead titanate transducer element was purchased from PI Ceramic (Lederhose, Germany).
The transducer element was made from a proprietary ceramic (PIC181), suited to high power acoustic applications and applications in resonance mode \cite{piceramic}.
The element was coated with a fired silver electrode (thick film, typically 10\,$\mu$m) covering both inner and outer surfaces, and was radially polarized with the +pole on the outer radial surface.
The outer diameter of the element was 76\,mm, the inner diameter 66\,mm, and the height was 50\,mm.
The radial resonant frequency of the transducer element was first approximated by:
\begin{equation}\label{eq:freq}
    f = \dfrac{N_t}{TH},
\end{equation}
where $TH$ is the thickness of the transducer element (5\,mm) and $N_t$ is the acousto-mechanical frequency coefficient in Hz$\cdot$m (2110\,Hz$\cdot$m for PIC181\cite{piceramic}).
Equation \ref{eq:freq} estimates a resonant frequency of 422\,kHz, but the accuracy of Eq. \ref{eq:freq} decreases when the tube height is less than the tube outer diameter, as is the case for this transducer.
A BNC cable was soldered to the element electrodes and the actual resonant frequency was obtained by maximizing the transmitted:reflected power ratio, giving a resonant frequency of 444\,kHz.
All experiments were thus performed at 444\,kHz, without electrical matching.

\begin{figure}
    \centering
    \includegraphics[width = \columnwidth]{Figure2.eps}
    \caption{The rewarming setup. a) Transducer dimensions. b) Nylon 6.6 mount dimensions. c) The mounted transducer and a cryovial half-inserted into the nylon 6.6 mount.}
    \label{fig:Tx_Nylon}
\end{figure}

A 71\,mm diameter natural nylon 6.6 rod was machined to fit inside the tubular transducer with a 0.25\,mm tolerance (See Fig. \ref{fig:Tx_Nylon}).
The 0.25\,mm radial tolerance allows for a sliding fit and avoids damaging the transducer via thermal expansion of the nylon 6.6 rod (coefficient of linear thermal expansion: 80 $\mu$m/(m$\cdot$K)) and thermal contraction of the transducer (coefficient of linear thermal expansion -4 to -6 $\mu$m/(m$\cdot$K) in the radial direction \cite{piceramic}).
The 0.25\,mm gap was filled with deionized water to couple the transducer to the nylon 6.6 rod.
The gap allowed for thermal expansion of the nylon 6.6 and transducer thermal contraction (up to 80$^{\circ}$C), or swelling of the nylon 6.6 due to water absorption.
A 12.5\,mm diameter hole was machined in the center of the rod for the placement of a cryovial.
Air-backing the tubular transducer further increases the heating efficiency by directing the acoustic energy toward the cryovial and forming a resonant cavity.

\subsection{Biological Medium}

Ground beef was chosen as the biological medium to be rewarmed, as the thermal and acoustic properties have been well characterized in work published in the food sciences literature for temperatures ranging from -20$^{\circ}$C to 20$^{\circ}$C \cite{barrera1983thermal,miles1974changes,shore1986attenuation,tavman2007apparent}.
The temperature-driven changes in beef sound speed\cite{miles1974changes}, attenuation\cite{shore1986attenuation}, thermal conductivity\cite{barrera1983thermal}, and specific heat\cite{tavman2007apparent} are displayed in Fig. \ref{fig:BeefAcousticThermalProps}.
Beef has also been used in prior work in ultrasonic thawing, allowing for comparisons of rewarming rates\cite{miles1999high}.
Ground beef was used to increase the ease of handling and shaping.

\begin{figure}
    \centering
    \includegraphics[width = \columnwidth]{Figure3.eps}
    \caption{Acoustic and thermal properties of beef as a function of temperature: a) sound speed\cite{miles1974changes}, b) attenuation\cite{shore1986attenuation} fit with splines,  c) thermal conductivity \cite{barrera1983thermal}, and d) specific heat\cite{tavman2007apparent}.}
    \label{fig:BeefAcousticThermalProps}
\end{figure}

The changes in beef acoustic and thermal properties over this temperature range present challenges and opportunities for ultrasound-based rewarming.
Thermal conductance decreases threefold as the beef thaws, reducing the diffusion of hot spots.
The sound speed decreases by a factor of nearly two, changing the acoustic impedance mismatch between the beef and the coupling medium.
The attenuation decreases several-fold, reducing the volume rate of heat deposition in thawed samples (this is advantageous for avoiding thermal runaway, but still reduces mean heating rate of a sample).
The specific heat curve shown in Fig. \ref{fig:BeefAcousticThermalProps} incorporates the latent heat of fusion and increases by over an order of magnitude as the beef temperature increases towards zero, reducing the heating rate but potentially increasing the homogeneity of heating within the frozen beef.

\subsection{Coupled acoustic-thermal simulations}

Acoustic and thermal simulations were performed to help design and predict the rewarming behaviour of the experimental setup.
Continuous wave ultrasound propagation was simulated using the k-Wave fluid 2D model (\texttt{kspaceFirstOrder2D}) from version 1.3 of the open-source k-Wave toolbox \cite{treeby2010k,treeby2012modeling} and using Matlab 2019a.
Rewarming was simulated using the k-Wave \texttt{kWaveDiffusion} code, which solves the Pennes Bioheat Equation\cite{pennes1948analysis,treeby2015contribution}.
Thermal conduction was included in the thermal model, allowing heat to flow from the nylon 6.6 to the beef.
The simulation domain was created using the \texttt{makeDisc} function to generate masks of the nylon 6.6 holder and the frozen beef, and the \texttt{makeCircle} function to generate the tubular source.
The initial temperature distribution was set as 20$^{\circ}$C in the surrounding air and nylon 6.6 and -20$^{\circ}$C in the frozen beef.
The output metrics of the simulation were the spatial compressional pressure profiles in the simulation domain in Pa and the corresponding volume rates of heat deposition $Q$ in W/m$^3$.
$Q$ was calculated using the plane wave assumption, i.e. $Q = \alpha p^2/(\rho c)$ , where $\alpha$ is compressional attenuation in Np/m, $p$ is pressure in Pa, $\rho$ is density in kg/m$^3$, and $c$ is compressional sound speed in m/s.

The acoustic and thermal parameters of frozen beef, nylon 6.6, and air used in simulation are listed in Table \ref{tbl:MediumProps}.
The acoustic and thermal properties of beef displayed in Table \ref{tbl:MediumProps} are those at -20$^{\circ}$C.
The compressional sound speed and attenuation in the nylon 6.6 was measured experimentally at 20$^{\circ}$C.
The acoustic and thermal properties of room temperature air were obtained from the Engineering Toolbox website\cite{ETB}.
The density of air was multiplied by 100 to reduce the computation grid requirements while maintaining effectively the same interfacial behaviour as expected with the physical air density value.

\begin{table}
\caption{\label{tbl:MediumProps}The acoustic and thermal properties (density $\rho$, sound speed $c$, attenuation coefficient $\alpha$, specific heat $C$, and thermal conductivity $k$) of beef at -20$^{\circ}$C and nylon 6.6 and air at room temperature (500\,kHz).}
\begin{ruledtabular}
\begin{tabular}{ccccc}
Property                        &  Beef                      & Nylon 6.6               & Air \\
\hline
$\rho$ [kg/m$^3$]               & 1073 \cite{barrera1983thermal}   & 1120 \cite{Onda}  & 116* \\
$c$ [m/s]                       & 2996 \cite{miles1974changes}     & 2586        & 343 \cite{ETB}\\
$\alpha$  [dB/(cm$\cdot$MHz)]         & 1.21 \cite{shore1986attenuation} & 0.9       & 0\\
$C$ [J/(kg$\cdot$K)]                  & 902 \cite{barrera1983thermal}    & 1700 \cite{ETB}    & 1000 \cite{ETB}   \\
$k$ [W/(m$\cdot$K)]                   & 1.4908 \cite{tavman2007apparent} & 0.25 \cite{ETB}   & 0.03  \cite{ETB}  \\
\end{tabular}
\end{ruledtabular}
\end{table}

The accurate representation of changes in acoustic and thermal properties of the beef with rewarming requires coupling of acoustic and thermal distributions in simulation.
This was approximated with consecutive acoustic $\rightarrow$ thermal $\rightarrow$ acoustic ... simulations with an update frequency `f' of the acoustic properties and fields and thermal properties.
The change in beef density was small relative to the changes in other acoustic and thermal properties and was consequently fixed at the values reported in Table \ref{tbl:MediumProps}.
The model geometry includes the air surrounding the transducer and coupling medium but convection and air flow around the device are not modelled.
This model should be an acceptable approximation of the rewarming experiments, as the heat transfer coefficients between the nylon 6.6 holder and the surrounding air is low.

The coupled acoustic and thermal simulations were performed in 2D, as a horizontal cross-section of the experimental setup.
The grid size for the 2D simulations was 432 x 432 (83.4324\,mm x 83.4324\,mm, dx =  0.193\,mm), and the acoustic simulations were run for 344\,$\mu$s of acoustic wave propagation (5.6\,ns time step, CFL = 0.01, 20 grid point perfectly matched layer).
The computation time for one acoustic simulation was approximately 30\,s using an NVIDIA GeForce Titan X GPU in an Intel(R) Xeon(R) CPU E5-2650 0 @ 2.00GHZ server and the \texttt{kspaceFirstOrder2DG} k-Wave code.
The grid size for the 2D rewarming steps was the same (432 x 432), the temporal discretization was 1\,ms, and simulating 200 time steps took approximately 2\,s using the \texttt{kWaveDiffusion} code on the server.
The total computation time, including pre- and post-processing steps, for 100\,s of simulated ultrasonic rewarming was 7.7 hours.
In the thermal simulations, the initial temperature distribution was generated by setting the ground beef temperature to -20$^{\circ}$C and the surrounding nylon 6.6 and air temperature to 20$^{\circ}$C.
This ensured that thermal conduction-based rewarming was incorporated into the model, mimicking the experimental rewarming conditions.
During the first 8 seconds of simulated rewarming, $Q$ was set to zero, representing the time taken to transfer the cryovials from the freezer to the rewarming setup and increase the acoustic power of the transducer to the desired value.

A set of coupled ultrasound and thermal simulations were performed, where the acoustic and thermal properties and fields were updated sequentially.
The update frequency of acoustic and thermal properties and acoustic field was 5\,Hz, obtained via a convergence study that found a $<$1.5$^{\circ}$C/s maximum absolute difference relative to a further doubling of the update frequency.
Figure \ref{fig:updateconvergence} shows the effect of the update rate on temperature evolution: fixing the acoustic and thermal properties at the values shown in Table \ref{tbl:MediumProps} for 10\,s (f = 0.1\,Hz) of heating resulted in an average rewarming rate nearly 400$^{\circ}$C per minute, an over-prediction of realistic rewarming rates.
Further increases in the field update frequency may be necessary to accurately simulate faster rewarming and regions with larger thermal gradients.

\begin{figure}
    \centering
    \includegraphics[width = \columnwidth]{Figure4.eps}
    \caption{The simulated heating rate decreases as the simulation update frequency (f) increases. Convergence to under a degree difference per time point occurs at f = 5\,Hz.}
    \label{fig:updateconvergence}
\end{figure}

\subsection{3D axisymmetric simulations}

The 3D axisymmetric code (\texttt{kspaceFirstOrderAS}) was used to estimate the spatial pressure magnitude distributions\cite{treeby2020nonlinear}.
The simulation domain for the 3D axisymmetric simulation was created using a 2D plane and 2D masks (See Fig. \ref{fig:AxisymmetricFields}) and includes the beef, radially surrounded by the nylon 6.6 holder, which was vertically capped and backed by air.
The source was created using the \texttt{makeLine} function, which becomes a tube when the 2D plane is transformed to axisymmetric coordinates.
Due to the computational cost, the 3D code was not used to perform consecutive acoustic and thermal rewarming simulations.
Running the axisymmetric code with the same CFL number and spatial discretization as the 2D simulations (input grid 324 x 243, 62.5743\,mm by 46.9307\,mm) takes approximately 50\,minutes.
With a similar increase in the computation time for the 3D thermal simulation,  full 3D acoustic and thermal simulations of 100\,s of experimental rewarming would take close to 800 hours.

\subsection{Experiments}

The acoustic and rewarming capabilities of the tubular transducer were evaluated experimentally.
The acoustic fields generated by the tubular transducer were characterized using hydrophone line and plane measurements, while the rewarming rates were characterized with thermocouples embedded in beef samples.
All analysis was performed in Matlab 2019a.

\subsubsection{Acoustic Characterization Experiments}

A set of ultrasound field characterization measurements were performed for the tubular transducer and nylon 6.6 holder.
The tubular transducer was coupled to the nylon holder with deionized water and sealed at the base of the connection between the transducer and nylon 6.6 holder.
The 12.5\,mm nylon 6.6 tube hole was also filled with deionized water for the field characterization measurements.
The ultrasound field was measured inside the tube using a tapered fibre-optic hydrophone (Precision Acoustics, Dorchester UK) mounted to the 3-axis scanning arm of a UMS automated scanning tank (Precision Acoustics, Dorchester UK).
The tapered fibre-optic hydrophone is small (diameter of approximately 120$\mu$m and active element diameter of 10$\mu$m) and the sensitivity is nearly invariant with respect to incident wave angle for frequencies below 1\,MHz, so a correction for the 90$^{\circ}$ incidence angle was not made\cite{morris2005development}.
The transducer was driven with a pulsed (2 cycle) or quasi-continuous signal (100 cycles) generated by a Keysight 33500B Series Waveform Generator and amplified by an  E\&I 1020L, 200\,W amplifier to 50\,V (pulse repetition frequency: 250\,Hz).
The fibre-optic hydrophone voltages (32 averages, 50\,Ms/s, 5000 samples) were digitized with a Tektronix DPO5034B oscilloscope and recorded.
The UMS system was used to perform 1D scans (axial scan dimensions: 40\,mm range, 0.2\,mm steps) and 2D scans (axial dimensions: 20\,mm, 0.5\,mm steps and lateral dimensions: 6\,mm range, 0.2\,mm steps) in the nylon 6.6 tube.
The 1D and 2D scan positions are shown relative to the nylon 6.6 tube in Fig. \ref{fig:EXPfieldprofile}d).
1D scans were repeated across the maximum pressure location with increasing driving voltages (25\,V to 135\,V in approximately 18\,V increments) to generate a pressure vs. voltage calibration.
The fibre-optic hydrophone was calibrated via comparison with a calibrated 0.2\,mm needle hydrophone using measurements of the field of a 270\,kHz source and a 400\,kHz source, and linearly extrapolating to 444\,kHz.
Voltage waveforms were then converted to pressure using this calibration.
Fibre-optic hydrophone measurements were not feasible at higher driving voltages and pressures; at higher pressures the optical fibre flexibility allowed the whole fibre to vibrate with the pulse-repetition frequency of the transducer, distorting the recorded waveforms by changing the position of the active element of the fibre optic hydrophone.

\subsubsection{Rewarming Experiments}

Rewarming experiments were performed with ground beef in 2\,mL cryovials.
Type T thermocouples were used to record the temperature during rewarming experiments (5SRTC-TT-TI-40-1M, Omega, Norwalk, CT, USA).
The thermocouples had a conductor diameter of 0.08\,mm, a bead diameter of approximately 0.24\,mm, and a cable outer diameter of 0.4 x 0.7\,mm.
The thermocouple cable insulation was stripped to expose the conductors, reducing the thermocouple wire diameter exposed to the ultrasound field.
The thermocouple wire diameter was approximately 1/40th of the source wavelength in frozen beef and was sufficiently small relative to the source wavelength to avoid heating artifacts by minimizing the scattering cross-section of the thermocouple and by minimizing the heat capacity and consequent thermal lag of the thermocouple \cite{fry1954determination}
The thermocouple may generate a viscous heating artifact when the ground beef is in a thawed state \cite{fry1954determination,morris2008investigation}.
However, thermocouple-based viscous heating is unlikely to contribute to rewarming in a frozen biological medium.
Ice is much more viscous than soft tissue (at least nine orders of magnitude)\cite{fowler1997glaciers} and calculating the viscous heating artifact for an ultrasound field intercepting a perpendicular thermocouple embedded in ice gives a negligible rate of heat deposition\cite{fry1954determination,tiennot2019numerical}, however, the viscous heating artefact has not been widely explored in solid materials.

Corning 2\,mL non-self standing cryovials were used for the rewarming experiments.
Thermocouples were strung between holes punched at the centers of the bottoms of the cryovials and the cryovial caps.
Thermocouple measurements were recorded at 1\,Hz using an NI-TC01 (National Instruments) temperature logger, with the recording starting prior to cryovial loading into the rewarming setup.
Ground British beef (5\% fat content) was loaded into the cryovials and around the thermocouple wires.
The thermocouple beads were pulled back 1\,cm into the cryovial, embedding the bead in the beef (Fig. \ref{fig:ThermocouplePlacement}).
The filled cryovials and embedded thermocouples were then loaded into a -20$^{\circ}$C freezer (LEC-medical) and left to freeze overnight.

\begin{figure}
    \centering
    \includegraphics[width = \columnwidth]{Figure5.eps}
    \caption{Thermocouple placement a) in the cryovials and b) in the ground beef-filled transducer (r denotes the distance from the center of the transducer).}
    \label{fig:ThermocouplePlacement}
\end{figure}

The cryovials were loaded individually into the rewarming setup.
The transducer was then driven with a 220\,V amplitude continuous wave signal, corresponding to a 160\,W output from the amplifier.
The signal amplitude at the signal generator was adjusted incrementally to maintain the amplifier power output as the electrical impedance of the transducer changed due to transducer warming.
The nylon 6.6 holder was cooled in room-temperature water for at least 10 minutes before each experiment to minimize differences in rewarming due to thermal conduction.

Two additional rewarming experiments were performed to demonstrate ultrasonic rewarming with negligible thermal-conduction based heating.
One experiment was performed with the transducer fully filled with ground beef.
Four thermocouples were embedded, one at the center of the transducer, and the other three spaced 1\,cm apart, extending radially from the first thermocouple (Fig. \ref{fig:ThermocouplePlacement}b).
The four thermocouples were connected to a TC-08 datalogger (Pico Technologies, Corby, UK), connected to a workstation for temperature logging.
The transducer ground-beef thermocouple assembly was frozen overnight in the -20$^{\circ}$C freezer.
The assembly was then removed from the freezer and rewarmed using the same heating protocol that was used for the cryovial experiments.
A second experiment was performed with a large volume nylon 6.6 holder, identical to the holder used for the cryovials but with a hole with an internal diameter of 46\,mm and a depth of 50\,mm in order to isolate the transducer from the ground beef and replicate the interfaces of the cryovial experiments.
The large volume holder was filled with ground beef (total volume of approximately 80\,mm$^3$), and two thermocouples were embedded (one at the center and one placed 6\,mm radially).
The nylon - ground beef assembly was also frozen overnight in the -20$^{\circ}$C freezer, and in this case the transducer was left at room temperature.

\section{Results}

\subsection{Simulations of the acoustic field in the nylon tube holder}

Simulations of the experimental setup were performed with the k-Wave axisymmetric code to generate 2D slices of the 3D pressure distribution.
The continuous-wave fields are shown for frozen beef (Fig.\ref{fig:AxisymmetricFields}a) and water (Fig.\ref{fig:AxisymmetricFields}b) in the nylon 6.6 tube.
The simulations show that the source generates a centered focus with a radial distribution of lobes, with some vertical variation in the pressure distribution in the beef but much more in the water distribution.
The water pressure field varies axially due to the low attenuation of ultrasound in water, giving the ultrasound waves longer to interact with the tube boundary conditions and develop axial standing wave patterns.
The combination of standing waves in the water pressure distribution and the higher attenuation of frozen beef results in higher pressures in the water simulation; the maximum pressure in water in Fig. \ref{fig:AxisymmetricFields} is 2.1 times the maximum pressure in the equivalent frozen beef simulation.

\begin{figure}
    \centering
    \includegraphics[width = \columnwidth]{Figure6.eps}
    \caption{Continuous-wave fields in the nylon 6.6 tube holder and a) the tube filled with frozen beef, or b) the tube filled with water. Magnitudes normalized by the water maximum.}
    \label{fig:AxisymmetricFields}
\end{figure}

\subsection{Experiments}

\subsubsection{Transducer Characterization}

Pressure fields were generated by the transducer with either short bursts or pseudo-continuous wave (CW) signals.
Figure \ref{fig:EXPfieldprofile}a) shows 1D measurements of the burst and CW conditions, taken along the line depicted in Fig. \ref{fig:EXPfieldprofile}d).
Figure \ref{fig:EXPfieldprofile}b) and c) show 2D measurements of the burst and CW conditions, respectively.
The 2D measurement planes are shown in Fig. \ref{fig:EXPfieldprofile}d).

\begin{figure}
    \centering
    \includegraphics[width = \columnwidth]{Figure7.eps}
    \caption{a) Measured axial field profile at the geometric center of the tube, with the transducer driven with a short burst or with a quasi-continuous wave (CW). Measured 2D field profiles for a b) short burst signal and c) quasi-CW signal. d) 1D and 2D scan positions in the water-filled tube.}
    \label{fig:EXPfieldprofile}
\end{figure}

\begin{figure}[ht]
    \centering
    \includegraphics[width = \columnwidth]{Figure8.eps}
    \caption{Amplified driving voltage versus pressure in water, extrapolated using simulation to frozen beef.}
    \label{fig:Voltage_vs_Pressure}
\end{figure}

Figure \ref{fig:Voltage_vs_Pressure} shows the relationship between the amplified driving voltage and the maximum recorded pressure in the nylon 6.6 tube.
Figure \ref{fig:Voltage_vs_Pressure} extends the linear fits to higher driving voltages corresponding to those used in the rewarming experiments; the fibre-optic hydrophone was unable to make accurate pressure measurements at these higher driving voltages because the hydrophone began resonating in the acoustic field.
The simulated water:frozen beef pressure magnitude ratio (2.1) was used to estimate the maximum pressure magnitude in frozen beef for a given driving voltage.
The estimated pressure amplitude at the driving voltage used for the rewarming experiments ($\approx$220\,V) was approximately 0.8\,MPa.
Assuming a linear signal, the mechanical index $in$ $situ$ is given by $MI = 0.8 / \sqrt{0.444} = 1.2$, below the threshold of 1.9 necessary to avoid cavitation in the presence of microbubbles.
This $in$ $situ$ mechanical index indicates that the sonications should not result in mechanical damage to the biological medium encased in the cryovial.

\subsubsection{Rewarming Experiments}

The cryovial temperatures for the four thermal conduction-only rewarming trials and the four ultrasound + thermal conduction rewarming trials are shown in Fig. \ref{fig:RewarmingEXPs}a), along with the mean temperatures for each experiment type.
The rewarming rates were obtained by taking temperature time derivatives of the mean temperatures and are displayed in Fig. \ref{fig:RewarmingEXPs}b).
Rewarming with ultrasound was at least twice as fast and on average over 10 times faster than with thermal conduction alone, despite the relatively fast rewarming induced by the large temperature gradient between the cryovial and nylon 6.6 holder.
The thawing time (defined as the time to rewarm from -2$^{\circ}$C to 2$^{\circ}$C) was over an order of magnitude faster when ultrasound was applied.
Rewarming rates are reported for the thermal conduction-only and thermal conduction + ultrasound experiments in Table \ref{tbl:RewarmingRates}.

\begin{figure}
    \centering
    \includegraphics[width = \columnwidth]{Figure9.eps}
    \caption{a) Rewarming temperature rise versus time with thermal conduction-only or ultrasound + thermal conduction. b) Heating rates versus time for the two types of experiments. }
    \label{fig:RewarmingEXPs}
\end{figure}

\begin{table}
\caption{\label{tbl:RewarmingRates}Rates of temperature rise between -20$^{\circ}$C and 5$^{\circ}$C for cryovial rewarming experiments.}
\begin{ruledtabular}
\begin{tabular}{ccc}
                    & Thermal             & Thermal Cond. \\
                    & Conduction          & + Ultrasound \\
\hline
Max. Heating Rate   & 23$^{\circ}$C/min   & 57$^{\circ}$C/min           \\
Mean Heating Rate   & 2$^{\circ}$C/min    & 23$^{\circ}$C/min           \\
Min. Heating Rate   & 0$^{\circ}$C/min    & 5$^{\circ}$C/min           \\
Thaw time (-2°C to 2°C) & 290$\pm$30\,s     & 24$\pm$2\,s               \\
\end{tabular}
\end{ruledtabular}
\end{table}

An additional rewarming experiment was performed with the tubular transducer filled with ground beef and thermocouples (Fig. \ref{fig:ThermocouplePlacement}).
The four thermocouple temperature recordings for this experiment are shown in Fig. \ref{fig:BurgerRewarming}a).
The rewarming rate was highest at the central thermocouple, which averages a rewarming rate of 50$^{\circ}$C per minute from 4\,s to 15\,s.
The approximately 4\,s delay corresponds to the time it took to turn on the source and increase the driving voltage to 220\,V.
The average rewarming rate at the central thermocouple was faster than the average rewarming rate seen in the cryovial experiments for the same initial rewarming period.
The faster rewarming rate may be due to fewer interfaces reflecting sound away from the center of the large volume beef sample.
This experiment demonstrates ultrasound-based rewarming with minimal influence from thermal conduction and demonstrates the absorption of ultrasonic energy into the frozen beef.

\begin{figure}
    \centering
    \includegraphics[width = \columnwidth]{Figure10.eps}
    \caption{The rewarming curves for a) four thermocouples embedded in the transducer-frozen beef assembly, and b) two thermocouples embedded in the large bore nylon 6.6 tube-beef assembly (thermocouples at radii r from the tube center).}
    \label{fig:BurgerRewarming}
\end{figure}

The experiment was repeated with the transducer coupled to a large bore nylon 6.6 tube containing 80\,mL of ground beef.
The temperature recordings from a central thermocouple and a thermocouple placed 6\,mm radially are shown in Fig. \ref{fig:BurgerRewarming}.
Similar behaviour is shown in both Fig. \ref{fig:BurgerRewarming}a) and b); the central thermocouple rewarms faster than the thermocouple placed radially.
The sub-zero rewarming rate at the central thermocouple was approximately 3 times slower slower than the central thermocouple in the transducer-frozen beef assembly rewarming experiment (which eliminated several reflective interfaces), and approximately twice as slow as in the cryovials, where thermal conduction sped up the early rewarming.

\subsection{Comparison of simulated and experimental rewarming}

Rewarming was simulated with a simplified 2D model and compared with the experimental temperature recordings, as shown in Fig. \ref{fig:USsim_vs_USexp}a).
The simulation replicates the experimental rewarming below 0$^{\circ}$C accurately, but the behaviour of the simulation above 0$^{\circ}$C deviates from the experimental measurements.
Figure \ref{fig:USsim_vs_USexp}a) shows a brief period at 35\,s where the simulated maximum temperature exceeds 0$^{\circ}$C, then cools again.
In simulation, this fast heating at 0$^{\circ}$C results from a brief but drastic increase in pressure and the corresponding volume rate of heat deposition.
The increase in pressure may be due to the beef sound speed reaching the perfect distribution to create a 444\,kHz resonant chamber.
This effect is not seen in experiment, possibly due to the assumptions and simplifications made in the simulation.
The simulation does not incorporate ground beef inhomogeneity, the acoustic properties of the beef around 0$^{\circ}$C may not be accurate, and shear waves that may arise in the experiment are not incorporated into the simulation.
The simulation also shows the beef rewarming to above 0$^{\circ}$C slightly earlier than was measured in experiment.
This may be the result of using a 2D simulation to reproduce a 3D experiment; the experimental rewarming setup uses a finite length tube, while the 2D simulation assumes the tube length is infinite and that heat only diffuses across the 2D plane, confining the heat in the third dimension.

\begin{figure}
    \centering
    \includegraphics[width = \columnwidth]{Figure11.eps}
    \caption{a) experimental vs. simulated rewarming (max., mean, and min. in a 1\,mm$^2$ area centered in the tube) b) 1D slices through the simulated cryovial for 80\,s of rewarming. The nylon 6.6 and beef edges are depicted with blue lines.}
    \label{fig:USsim_vs_USexp}
\end{figure}

The temperature slices from the rewarming simulations (Fig. \ref{fig:USsim_vs_USexp}b) show that below 0$^{\circ}$C, the rewarming of the periphery from thermal conduction and the central ultrasonic rewarming act together to reduce the temperature differential in the cryovial.
Before 35\,s, the maximum temperature differential within the cryovial was between 1$^{\circ}$C and 3$^{\circ}$C.
At 35\,s the simulated fast rewarming at the center of the cryovial results in a maximum temperature differential of 7$^{\circ}$C.
After 35\,s, when the mean temperature exceeds 0$^{\circ}$C, the absorption of ultrasound decreases which results in cooling at the centre of the cryovial and further rewarming towards the edges.
After 60\,s, the entire cryovial has thawed and the center of the cryovial again rewarms faster than its surroundings.

\section{Discussion}

The work presented here experimentally demonstrates the principle of ultrasonic rewarming and provides a foundation for further investigation and development of the technique.
Ultrasonic rewarming of frozen ground beef was demonstrated in 2\,mL cryovials at average rates over ten times faster than with thermal conduction alone.
The ultrasound system developed in this work to rapidly rewarm 2\,mL cryovials may be a useful tool to study the recovery of cryopreserved cells commonly stored in these sterile containers.
The transducer and nylon 6.6 setup was characterized with fibre-optic hydrophone measurements in the nylon tube cavity.
Simulation was then used to estimate the water to frozen beef pressure ratio, and it was estimated that the pressure amplitude at the centre of the cryovial during rapid rewarming of the beef was approximately 0.8\,MPa.
Ultrasonic rewarming was also demonstrated in two larger beef volumes where thermal conduction from the surroundings was minimized and faster ultrasonic rewarming was demonstrated at the transducer focus than at radial positions.
A hypothesized advantage of ultrasonic rewarming over radiofrequency-based methods is the decrease in absorption as the medium thaws, reducing the hot spot problem that can occur with radiofrequency-based methods.
The simulations in this work support this, showing a small temperature differential across the cryovial during the transition through sub-zero temperatures; this will be confirmed by future experimental work.

Prior work in ultrasonic rewarming has been limited to long timescales (tens of minutes to hours)\cite{miles1999high} or to simulations of rewarming below -20$^{\circ}$C that neglected the temperature-dependence of the thermal and acoustic parameters of biological tissues\cite{olmo2020use}.
This work experimentally demonstrates fast ultrasonic rewarming at high sub-zero temperatures, through the transition zone from frozen to thawed.
The mean rewarming rate measured at the centre of the cryovials was 23$^{\circ}$C per minute, and simulations suggest that the temperature gradient across the sample is of the order of a few degrees.
For slow frozen materials, slow warming rates of the order of 1$^{\circ}$C per minute are feasible, but increases in post thaw viability have been demonstrated where small volumes cell samples were warmed at faster rates in a water bath ($\sim$15$^{\circ}$C/min vs 0.6$^{\circ}$C/min)  \cite{kilbride2017cryopreservation}.
The faster warming rates measured here may further improve cell viability in these type of materials.
Previously, larger volumes of several litres have been rewarmed slowly (0.6$^{\circ}$C per minute) to reduce temperature differentials with acceptable 72\,h post-rewarming cell viability\cite{kilbride2017cryopreservation}. The addition of ultrasound could increase the rewarming rate at the center of these volumes, improving the bulk rewarming rate without the risk of generating a large temperature differential within the volume, consequently reducing the time needed for the cells to regain normal functionality.

For successful rewarming of vitrified media, warming rates of 100$^{\circ}$C/min or more (an order of magnitude faster than critical cooling rates) are commonly cited as necessary to avoid ice nucleation and growth\cite{olmo2020use,robinson2002electromagnetic,fahy2015principles}. However, this may be an underestimate and the critical rewarming rate may be multiple orders of magnitude faster than the critical cooling rate\cite{baudot2004thermal,fahy1987biological,fahy2015principles}.
With further optimisation of the ultrasonic warming device, it may be possible to increase rewarming rates towards the 50$^{\circ}$C/min generated by magnetic nanoparticle heating or further towards 100$^{\circ}$C/min which may be in the lower limit of rates required to warm vitrified materials.

To facilitate faster warming, transducer acoustic power output will be increased by optimising the efficiency of the source e.g. by electrical impedance matching and addition of an acoustic matching layer between the transducer and the nylon 6.6 holder.
The transducer frequency used here was chosen by balancing the focal dimensions with attenuation and lies within the 0.4-1\,MHz range of ultrasound frequencies previously found to be capable of slow ultrasound-based thawing\cite{miles1999high}. Together with increased knowledge of the frequency dependence of absorption in frozen materials, further optimisation of the source frequency may increase warming rates or ensure uniform heating of larger volumes.
The time course of energy delivery could also be adjusted to deliver a more constant rate of warming through the temperature range, reduce temperature gradients and reduce warming above 0$^{\circ}$C.
The use of nylon 6.6 to couple the transducer to the cryovials generated fast rewarming rates, but a deformable or conformal coupling layer may be needed to rewarm larger volumes where there are significant shape or volume changes during rewarming.
The tubular transducer geometry is well-matched to the cryovial geometry, but larger-volume ultrasonic rewarming will require focal steering achieved though the development of a phased array transducer\cite{olmo2020use}.

The relative contributions of scattering and absorption to attenuation in frozen biological materials are not known, but the experiments presented here demonstrate that there is some contribution of absorption to the attenuation that results in warming of the samples.
The rewarming simulations were performed with the acoustic absorption coefficient set to the measured attenuation coefficient of beef.
This is unlikely to represent the true level of absorption although the the simulation matches the experiment well.
This may be partly because the rewarming setup is effectively a closed acoustic system meaning that scattered sound will eventually be absorbed.
Further investigation of the contributions of absorption and scattering to acoustic attenuation of frozen biological materials is required for the further development of this technology.
Further characterisation of the acoustic and thermal properties of different cell and tissue types at sub-zero temperatures, as well as their dependence on cooling rate and mechanism (e.g. progressive solidification, network solidification, vitrification), and presence of cryoprotective and ice nucleating agents is also critical in understanding and optimising the delivery of energy for rewarming.
This information will enable more accurate modelling for development of more sophisticated ultrasonic rewarming platforms and optimisation of rewarming rates.

Ultrasonic rewarming research is in its early stages and there are many remaining parameters to investigate.
Rewarming of different samples of greater clinical relevance will be investigated, for example by examining rewarming rates in different cell types and their dependence on cooling rate and cryoprotectant type and concentration.
For example, DMSO, glycerol, and other cryoprotectants are used in concentrations between 10-60\% w/v depending on the cooling method, and are likely to influence the acoustic properties and consequent ultrasonic rewarming rates that can be achieved.
In future work, warming from the low (-140$^{\circ}$C to -80$^{\circ}$C) temperatures used for long-term cryopreservation will be investigated. It is likely that attenuation and absorption remain relatively constant at lower sub-zero temperatures, suggesting that fast ultrasonic rewarming can also be achieved at over this range.
Cell survival studies will be performed to evaluate the effect of rewarming rate on cell viability which is critical for optimising the rewarming and to investigate any other effects of exposure to ultrasound.
Lower `power' ultrasonic frequencies (20-100\,kHz) can generate ice nucleation\cite{chow2005study,dalvi2017review}, which when performed without the formation of large ice crystals can improve cell viability \cite{morris2013controlled}.
Further work is needed to determine whether higher frequency ultrasound generates ice nucleation and whether there is time for the ice crystals to grow during rapid ultrasonic rewarming.

\section{Conclusion}

A custom tubular ultrasonic transducer setup was built to demonstrate ultrasonic rewarming in 2\,mL cryovials.
The rewarming was performed from -20$^{\circ}$C to above 5$^{\circ}$C.
Experiments demonstrated the rewarming of frozen samples with sub-megahertz pressures, at a maximum rate of 57$^{\circ}$C per minute and at an average rate of over ten times faster than with thermal conduction alone.
The experiments demonstrated that absorption forms a part of frozen beef attenuation, and that ultrasound may be a viable method for rewarming cryopreserved samples.
Acoustic and thermal simulations matched well with experiment and suggest warming takes place with low thermal gradients across the sample.
This work presents a step towards using ultrasound to rewarm samples cryopreserved in cryovials, and develops the knowledge and tools that will be required to rewarm larger cryopreserved volumes.

\section{Acknowledgements}

This work was supported in part by a UKRI Future Leaders Fellowship [grant number MR/T019166/1], and in part by the Wellcome/EPSRC Centre for Interventional and Surgical Sciences (WEISS) (203145Z/16/Z).
For the purpose of open access, the author has applied a CC BY public copyright licence to any Author Accepted Manuscript version arising from this submission.
The authors wish to thank Professor Barry Fuller and Professor Clare Selden for useful discussions on this work.

Medical Physics (physics.med-ph); Biological Physics (physics.bio-ph)


X-ray emission from isolated neutron stars

Abstract: X-ray emission is a common feature of all varieties of isolated neutron stars (INS) and, thanks to the advent of sensitive instruments with good spectroscopic, timing, and imaging capabilities, X-ray observations have become an essential tool in the study of these objects. Non-thermal X-rays from young, energetic radio pulsars have been detected since the beginning of X-ray astronomy, and the long-sought thermal emission from cooling neutron star's surfaces can now be studied in detail in many pulsars spanning different ages, magnetic fields, and, possibly, surface compositions. In addition, other different manifestations of INS have been discovered with X-ray observations. These new classes of high-energy sources, comprising the nearby X-ray Dim Isolated Neutron Stars, the Central Compact Objects in supernova remnants, the Anomalous X-ray Pulsars, and the Soft Gamma-ray Repeaters, now add up to several tens of confirmed members, plus many candidates, and allow us to study a variety of phenomena unobservable in "standard'' radio pulsars.

\section{Introduction}
\label{sec:introduction}

With more than 1800 detections, rotation powered  pulsars (RPP) constitute
by far the largest class of isolated neutron stars (INS),
despite  only one out of $\sim$10 radio pulsars is visible because of beaming.
Accounting for the selection effects of radio observations,
a total population of $\sim10^6$  RPP is estimated for  the whole
Galaxy \cite{fau06}.
Observations at gamma-ray energy, where pulsar beaming angles are larger,
are now contributing to increase the number of known RPP \cite{abd09_blind}.
About one hundred  RPP have  been detected also at X-ray energies \cite{bec09}:
they include  the youngest and more energetic pulsars (like the Crab),
a few older neutron stars at  small distances (e.g. Vela and Geminga)
and several tens of recycled millisecond pulsars (most of which are found in globular clusters \cite{gri09}).

All kinds of INS, not only the RPP, are X-ray emitters.
X-ray observations have been crucial to discover other manifestations of INS,
that for various reasons were missed in the standard searches for
radio pulsars.  The nearby X-ray Dim Isolated Neutron Stars (XDINS),
the Central Compact Objects (CCOs) in supernova remnants, the Anomalous
X-ray Pulsars (AXPs), and the Soft Gamma-ray Repeaters (SGRs) are examples of these
new classes, which, despite totalling  only a few tens of sources, are
particularly  interesting because they offer a different view on a
variety of phenomena unobservable in "standard'' radio pulsars.

This review is focussed on the  main properties of the  X-ray emission from these
new classes of INS. For an excellent review of the X-ray emission from RPP see \cite{bec09}.
The location  on a
P--$\dot \textrm{P}$ plot diagram of some of the  objects discussed here is
shown in Fig. \ref{fig:ppdot}. The figure also gives  lines
of constant magnetic field (computed assuming dipole braking,
B = 3.2 10$^{19}$ (P $\dot \textrm{P}$)$^{1/2}$ G), characteristic age ($\tau_c$=$\frac{1}{2}$P/$\dot \textrm{P}$), and
spin-down luminosity (L$_{SD}$ = I $\Omega$ $\dot \Omega$ =
4 10$^{46}$ $\dot \textrm{P}$ / P$^3$ erg s$^{-1}$).

\begin{figure}[b]
\includegraphics[scale=.37, angle=0]{PPdot_revised.ps}
\caption{P--$\dot \textrm{P}$ diagram for different classes of INS:
XDINS ($\times$), CCO ($\diamond$), AXP (squares), SGR ($\bigtriangleup$),
RRAT (+). The vertical bars indicate the range of $\dot \textrm{P}$ variations
observed in SGR and AXP.
Lines of constant magnetic field (dashed), characteristic age (dotted), and spin-down
luminosity (dash-dotted) are also indicated.}
\label{fig:ppdot}
\end{figure}

\section{Origin of the X-ray emission in isolated neutron stars}
\label{sec:ori}

The X-ray emission observed from INS can be powered by internal heat, rotational energy, accretion, and
magnetic field decay. The relative importance of these energy sources,
that can also operate at the same time, depends on the age and  physical properties of the neutron star.

Neutron stars have internal temperatures of $\sim10^{11}$ K at birth,
that rapidly drop to $\sim10^{9}$ K.
For the following $\sim10^5-10^6$ years  the dominant cooling mechanism
is neutrino emission from the star's isothermal core.
This leads to surface temperatures of several $10^5$ to $10^6$ K, with  thermal emission peaking in the soft X-ray band.
Temperature gradients on the star's surface generally produce observable modulations at
the rotation period. Thermal X-ray emission
can be observed in INS with ages of $\sim10^4-10^6$ years,
provided they are sufficiently close and not too absorbed by the interstellar medium.
Older neutron stars are too cool to significantly emit  X-rays,
while in the youngest pulsars the thermal radiation is difficult to detect
because it is outshined  by the brighter non-thermal emission. A recent review on the thermal emission
from neutron star is given by \cite{zav09}.

Non-thermal emission, that extends over a broad energy range, originates from
charged particles accelerated in the NS magnetospheres at the expense
of rotational energy (see, e.g., \cite{che09}).
Non-thermal  X-rays are characterized by power-law spectra and  strongly
anisotropic emission patterns, giving rise to large pulsed fractions.
The pulse profiles often show narrow (double) peaks, but in many cases nearly
sinusoidal profiles are observed.
The most luminous RPP are the Crab and two young pulsars in the
Large Magellanic Cloud
(the only three pulsars with L$_{SD}>10^{38}$ erg s$^{-1}$).
The efficiency with which the rotational energy is converted to non-thermal luminosity
is about 10$^{-3}$ \cite{bec97,li08}, but there is a
large dispersion around the average value \cite{pos02}, as expected because of
different viewing orientations and, possibly,  also other effects.
Studies of the L$_X$--L$_{SD}$ relation with a large sample are complicated  by the
fact that a significant fraction of the rotational energy
loss powers pulsar wind nebulae,
which are difficult to disentangle in the more distant and/or fainter objects
without adequate spatial resolution.

Accretion is a well established process in X-ray binaries (e.g., \cite{fra02}).
In the lack of a companion star,
accreting matter could originate directly from the interstellar medium.
However, the relatively large space velocity of neutron stars, and the low density of the
interstellar medium, make this process unable to provide sufficiently high luminosities.
Alternatively, the  matter could  be supplied by a debris disk
formed by fall-back in the supernova explosion that produced the neutron star.
Although this possibility  seems more promising,
no unambiguous  evidence for an INS powered by accretion has been found yet.

The relevance of magnetic energy in powering the emission from  neutron stars has
been recognized only recently, with the observation of SGRs and
AXPs (see \cite{woo06, mer08} for reviews). High magnetic fields (B$\sim$10$^{14}$--10$^{15}$ G)
were first invoked to explain the SGRs and, in particular, the unique
properties of the exceptional event of March 5, 1979 from SGR~0526$-$66\ \cite{pac92, dun92}.
Objects in which magnetic field decay is the
dominant energy source have been named Magnetars. They are particularly interesting
because they offer the unique possibility to study physical processes
in magnetic fields of unequalled strength.

\begin{table}
\caption{Isolated Neutron Stars}
\label{tab:list}
\begin{tabular}{lccccl}
\hline\noalign{\smallskip}
Name  & P & $\dot \textrm{P}$       & D$^{(a)}$              &  Association & Notes$^{(b)}$  and \\
      &  (s)  &  (s s$^{-1}$)  & (kpc) &   &    References  \\
\noalign{\smallskip}\svhline\noalign{\smallskip}
 \multicolumn{6}{c}{XDINS} \\
 \hline
RX~J0420.0$-$5022\ & 3.45  &  -    &  0.345   & - &  \cite{hab04b} \\[3pt]
RX~J0720.4$-$3125\ & 8.39  & 7.01 10$^{-14}$ &  0.36    & - &   G? \cite{hab97,kap07,van07,hoh10} \\[3pt]
RX~J0806.4$-$4123\ & 11.37 & (5.5$\pm$3.0) 10$^{-14}$  &  0.25    & - &  \cite{hab04b,kap09b} \\[3pt]
RX~J1308.6$+$2127\  & 10.31 & 1.120 10$^{-13}$  & 0.5     & - & RBS1223 \cite{ham02,kap05a} \\[3pt]
RX~J1605.3$+$3249\  &   -   &   -     & 0.39   & - & RBS1556 \cite{mot99,zan06} \\[3pt]
RX~J1856.5$-$3754\  & 7.06 & (2.97$\pm$0.07) 10$^{-14}$ & 0.16   & - &  \cite{wal01,tie07,van07a,van08} \\[3pt]
RX~J2143.0$+$0654\  &  9.43 & (4.1$\pm$1.8) 10$^{-14}$  & 0.43 & - & RBS1774 \cite{zam01,kap09d} \\[3pt]
2XMM J104608.7--594306    &  -     &  -      &  2   & - & candidate \cite{pir09} \\[3pt]
RX J1412.9+7922 &  -    &  -       &  3.6 &   - & candidate, Calvera \cite{rut08,she09}\\[3pt]
\hline
 \multicolumn{6}{c}{CCO} \\
 \hline
RX J0822.0--4300 & 0.122 & $<$8 10$^{-15}$& 2.2 & Puppis A & \cite{pet96,got09} \\[3pt]
CXOU J085201.4--461753 & -      & -       &   1  & G266.1--1.2  & \cite{pav01} \\[3pt]
1E 1207.4--5209 & 0.424 & $<$2.5 10$^{-16}$& 2 & G296.5+10.0 &  \cite{mer02,del04,got04} \\[3pt]
CXOU J160103.1--513353 &  -  &  - & 5  & G330.2+0.1 &  \cite{par06} \\[3pt]
RX J1713.4--3949  & - & - & 1.3 & G347.3--0.5 & \cite{laz03} \\[3pt]
CXOU J185238.6+004020  & 0.105 & 8.7 10$^{-18}$ & 7.1 & Kes 79 & \cite{hal07,hal10} \\[3pt]
1E 161348--5055.1  & -  & -  & 3.3  & RCW 103 &   T \cite{del06}\\[3pt]
CXOU J232327.8+584842  &  - &  - &  3.4 & Cas A & \cite{cha01,mer02c,pav09} \\[3pt]
XMMU J172054.5--372652  & - & - &   4.5 &  G350.1--0.3 & candidate \cite{gae08} \\[3pt]
XMMU J173203.3--344518  & - & - &   3.2 & G353.6--0.7 & candidate \cite{tia10}\\[3pt]
CXOU J181852.0--150213  & - & - &  8.5 & G15.9+0.2 & candidate \cite{rey06b} \\[3pt]
\hline
 \multicolumn{6}{c}{AXP and SGR} \\
 \hline
CXOU~J010043.1$-$721134  & 8.02 & 1.9 10$^{-11}$ &60 & SMC &  AXP \cite{lam02,mcg05,tie08} \\[3pt
4U~0142$+$61  & 8.69 & 2 10$^{-12}$&  3.6     &  - &    B, G? \cite{isr94,dib07b} \\[3pt]
1E~1048.1$-$5937 & 6.45 & (1--10) 10$^{-11}$ &  8 & - &   B, G \cite{sew86,tie05a,gav02,dib09} \\[3pt]
1E 1547.0$-$5408  & 2.07 & 2.3 10$^{-11}$ &5 & G327.2--0.1 &   T, B, R \cite{cam07c,hal08,esp08b,mer09}     \\[3pt]
PSR J1622--4950 & 4.33 & 1.7 10$^{-11}$ & 9 & - &   R  \cite{lev10} \\[3pt]
CXOU~J164710.2$-$455216  & 10.6 & 9.2 10$^{-13}$ & 3.9 & Westerlund 1 &    T, B, G  \cite{mun07,isr07b} \\[3pt]
1RXS~J170849.0$-$400910  & 11.0 & 2.4 10$^{-11}$  & 5& - &   G \cite{sug97,isr99,dal03,dib07a} \\[3pt]
XTE~J1810$-$197 & 5.54 &(0.8--2.2) 10$^{-11}$&3.1 & - &   T, B, R \cite{ibr04,cam06,woo05,ber09} \\[3pt]
1E~1841$-$045 & 11.77& 4.1 10$^{-11}$ & 8.5 & Kes 73 &   G \cite{vas97,got02,mor03,dib07a} \\[3pt]
1E~2259$+$586     & 6.98 & 4.8 10$^{-13}$&4 &CTB 109&    B, G \cite{fah81,kas03} \\[3pt]
AX J1844.8--0256 & 6.97 & - & 8.5 & G29.6+0.1 & candidate, T \cite{vas00,tam06}\\[3pt]
\hline
SGR~0418$+$5729  &  9.1 & $<$1.1 10$^{-13}$     &   2   & -  &   T, B \cite{van10,esp10} \\[3pt]
SGR~0501$+$4516  &  5.76 & 7.1 10$^{-12}$   & 1.5 &   - &   T, B \cite{rea09,apt09,eno09} \\[3pt]
SGR~0526$-$66  & 8.05 & 6.5 10$^{-11}$ & 55 & LMC, N49&  B, GF \cite{maz79,maz99b,tie09} \\[3pt]
SGR~1627$-$41 & 2.59 & 1.9 10$^{-11}$  &    11 & - &   T, B \cite{woo99c,mer06a,esp09a,esp09b}      \\[3pt]
SGR~1806$-$20&7.6&(8--80) 10$^{-11}$&8.7&star cluster&  B, GF  \cite{kou98,pal05,mer05c,woo07} \\[3pt]
SGR 1833--0832 &  7.6  & 7.4 10$^{-12}$  &   10    & -  &   T, B \cite{gog10,esp10b} \\[3pt]
SGR~1900$+$14&5.2 &(5--14) 10$^{-11}$ &15&star cluster &  B, GF, G? \cite{maz79b,woo99b,mer06b}\\[3pt]%
\hline
\noalign{\smallskip}\hline\noalign{\smallskip}
\end{tabular}

$^{(a)}$ in several cases the distances have large uncertainties; the values
adopted in the figures of this work are indicated here.

$^{(b)}$ B = bursts, G = glitches, GF = giant flares, R = radio emission, T = transient

\end{table}

\section{The X-ray Dim Isolated Neutron Stars}
\label{sec:xdins}

A handful of X-ray sources with
large X-ray-to-optical flux ratios,  F$_{x}$/F$_{opt}\sim$10$^4$--10$^5$,
typical of INS, were discovered in the ROSAT satellite All Sky Survey.
Their INS nature was confirmed by the measurement of a large proper
motion in the brightest source, RX~J1856.5$-$3754\ \cite{wal01}, and by the discovery
of pulsations at few seconds in other sources \cite{hab97,ham02}.
The original members of this class\footnote{they have been nicknamed ''The Magnificent Seven'' (M7, hereinafter)} and two new candidates
are listed in Table \ref{tab:list}. Recent reviews on the XDINS are given in  \cite{hab07,tur09}.

XDINS have very soft thermal spectra (blackbody temperatures T$_{BB}\sim$40--110 eV), X--ray luminosities
L$_X\sim10^{30}-10^{32}$ erg s$^{-1}$, spin periods in the 3--12 s range,
faint optical counterparts (V$>$25), and no radio emission.
Period derivatives of the
order of 10$^{-14}$--10$^{-13}$ s s$^{-1}$ have been   measured for several XDINS through
phase connected timing, but in some cases these values are   still
poorly constrained   (see Table \ref{tab:list}).
The  temperatures of XDINS
are consistent with neutron stars cooling curves if they have ages
of 10$^5$--10$^6$ years and  the effect of their
strong magnetic field  is taken into account \cite{agu08}.
It is thus generally believed that the XDINS  are powered by residual thermal energy.

X-ray sources with such soft spectra can only observed if the interstellar
absorption is small:  all the M7 are closer than $\sim$0.5  kpc and
have  N$_H<$4 10$^{20}$ cm$^{-2}$.
Their proximity allowed in several cases the  measurement of parallax and
and proper motion
showing a distribution of  transverse velocities  consistent
with that of radio pulsars \cite{mot09}.
For these velocities, and typical ISM densities, accretion from the ISM, originally invoked
to explain the XDINS emission, cannot  provide the observed luminosity.
Although the timing parameters give spin-down ages
$\tau_c$$\sim$ 1--2 million years,
the measured spatial velocities and likely birth places in nearby OB associations
strongly suggest that their true ages are smaller (see Table \ref{tab:ages}).

\begin{table}
\caption{Comparison of estimated and spin-down ages of INS} \label{tab:ages}
\begin{tabular}{lccccl}
\hline\noalign{\smallskip}
Source  & Spin-down   &  Estimated      &  Method  \\
      &  age ($\tau_c$ )  &      age &     \\
\noalign{\smallskip}\svhline\noalign{\smallskip}
RX~J0720.4$-$3125\  & 1.9 Myr  &  (0.7$^{+0.3}_{-0.2}$) Myr  & proper motion  \\[3pt]
RX~J1308.6$+$2127\  & 1.5 Myr  & (0.5--1.4) Myr & proper motion  \\[3pt]
RX~J1856.5$-$3754\ & 3.8 Myr & 0.4 Myr  & proper motion \\[3pt]
\hline
1E~1841$-$045    &  4.5 kyr   &  0.5--1 kyr   & SNR age      \\[3pt]
1E~2259$+$586     & 230 kyr & (11.7$\pm$1.2) kyr   & SNR age    \\[3pt]
SGR~0526$-$66    &  2 kyr   &  5 kyr  & SNR  age  \\[3pt]
\hline
CXOU J085201.4--461753 & $>$240 kyr &  3.7 kyr &  SNR  age   \\[3pt]
1E 1207.4$-$5209 & $>$27 Myr & 7 kyr &   SNR  age \\[3pt]
CXOU J185238.6+004020  & 190  Myr  &  7 kyr & SNR  age  \\[3pt]
\hline
\noalign{\smallskip}\hline\noalign{\smallskip}
\end{tabular}

\end{table}

Despite various attempts, none of the XDINS has been detected
in the radio band. The most recent  upper limits  of  0.14--5 $\mu$Jy kpc$^2$,
obtained at 1400 MHz \cite{kon09}, correspond to pulsed luminosities
well below those of the  faintest observed radio PSRs.
However, given the small sample and considering the anticorrelation
between beaming aperture angle and pulse period, it cannot be excluded that the lack of radio emission be
simply due to unfavorable orientations.

The X-ray light curves of XDINS are nearly sinusoidal, with  pulsed fractions,
ranging from 1.5\% to $\sim$20\%, and  moderate energy dependence. Only
in the  case of RX~J1308.6$+$2127\ two peaks are clearly seen. These properties indicate that
we are likely seeing emission from a large fraction of the star's surface.

The first  spectra of XDINS, obtained with ROSAT, could be well described with blackbody curves.
With the better data of XMM and Chandra, significant deviations from this simple model have been detected in most XDINS (the only  exception is  RX~J1856.5$-$3754\ \cite{bur03}).
The features consist of broad absorption lines  in the energy range $\sim$0.2--0.8 keV
(see \cite{van07a,hab07} and references therein).
Multiple lines, with energy spacings  consistent with harmonic ratios,
are present in nearly half of the sources.
They have been interpreted as proton cyclotron lines or
as H and/or He atomic transitions, falling in the soft X-ray range
due to the effect of a high magnetic field on the binding energies.
Both interpretations require   B$\sim10^{13}-10^{14}$ G,
broadly consistent with the values of 2.5 and 3.4 10$^{13}$ G derived from the timing
parameters of  RX~J0720.4$-$3125\ and RX~J1308.6$+$2127\ assuming  dipole
braking\footnote{Note that the lines measure the surface field, which can be higher
than the dipole component  dominating at large radii and responsible for the spin-down.}.

Variability on a time scale of few years was discovered in
RX~J0720.4$-$3125\ \cite{dev04}, a target regularly observed for XMM calibrations
because it was believed to be a constant source.
If we exclude AXPs and SGRs, this has been the first INS  to show
significant changes in its X-ray flux, spectrum and pulse profile.
Based on the first observations,
it was suggested that the variations
had  a periodicity of 7.1$\pm$0.5 years  and could be caused by precession of the
neutron star \cite{hab06}.
However, the most recent analysis, which include new observations and
take into account the variable phase shifts between hard and soft energy bands, do not
support the periodicity  \cite{hoh09}.
Most of the spectral variation occurred in a relatively short timescale
of half a year \cite{van07}, in coincidence with a possible glitch.
A change in the temperature and/or composition of part of the NS surface might
have been caused by the energy released in the glitch, or alternatively
by a sudden accretion episode.

The other XDINS seem to be constant in flux, but small variations might have been unnoticed
since, with the exception of RX~J1856.5$-$3754\  ,  they have not been observed very often with sensitive instruments.
On the other hand, RX~J1856.5$-$3754\ is twice as bright as RX~J0720.4$-$3125\ and  has  been observed repeatedly,
but no variability has been reported. Interestingly this is the only XDINS without spectral lines
\cite{bur03} and with the smallest pulsed fraction ($\sim$1.5\% \cite{tie07}).

Until recently, attempts to enlarge the XDINS sample had little success.
Ongoing searches in the  ROSAT Bright Sources Catalogue are reported in
\cite{tur10}. Looking for new INS among the tens of thousands serendipitous
XMM and Chandra sources requires a long and difficult multi-wavelength effort,
also considering that they are fainter and more distant than the bright M7.
One promising source has been   selected
among XDINS candidates discovered with XMM \cite{pir09}, but its temperature
(kT$_{BB}$=117 eV) and  estimated luminosity ($\sim10^{33}$ erg s$^{-1}$) are larger
than those of the M7.
The same is true for another candidate, RX J1412.9+7922  \cite{rut08},
for which other interpretations  (e.g. a nearby millisecond pulsar) cannot be excluded yet.
Its blackbody temperature is kT$_{BB}\sim$200 eV. A   hydrogen
atmosphere model (kT$\sim$122 eV), with a possible emission line at 0.53 keV,
gives a better fit to the Chandra data \cite{she09}.
If the emission comes from the whole star's surface, its distance would be 3.6 kpc,
while a smaller emission region  and  distance  should give rise to
pulsations.  Only the detection of a periodicity could clarify the real
nature of these new XDINS candidates.

\section{ Central Compact Objects in Supernova Remnants}
\label{sec:cco}

The group of Central Compact Objects (CCOs) consists of several radio quiet
X--ray sources located at the center of shell-like
SNRs and with  high F$_x$/F$_{opt}$\cite{pav04,del08}.
The association  with SNRs (see Table \ref{tab:list})
imply ages at most of a few tens kyrs, and gives the possibility  to know their distances.
Three CCOs show X-ray pulsations with periods in the 0.1--0.4 s range,
consistent with  young neutron stars,
but with unexpectedly  small  spin-down rates  \cite{got07,hal10}.
The extremely small $\dot \textrm{P}$  measured for the  CCO in  Kes 79
implies a dipole magnetic field of only 3.1 10$^{10}$ G, while the
upper limits for the sources in G296.5+10.0 and Puppis A
give B$<$3.3  10$^{11}$ G and  B$<$10$^{12}$ G,   respectively.
The  spin-down ages of these three CCOs exceed by orders of magnitude the  ages
of their associated SNRs (see Table \ref{tab:ages}), indicating that these INS were probably
born with spin periods very close to the current values.
It is likely that their relatively long initial spin period and low
magnetic field be causally connected. This has led to the ''anti-magnetar''
nickname for these objects.
The three ''anti-magnetars'' show  a  variety of different  pulse profiles: the CCOs
in Kes 79 has a pulsed fraction of  $\; \buildrel > \over \sim \;$ 60\%,  one of the highest among thermally
emitting INS \cite{hal07}. A more standard value of \ltsima10\%
is shown in 1E 1207--5209, but with the
peculiarity that most of the pulsation can be attributed to phase variations in the
shape of its prominent absorption lines at 0.7, 1.4, 2.1 and possibly 2.8 keV \cite{del04}.
The pulsations in the Puppis A CCOs eluded an earlier detection because
of a phase shift of 180$^{\circ}$ between nearly sinusoidal profiles in the
soft ($<$1.2 keV) and hard bands \cite{got09}.

Most CCOs have thermal-like X-ray emission, with  blackbody temperatures in the
range $\sim$0.2--0.5 keV,
and no evidence for additional non-thermal components.
The rotational energy loss of the three pulsed COOs is certainly too small to give
a detectable contribution to their observed X-rays, assuming a typical efficiency
of 10$^{-3}$. The emission could be due to residual cooling and/or weak accretion
from a residual disk.

When data of good statistical quality are available,
single-temperature blackbody  spectra do not provide acceptable fits.
Better fits are obtained with the sum of two blackbodies  with kT$_1$
in the range 0.16--0.4 keV and kT$_{2}\sim$2 kT$_{1}$.
This simple model can be thought as
a first approximation of a non-uniform temperature distribution. The hotter blackbody
is energetically significant:
the bolometric luminosities of the two components are in
the ratio  L$_2$/L$_1$$\sim$0.2--0.7.
The two-blackbody model implies small X-ray emitting areas (R$_1$$\sim$0.4--4 km,
R$_2$$\sim$0.06--0.8 km),
which are difficult to explain in weakly magnetized neutron stars.
The situation does not  significantly change if  atmosphere models, which
give lower temperatures and larger radii, are used instead of the blackbody fits.

It is unclear whether also the CCOs with unknown spin period
can  be interpreted as ''anti-magnetars''.
The lack of pulsar wind nebulae and of radio/gamma emission suggests that
also these  COOs have a small L$_{SD}$, but this can result equally well  from
a small $\dot \textrm{P}$ or from a long spin period.
Thus it cannot be excluded that CCOs constitute a
heterogenous class comprising objects of different kinds.
The source
at the center of the $\sim$330 years old  Cas A  SNR is one of the most studied CCOs.
It has a very soft  spectrum and   most spectral fits, with a variety of thermal models, indicate
an emitting area smaller than the whole surface of a neutron star.
No pulsations have been detected, down to quite stringent limits \cite{mer02c,pav09,hal10},
but an anti-magnetar interpretation remains plausible, considering the possibility of an unfavorable
orientation and the variety of pulsed fractions seen in other INS.

On the other hand,  the CCO  in RCW 103, with a clear periodicity
at 6.67 hours and long term variations spanning more than two orders of magnitude
\cite{del06}, has unique properties
requiring an ad-hoc explanation. The periodicity, if interpreted as an
orbital motion, suggests a low mass binary \cite{del06,bha09},
but the deep optical limits pose severe constraints on the companion star \cite{del08b}.
Models involving magnetars, isolated or in a synchronous binary, and
fall-back disks have also been proposed \cite{piz08,li07}.

\section{The Magnetar candidates: Anomalous X-ray Pulsars and Soft Gamma-ray Repeaters}
\label{sec:magnetars}

Anomalous X-ray Pulsars (AXPs) and  Soft Gamma-Ray Repeaters (SGRs) are
spinning-down pulsars characterized by a luminosity in the soft and hard X-rays
larger than their available rotational energy loss and by the emission of bursts
and flares.

The first AXPs, discovered as bright pulsars in the soft X-ray range ($<$10 keV),
were initially classified as X-ray binaries powered by accretion.
Further X--ray data, coupled to deep  searches for optical/IR
counterparts, revealed their different nature \cite{mer95}. In
particular, these observations showed that the narrow period
distribution, long-term spin-down, and soft spectrum of the AXPs
were at variance with the properties of the larger
population of pulsars in massive binaries. Furthermore, no signs of binary companions
could be found in the AXPs \cite{mer98}.

The SGRs were instead discovered through the detection of bright and short bursts in the hard X-ray/soft
gamma-ray range, and initially considered  a particular subclass of gamma-ray bursts
\cite{lar86,att87}, with the notable property of ``repeating''
from the same sky direction. When accurate localizations became
available, it was possible to identify the X-ray counterparts of
SGRs, finding that they are pulsating sources very similar to the AXPs.

Although historically divided in two classes, many similarities indicate
that AXPs and SGRs are probably the same kind of astrophysical  objects
(see \cite{woo06,kas07,mer08} for reviews).
The most successful model to explain AXPs and SGRs  involves
highly magnetized neutron stars, or ``magnetars'' \cite{dun92,
tho95, tho96}, but other possibilities have been
proposed, e.g., models based on INS accreting
from residual disks \cite{cha00a,per00b,alp01,ert09}, or different kinds of
quark stars \cite{xu07,hor07,ouy07a,cea06}.

Magnetars differ from the normal pulsars, not only
for their higher field intensity, but also because their magnetosphere is not
dipolar. It is believed to consist of a twisted dipole, i.e. a field with  a significant
azimuthal component. This causes the presence of large-scale currents
with high charge density flowing in the magnetosphere and
affecting the emerging spectra by resonant scattering.

At energy below 10 keV, the X-ray spectra of AXP and SGRs have been traditionally described with
a two-components model consisting of a power-law plus a blackbody with kT$_{BB}\sim$0.5 keV.
The steepness of the power-law component (especially in the AXPs,
photon index $\sim$3--4) requires  N$_H$ values higher
than those independently estimated in other ways and leads to
overestimate the flux of the near infrared and optical counterparts (unless a
drastic, and possibly un-physical, spectral cut-off is invoked).
This problem is solved by adopting a two-blackbody model, that
gives equally good fits \cite{hal05}. XMM observations of CXOU~J010043.1$-$721134\ in  the Small
Magellanic Cloud, less affected by the interstellar absorption than the
Galactic AXPs/SGRs, have clearly demonstrated  that the power-law plus
blackbody model can be rejected with high confidence, while a good
fit to the soft X-ray spectrum is obtained with the sum of two blackbodies \cite{tie08}.

Magnetars have  been detected also in the hard X-ray range ($\; \buildrel > \over \sim \;$ 20 keV), with INTEGRAL and RXTE.
Hard power-law tails, extending to $\sim$150 keV, have been seen both in AXPs \cite{kui04,den08a,den08b}
and SGRs \cite{mer05a,goe06b,esp07,rea09}.
The hard X-ray emission is pulsed, and gives a luminosity  comparable to, and in some cases higher than,
that in the soft X-ray band.

Considerable effort is ongoing in the analysis and interpretation of the AXPs and SGRs spectra,
in order to progress from purely phenomenological fits  to more physical models.
The presence of a relatively dense plasma in magnetospheres with a twisted configuration \cite{tho02}
is expected to affect the emergent spectrum through resonant cyclotron
scattering  (RCS) of the thermal photons emitted by the underlying
neutron star's surface. A simplified model, based on  a
semi-analytical treatment of these effects \cite{lyu06b}, has been
systematically applied to several magnetar candidates
\cite{rea08}. This showed that at E$<$10 keV the RCS model can replace
the blackbody plus power-law model (in most sources an additional
power-law is still required to fit the hard X-ray tails). More
realistic 3-D simulations have also been performed to compute the
magnetar's spectral models \cite{fer07,nob08a,nob08b} and quite
successfully applied to the spectra of several magnetar candidates
\cite{zan09}.

\begin{figure}[b]
\includegraphics[scale=.4, angle=-90]{magnetar_outbursts_last.ps}
\caption{X-ray light curves of  outburst decays observed
in several AXPs and SGRs (from \cite{rea09}).}
\label{fig:dec}
\end{figure}

Two possibilities have been proposed to explain the
high-energy tails in the context of the twisted
magnetosphere model  \cite{tho05}:
bremsstrahlung from a thin turbulent layer of
the star's surface heated to kT$\sim$100 keV by magnetospheric
currents and synchrotron emission from mildly relativistic pairs
produced at a height of $\sim100$ km above the neutron star.
Compton up-scattering of soft thermal photons, emitted from the
star surface, by relativistic electrons threaded in the
magnetosphere has  also been
considered \cite{har91,gon00,bar07}.

Compared to all the other INS, AXPs and SGRs stand apart for their
striking variability properties.
I refer not only to the  short bursts and  (giant) flares, but
also to the so called ''persistent'' emission,
which for historically bright sources varies around  10$^{36}$ erg s$^{-1}$,
while for newly-discovered transients has been seen to span up to five
orders of magnitude.
About  half of the AXPs/SGRs have been always seen at relatively high
flux levels since the time of their discovery.
They were thought to be steady sources, but, as more sensitive imaging observations
became available, it appeared that variations  of $\sim$50\%
on timescales from days to months are quite common.
Long term variations are sometimes correlated with changes in the spin-down properties.
This is not unexpected in the magnetar model involving variations in the magnetospheric
twist angle.
For example, the luminosity variations, the correlation between spectral hardness and spin-down
rate, and the increase in the bursting activity observed before the 2004
giant flare of SGR~1806$-$20\ \cite{mer05c}
are  consistent with a gradually increasing magnetospheric twist angle.
This gives rise to a larger optical depth for resonant scattering (causing a
hardening of the X-ray spectrum),  a  higher rate of crustal fractures
(responsible for the bursts), and  a faster spin-down
(because of the larger fraction of field lines that open out across
the speed of light cylinder).
The spectral softening and the decrease in the flux
and spin-down rate observed after the giant
flare \cite{rea05b,tie05b} are also consistent with this picture, suggesting a
twist angle reduction as consequence of the global magnetic rearrangement
following the giant flare.

The first transient AXP, XTE~J1810$-$197 ,  was serendipitously discovered during the decay phase
of an outburst in 2003 \cite{ibr04}. Comparison of the highest observed luminosity
with archival data, showing a  faint soft source with L$_{x}\sim10^{33}$ erg s$^{-1}$,
indicated a dynamic range of two orders of magnitude.
Several other transient AXP/SGR have been discovered after XTE~J1810$-$197\ (see Table \ref{tab:list}).
When bright, their spectral and timing properties are similar to those of the
persistent sources.
During their low (or ''quiescent'') states  they have luminosity
of $\sim10^{32}$ erg s$^{-1}$ and soft thermal spectra, that make them similar to the CCOs.

Transient AXPs/SGRs now outnumber the persistent members of the class, and their
number will certainly increase in the future.
It is  difficult to estimate how many AXP/SGRs are present
in the Galaxy, because the duty cycle of transient magnetars is poorly known.
At least two sources have shown more than one outburst
(SGR~1627$-$41\ in 1998 and in 2008 \cite{esp08};
1E 1547.0$-$5408\ in Summer 2007 \cite{hal08} and 2008--2009 \cite{mer09}),  but other events might have been missed.
In fact the rate of discoveries of new transients has increased in the latest years,
thanks to the presence of satellites with large field of view and continuous sky monitoring,
like Fermi and Swift.  These new data revealed a variety of behaviors for what concerns
the outburst evolution of the different sources, as shown by the examples plotted in Fig. \ref{fig:dec}.
Poorly sampled light curves can be roughly described by power-law decays  with a variety of slopes,
while the well sampled  light curves show that the time decay is initially steeper and
flattens after about one day. This could indicate the presence of two different mechanisms operating
on different timescales, e.g. dissipation of magnetospheric currents followed by cooling of the star's
crust. Re-brightenings, or ''bumps'', superimposed on the decaying trends have also been observed, for
example in SGR~1900$+$14\ \cite{fer03}.

The observation of outburst decays, and, more in general, of the spectral and flux variability of transient
AXPs/SGRs, has the potential to greatly advance our understanding of neutron stars.
For example,  detailed studies of the long term spectral evolution have been reported
for XTE~J1810$-$197\ \cite{got05,ber09}.   However some
caveats should be remembered when interpreting the observational data.
Often the observations have a poor and irregular time sampling, especially in the initial and possibly
more variable phases, and/or are obtained with different instruments.
This reduces the possibility of determining correlations between the different energy
ranges in a robust way.
At late times, when the fluxes are small, spectral uncertainties can be significant,
also in view of the bolometric corrections required by the  high interstellar absorption.
Furthermore, the observed fluxes most likely result from the superposition of
various physical components with different variability behaviors and difficult to disentangle.
Part of these problems will be reduced when broad-band, sensitive instruments
able to react in a short time to the discovery of new transient events will become available.

\section{Rotating Radio Transients}
\label{sec:rrat}

New techniques of analysis for transient radio signals, applied to data
of the Parkes Multibeam Survey,  have led to the discovery
of  Rotating Radio Transients (RRATs) \cite{mcl06}:  neutron stars which emit
short (2-30 ms) radio pulses at intervals of minutes
to hours. Their rotation periods, ranging from  0.4 to 7 s, have been inferred from the largest common divisors of the time intervals between bursts. RRATs are probably a large galactic
population, that remained undiscovered for a long time due to lack of adequate radio searches,
but their relation to other classes of INS is unclear (see \cite{mcl09b} for a review).

The radio properties of RRATs show analogies with those observed in some peculiar
radio RPP, and they might represent extreme cases of some of these phenomena.
On the other hand, their long rotation periods, and the fact that transient
radio emission has been observed in a couple of AXPs \cite{cam06,cam07c}, suggest a possible
relation with magnetars.
The  period derivatives, determined to date for only $\sim$6--7 of the nearly
twenty known RRATs \cite{mcl09}, indicate  magnetic fields above 10$^{13}$ G
only in about half of the cases (see Fig. \ref{fig:ppdot}).

PSR~J1819$-$1458\ is the only RRAT that has been detected at X-ray energies \cite{rey06}.
It has a luminosity of $\sim$(2--5) 10$^{33}$ erg s$^{-1}$  (for d=3.6 kpc) and its
X-ray emission is pulsed at the rotation period of 4.26 s. The   light curve is nearly
sinusoidal, with a pulsed fraction of about 35\%. The spectrum,   well fit by a
blackbody (T$_{BB}\sim$0.14 keV) plus  a possible feature around 1 keV
\cite{mcl07},  indicates that the emission is of thermal origin (the luminosity
is subject to the distance uncertainty, but most likely larger than L$_{SD}$=3 10$^{32}$ erg s$^{-1}$).
A Chandra observation has shown evidence for diffuse X-ray emission surrounding PSR~J1819$-$1458\  \cite{rea09b}.
If this nebula is powered by the rotational energy of the RRAT, an
efficiency  of $\sim$20\%, much higher than that of all the other pulsar wind nebulae, is required.
This discrepancy could be reduced by a smaller distance, or by invoking a
shock caused by a large spatial velocity of PSR~J1819$-$1458 , but  the possibility that
magnetic energy contributes to power the observed diffuse emission has also been proposed \cite{rea09b}.

X-ray  upper limits at a level much smaller than the flux of PSR~J1819$-$1458\ have been reported for other two RRATs \cite{kap09c}.
If the X-ray emission from RRATs is simply due to cooling, this could be caused by a larger
age of these two objects (their spin-down ages are 0.4 and 0.8 Myr,
compared to  $\tau_{SD}$=0.1 Myr   for PSR~J1819$-$1458\ ). However, the significant uncertainties on distance
and absorption prevent a firm conclusion.

\begin{table}
\caption{Rotation Powered Pulsars with B$>$4 10$^{13}$ G.}
\label{tab:highB}
\begin{tabular}{lccccl}
\hline
\noalign{\smallskip}
Name  & Period  & $\dot \textrm{P}$       & B   & L$_x$  &   \\
      &  (s)  &  (s s$^{-1}$)  & (G) &  (erg s$^{-1}$) &     \\
\noalign{\smallskip}\svhline\noalign{\smallskip}
 J1847--0130~~ &  6.7   & 1.3 10$^{-12}$ & ~~9.4 10$^{13}$~~ & $<$5 10$^{33}$ &   \\[3pt]
 J1718--3718~~ &  3.4   & 1.6 10$^{-12}$ & ~~7.4 10$^{13}$~~ &   6 10$^{33}$ &   \\[3pt]
 J1814--1744~~ &  4.0   & 7.5 10$^{-13}$ & ~~5.5 10$^{13}$~~ &  $<$2 10$^{33}$&   \\[3pt]
 J1734--3333~~ &  1.2   & 2.3 10$^{-12}$ & ~~5.2 10$^{13}$~~ &  --  &   \\[3pt]
 J1846--0258~~ &  0.3   & 7.1 10$^{-12}$ & ~~4.9 10$^{13}$~~ &  4 10$^{34}$& in SNR Kes 75  \\[3pt]
 J1119--6127~~ &  0.4   & 4.0 10$^{-12}$ & ~~4.1 10$^{13}$~~ &  3 10$^{33}$ & in SNR G292.2--0.5  \\[3pt]
\hline
\noalign{\smallskip}\hline\noalign{\smallskip}
\end{tabular}
\end{table}

\section{Rotation-powered pulsars with high magnetic field}
\label{sec:rpp}

Table \ref{tab:highB} lists all the RPP for which the timing parameters
indicate a dipolar field
greater than,  or close to, the quantum critical value B$_c$ = $\frac{m^2 c^3}{\hbar
e}$=4.4 10$^{13}$ G.  Most of them
have spin periods of a few seconds and in the P-$\dot \textrm{P}$ diagram
lie in   the same region   of the AXPs/SGRs.
However, these four pulsars do not show magnetar characteristics,
such as the production of bursts/flares  or particularly strong and variable  emission.
In fact their X-ray luminosity is smaller than L$_{SD}$ \cite{piv00},
as expected for standard RPP.

The situation is quite different for the young PSR J1846--0258, located in the
supernova remnant Kes 75.
In this 0.3 s pulsar, the high inferred B derives from its large $\dot \textrm{P}$,  the highest of all RPPs.
Its spin-down luminosity of 8 10$^{36}$ erg s$^{-1}$ is certainly sufficient  to
power the  X-ray  emission  observed from the pulsar, L$_X$= 2.6 10$^{34}$ (d/6 kpc)$^2$ erg s$^{-1}$,
and  from the associated pulsar wind nebula,
L$_{PWN}$= 1.4 10$^{35}$ (d/6 kpc)$^2$ erg s$^{-1}$ \cite{ng08}.
Thus the discovery of   magnetar-like activity
from PSR J1846--0258 was quite surprising, because this pulsars was considered a
\textit{bona-fide} RPP, despite the lack of a radio detection
(easily explained as due to  an unfavorable orientation).
Four  short bursts from  PSR J1846--0258 were seen with RXTE  on May 31, 2006 \cite{gav08b}.
The onset of  the bursting activity  was coincident with a large
spin-up glitch and with the beginning of an enhancement
in the pulsed X-ray flux which lasted about one  month \cite{liv10}.
Spectral and flux variations associated with this event  were seen also
in the hard X-ray range \cite{kui09a}.

PSR J1119--6127  \cite{gon05,saf08}, has timing parameters very similar to those of PSR J1846--0258, and, like the latter,
it is located in a young SNR (G 292.2--0.5) and is surrounded by a pulsar wind nebula.
Its X-ray emission shows a large pulsed fraction and a spectrum consisting of a thermal-like
component with temperature kT$_{BB}\sim$0.2 keV and a hard power-law.
The strong thermal emission, quite remarkable for a young pulsar, might be the result of some recent magnetar-like activity, similar to that observed in PSR J1846--0258.
It would be important to monitor this pulsar, and other RPP with similar
parameters, to look for bursts.

\begin{figure}[b]
\includegraphics[scale=.4, angle=0]{kTLum.ps}
\caption{Temperature-Luminosity diagram for different classes of INS
(XDINS ($\times$), CCO ($\diamond$), persistent AXP (squares)).
For the sake of uniformity, all the temperatures are  from blackbody fits
(model atmospheres generally yield values smaller by a factor $\sim$2--3).
For the CCO the luminosities and temperatures of the two-blackbody fits are indicated.
For RCW 103 the figure shows two intensity states.
In AXPs the contribution to the luminosity from the power-law
spectral components are neglected.
The lines indicate the blackbody relations for emitting radii of 1 km and 10 km.}
\label{fig:HR}
\end{figure}

\section{Conclusions}

X-ray observations have revealed  an unexpected variety of manifestations of INS.
The thermal components in these objects span more than four orders of magnitude in luminosity (see Fig. \ref{fig:HR}), with a dependence on temperature broadly consistent with the blackbody relation L$\propto$T$^4$ for objects of similar size. The observed spectral and timing properties indicate that different energy sources, besides the obvious cooling of the primeval internal heat, are
at play in powering such thermal-like components.

The effects  of a strong and dynamic magnetic field are especially evident in the paroxysmic
behavior of SGRs and AXPs, but it is also possible that the variability
observed in the XDIN RX~J0720.4$-$3125\ be due to changes in the magnetic field configuration, and
certainly the magnetic field plays an important role in the thermal
evolution of INS. This is well demonstrated by the relatively high surface temperatures of
XDINS.  It is now clear that at least a fraction of the enigmatic CCOs are   young
pulsars with very weak magnetic fields.  However, also in these sources, localized regions of
stronger field might be present and be responsible for the small emitting areas inferred
from the spectral fits.

The magnetar-like behavior of  PSR J1846--0258,
together with the recent discovery of PSR J1622--4950 \cite{lev10},
the first ''radio-selected'' magnetar,  point to a closer  connection between
radio pulsars and magnetars than previously thought.
It is also noteworthy that SGR 0418+5729, one of the most recently
discovered SGRs, has an inferred (dipolar) field smaller than $\sim$3 10$^{13}$ G.

The data collected in the latest years, besides enlarging the sample of INS,
are leading to recognize objects with ''intermediate'' properties.
While this makes the boundaries between the different INS classes  less
well defined, it gives the prospects of unifying in a global coherent picture
all the varieties  of  INS manifestations.

High Energy Astrophysical Phenomena (astro-ph.HE)