Get 20M+ Full-Text Papers For Less Than $1.50/day. Start a 14-Day Trial for You or Your Team.

Learn More →

Is a quantum biosensing revolution approaching?

Is a quantum biosensing revolution approaching? Is biosensing revolution approaching? Review on Biocompatible ODMR Techniques 1;2 2 2 2 3 3 2 2;4 G. Petrini , E. Moreva , E. Bernardi , P. Traina , G. Tomagra , V. Carabelli ,I. P. Degiovanni , M. Genovese . Physics Department, University of Torino, Torino, Italy Istituto Nazionale di Ricerca Metrologica (INRIM), Strada delle cacce 91, Torino, Italy Department of Drug and Science Technology and NIS Inter-departmental Centre Torino, University of Torino, Italy Istituto Nazionale di Fisica Nucleare (INFN) Sez. Torino, Torino, Italy Understanding the human brain is one of the most signi cant challenges of the 21st century. As theoretical studies continue to improve the description of the complex mechanisms that regu- late biological processes, in parallel numerous experiments are conducted to enrich or verify these theoretical predictions also with the aim of extrapolating more accurate models. In the eld of mag- netometry, among the various sensors proposed for biological application, NV centers are emerging as a promising solution due to their perfect biocompatibility and the possibility of being positioned in close proximity to the cell membrane, thus allowing a nanometric spatial resolution down to the nano-scale. Still many issues must be overcome to obtain both spatial resolution and sensitivity capable of revealing the very weak electromagnetic elds generated by neurons, or other excitable cells, during their ring activity. However, over the last few years, signi cant improvements have been achieved in this direction, thanks to the use of innovative techniques. In this review, the new results regarding the application of NV centers will be analyzed and the main challenges that must be a orded for leading to practical applications will be discussed. PACS numbers: I. INTRODUCTION Electromagnetic eld sensing is of the utmost importance for several applications in current scienti c research, fostering the search for novel high sensitivity sensors. Several innovative electromagnetic eld sensors emerged in 1{10 the last years , whose main goal is revealing less and less intense elds with an increased spatial resolution. In particular high sensitivity sensing coupled to high resolution is of the utmost relevance in biological research, especially, for instance, in studies of human brain cell currents, which are typically extremely faint. Localized 11{14 monitoring of neuronal elds would allow not only the investigation of brain currents during cognitive processes in order to improve neurological diagnostic systems, but also identifying the early stages of neurodegenerative disease, 15,16 like Parkinson's, Alzheimers's disease and other forms of dementia . Furthermore, since localized temperature gradients and heat dissipation occur within cellular microdomains, the 17,18 exploitation of sensing probes for multi-task applications would be very fruitful . Among the various sensing 19,20 devices that have emerged over the years , promising sensors for the detection of biological elds are color 21,22 23{30 centers in diamond . Color centers are impurities in the crystalline matrix that, when stimulated, emit uorescence. In particular, the nitrogen-vacancy (NV) complex is by far the most promising due to its level structure . This dependence allows the realization of techniques for optical initialization and spin readout by means 33 34 of the Optically Detected Magnetic Resonance (ODMR) technique . Furthermore, its spin energy levels are 35{41 42 sensitive not only to electromagnetic elds , but also to temperature variations . These exceptional properties, together with their photostability at room-temperature and the non-toxicity of diamond , promotes the NV complex 44{49 as a very promising candidate for biological application . The paper is structured as follows: section 2 reports a brief description of the theory of quantum sensing, section 3 analyzes the magnetic eld generated by mammalian neuronal cells and cardiac tissue, section 4 deals with experiments aimed at the detection of cells elds and nally in section 5 the experimental techniques, used to enhance the sensitivity of NV centers to be used as biosensors, are highlighted. II. THE THEORY OF QUANTUM SENSING WITH NV CENTERS The nitrogen-vacancy (NV) defect is a natural complex of impurities in diamond crystalline matrix. This complex is composed of a substitution nitrogen atom and a vacancy-type defect, located in adjacent reticular sites . This system has a pyramidal symmetry (C ) and it has, as axis of symmetry, the line that connects the nitrogen atom 3v with the vacancy (see Figure 1a). With respect to the tetrahedral structure of the diamond, there are 4 possible arXiv:2006.00831v2 [physics.app-ph] 1 Jul 2020 2 orientations of this defect, all equiprobable in conditions of conventional syntheses. Moreover, there are two charged states in which it is possible to nd the nitrogen-vacancy defects, distinguished by the number of electrons involved. The 3 carbon atoms surrounding the vacancy contribute to sharing 1 electron each to the complex, while nitrogen contributes with 2. If, in total, only these 5 electrons are present in the system, the center is electrically neutral and it is referred to as NV , with total electronic spin S = 1/2. Alternatively, the defect can trap 1 additional electron from the surrounding lattice, creating the NV center. In this case the total electronic spin becomes S = 1, with spin component along symmetry axis fjm = 0 >;jm = +1 >;jm = 1 >g. The most promising con guration for s s s quantum sensing exploits the spin property of the NV complex. Its sp orbitals linearly combine to form 4 molecular orbitals: the lowest energy state of the ground con guration that is the orbital singlet, spin triplet state A , the 3 1 1 electronic excited states that are orbital doublet, spin triplet E, spin singlet orbital singlet E and A . a) b) FIG. 1: a) Diamond crystalline structure with nitrogen-vacancy defect; b) NV radiative state transitions that occur during 3 3 1 laser pumping. Radiative optical transition E ! A with 637 nm zero phonon line (ZPL), and non optical transition E ! 1 3 1 1 A with 1042 nm ZPL. Non radiative intersystem crossing (ISC) transitions subsist between E ! A and between E ! 1 1 A . As we will discuss, due to its electronic levels the NV complex is the most promising con guration for quantum sensing. By irradiating the complex with a 532 nm pump laser (Figure 1b), the electronic state is excited in a non-resonant way and afterwards it relaxes to the fundamental state emitting a photon with wave lenght between 637 nm (zero phonon line) and 800 nm (phonon sideband). While the optical excitation from the jm = 0 > state is spin preserving, the transition from jm = 1 > has a nite branching ratio into the metastable singlet E, with a lifetime of 300 ns. This singlet state relaxes into jm = 0 > through a non-radiative processes and weak infrared emission peaking at 1042 nm, leading to a drop in uorescence output up to 30% for a single NV , or 1-2% for a large NV ensemble, with respect to the situation when the system is initialized in jm = 0 >. This allows the optical readout of the spin state. Since the main focus of the present work is on sensing, from now on we will refer to NV as NV for simplicity. A. NV ground electronic state 3 34,50 The Hamiltonian of A , the ground spin state of the NV system, can be written in the following form : gs ^ ^ ^ ^ ^ ^ = SDS + SAI + IQI (1) ^ ^ ^ ^ ^ ^ ^ ^ where S = (S ; S ; S ) and I = (I ; I ; I ) are the dimensionless electron and nitrogen nuclear spin operators, x y z x y z respectively. The rst term represents the ne structure splitting due to the electronic spin-spin interaction, coupled by the ne structure tensor D. The second term is generated by the hyper ne interaction between NV electrons and 14 15 the nitrogen nucleus (I = 1 for a N nucleus, while I = 1/2 for a N nucleus), with the hyper ne tensor A. Finally, the third term represents the nuclear electric quadrupole interaction, with the electric quadrupole tensor Q. It should be noted that, in this notation, the component z coincides with the NV axis of symmetry. Due to the symmetry of 51,52 the NV center, D, A, and Q are diagonal in the NV coordinate system and, in terms of the natural spin-triplet 3 basis fjm = 0 >;jm = +1 >;jm = 1 >g, the Hamiltonian can be written as: s s s gs 2 2 == ? 2 2 ^ ^ ^ ^ ^ ^ ^ ^ ^ ^ = D [S S =3] + A S I + A [S I + S I ] + Q [I I =3] (2) gs z z x x y y gs z gs gs z h | {z } | {z } | {z } electronic spinspin interaction nuclear spinspin interaction electronnucleus spin interaction == where D ' 2.87 GHz is the zero eld splitting, Q is the nuclear electric quadrupole parameter, A and A are gs gs gs gs 53,54 the axial and non-axial magnetic hyper ne parameters . The parameters values are reported in Table I. Hyper ne parameters Value Zero eld splitting D ' 2.87 GHz gs == Axial hyper ne term A ' -2.14 MHz gs; N == A ' 3.03 MHz gs; N Transverse hyper ne term A ' -2.70 MHz gs; N A ' 3.65 MHz gs; N Nuclear electric quadrupole term Q ' -5 MHz gs TABLE I: Hyper ne parameters for the NV defect determined at room temperature. B. The Optically Detected Magnetic Resonance technique One of the characteristics that makes NV centers so attractive and convenient as key element in various type of sensors is the possibility to discriminate the spin components of the electronic state. This is allowed by the di erent coupling of the jm = 0 > state with a metastable level, compared to the jm = 1 > states and results s s in a variation of the photoluminescence (PL) of the defect under laser non resonant excitation. Optically Detected Magnetic Resonance (ODMR) consists in the application of a microwave eld (MW) on the sample, simultaneously with its exposure to a non-resonant laser at a frequency higher than the resonant one, corresponding to the energy gap between the ground and the E level (e.g. 532 nm) (see Figure 2a). When the frequency of the MW reaches the ground state resonance D of the NVs, with a certain probability (depending on the MW power), those NV centers gs will be initialized in the states jm = 1 > rather than jm = 0 >. As mentioned, this corresponds to a reduction in s s photoluminescence of the NV centers, as it can be observed, e.g., in Figure 2b, where a typical ODMR spectrum is reported with the expected uorescence dip at the zero eld splitting frequency D . gs a) b) FIG. 2: a) NV radiative state transitions that occur during laser pumping and MW excitation. The coupling of the state jm = 1 > with the metastable level generates a statistically lower uorescent emission than the jm = 0 > initialized s s electronic state. b) Fluorescence collected from the NV center as a function of the MW frequency varies. A dip in correspondence of the zero eld splitting D (resonance frequency of the undisturbed NV center, at room temperature) can be observed. gs The coupling terms of a NV center with the electric, magnetic elds and local temperature variations will be analyzed in following. 4 C. Magnetic eld sensing A static magnetic eld produces the well-known Zeeman e ect , that it is described by: == V  g  g gs B gs B gs N N ^ ^ ^ ^~ = S B + (S B + S B ) + IB (3) z z x x y y h h h h | {z } | {z } Zeeman interaction nuclear Zeeman interaction == where  is the Bohr magneton,  is the nuclear magneton, g and g are the components of the ground state B N gs gs electronic g-factor tensor and g is the isotropic nuclear g-factor. In the presence of relatively weak magnetic elds, it is possible to approximate the almost diagonal g-factor tensor in a diagonal form, with constant g = 2.003 . As reported in Table II, the interaction of the magnetic eld with the nucleus is 2000 times smaller and, consequently, it is typically neglected . The presence of external elds eliminates the energy degeneracy of the levels jm = 1 >, B e whose splitting become B , where = (see Figure 3). e z e a) b) 3 14 15 FIG. 3: NV ground-state A scheme. Above: a) N hyper ne states and b) N hyper ne states. Below: schematic ODMR spectra. The spectra are shown considering Zeeman splitting and hyper ne splitting. If, instead of a single NV center, an ensemble of NV centers is considered, up to eight magnetic resonance dips can be observed, due to the four possible orientation of the NV axis in the diamond's crystalline matrix (see Figure 4). For certain directions of the magnetic eld, some resonances can be degenerate. A NV-based magnetometer can be realized, for example, by applying a bias eld along the NV axis, removing the degeneracy, so that changes in the magnetic eld projection along this axis a ect the resonance frequencies almost linearly. Another option is to use all four NV alignments; although the eight ODMR frequencies have more complicated dependence on B, this option yields information about the direction of magnetic eld . 38,39 35,40 The use of NV center as a magnetic eld sensor rstly was proposed in and demonstrated with single NV and NV ensembles in 2008. D. Electric eld sensing The Hamiltonian describing the interaction with the electric eld was derived from molecular orbit theory by Doherty et al. and it can be written in the following form: ^ ^ V S gs == 2 ? 2 2 ? ^ ^ ^ ^ ^ ^ ^ = d (E + F )[S ] + d (E + F )(S S ) + d (E + F )(S S + S S ) (4) z z x x y y x y y x gs z gs y x gs h 3 5 B≠0 B=0 |0> |-1> |+1> |±1> b) MW [GHz] a) MW [GHz] B≠0 c) MW [GHz] FIG. 4: ODMR spectra in the absence a) and in the presence b) of an external bias magnetic eld. The magnetic eld lifts the degeneracy of the jm = 1 > states and results in two separate dips in the ODMR spectrum. c) An example ODMR spectrum (excited at 532 nm) with a magnetic eld in an arbitrary direction for an ensemble NV centers in diamond. Each of the four NV alignments has a di erent magnetic eld projection along its quantization axis, leading to eight ODMR peaks (two for each NV alignment). For each dip a coupling with the nuclear spin of the N atom generates additional three hyper ne levels . == where d and d are respectively the axial and non-axial Stark shift components of the permanent electric dipole gs gs ? 57 ~ ~ moment d in the ground triplet state , E is the electric eld and F is the mechanical strain. gs 41,58 ~ ~ According to Equation (4) the e ect of the electric eld E plays the same role as mechanical strain F . The strain depends on the diamond material: in single-crystal samples, the mechanical strain eld is substantially negligible; while, in polycrystalline ones, a relatively high strain eld is induced by the growth conditions, leading to a splitting of the spin state jm = 1 > even in absence of external elds. The frequency shift caused by the electric eld is much smaller than the shift produced by the presence of a magnetic eld (see Table II). For this reason, in order to reliably measure this second-order e ect caused by the Stark shift, it is necessary to decouple it from the Zeeman shift. Brie y, the ne structure Hamiltonian of the NV ground state, describing the energy levels of the electronic spin ^ ~ ~ ~ states due to the spin (S) interaction with the static magnetic (B), electric (E), and strain (F ) elds, can be written in terms of the natural spin-triplet basis fjm = 0 > , jm = +1 >, jm = 1 >g in the following matrix form: s s s 0 1 B iB B +iB x y x y p p 0  g  g B e B e 2 2 B C B +iB x y ? ^ B C H =  g hD +  g B hd (P iP ) gs B e B e z x y gs @ 2 A B iB x y g hd (P + iP ) hD  g B B e x y B e z gs where it is possible to observe that the natural-spin basis vectors are eigenstates of the Hamiltonian only in the == presence of both the magnetic and electric eld aligned with the NV axis. In this condition, D = D + d P gs gs z describes the frequency shift of the resonance lines resulting from the zero- eld splitting and from the Stark e ect ~ ~ ~ associated with the component of the vector P = E +F . Otherwise, external elds not aligned to NV symmetry axis produce a non-diagonal matrix, and therefore energy levels of unde ned spin. In particular, the presence of additional transverse strain and electric- eld components P modi es the ground-state structure. The Hamiltonian assumes a quasidiagonal form considering a new spin basis fj0 >, j+ >, j >g, obtained by a eld-dependent mixing of the jm = +1 > and jm = 1 > spin states according to the following unitary operator: s s 0 1 1 0 0 B   C ^ i  i 2 2 U = 0 e sin( ) e sin( ) @ A 2 2 i i 2 2 0 e cos( ) e sin( ) 2 2 Fluorescence [a.u] Fluorescence [a.u] Fluorescence [a.u] 6 ? ? where tan() = P =P and tan() = (d P )=( g B ) are the eld-dependent phases de ning the spin state mixing. x y B e z gs The Hamiltonian takes the following form in the fj0 >,j+ >,j >g basis: 0 1 ? ? 0 c  g B c  g B 1 B e 2 B e B C ^0 ^ ^ y  ? H = UH U = gs gs @c  g B hD + W 0 A B e c  g B 0 hD W B e with ? ? 2 2 W = (hd P ) + ( g B ) (5) B e z gs The complex constants c and c represent the phase of the matrix elements and B is the transverse component 1 2 of the magnetic eld with respect to the NV axis. If B  0, the non-diagonal terms can be neglected and the Hamiltonian can be regarded as diagonal in the basis fj0 >, j+ >, j >g. The energy di erence between the j0 > and the j > states is hD  W , corresponding to ODMR resonances separated by 2W/h depending on the strengths of the magnetic, electric, and strain elds, as well as their orientations with respect to the axes of the NV center. Since the states j > are a coherent superposition of the states jm = 1 >, we underline that the ODMR resonance is observed also in this case as a reduction in the uorescence emission at the new MW resonance frequencies. E. Temperature sensing Another interesting feature of the NV complex is the temperature dependence of its spin levels . Indeed, the microscopic origin of D , also called the zero eld splitting (ZFS) parameter, is due to spin-spin interactions in gs the NV's orbital structures, and the value depends on the lattice length, which is strongly correlated to the local temperature. When the local temperature increases the diamond lattice spacing of the NV center increases as well, lowering the spin-spin interaction and reducing the ZFS parameter D . Under ambient conditions D ' 2:87 GHz gs gs and the temperature dependence is dD=dT ' 74kHz=K . In general, the ZFS parameter shows a non-linear dependence, and its value increases when temperature decreases . Property Coupling coecient e 1 Magnetic eld = ' 28 GHz T N N = ' 15 MHz T Electric eld d ' 3.5 mHz V m ==;gs d ' 0.17 Hz V m ?;gs Temperature @D =@T ' -74 kHz K gs TABLE II: Coupling coecient of the NV center with the external elds and temperature. To realize a NV-based temperature sensor, the most obvious solution is exploiting the D temperature dependence. gs ~ ~ ~ This requires that no external eld is present (B; E; F = 0), i.e. jm = 1 > is degenerate. In this case, an increase in temperature leads to a decrease in the resonance frequency, associated with a shift of the degenerate levels jm = 1 > towards the level jm = 0 >. Nevertheless, this is the simplest but not the optimal solution, since, even in the absence of applied elds, the sample may have an internal strain and may be a ected by the Earth's magnetic eld. Unless it is possible to nd a diamond sample with negligible F and to design an experimental set up able to reasonably compensate for the external magnetic eld (e.g. Helmholtz coils), the dips would not be perfectly overlapped because of the non perfect degeneracy of jm = 1 >, thus showing a larger full-width-at-half-maximum (FWHM) and therefore a lower resolution. A better solution is to apply an external magnetic eld in order to signi cantly separate the spin levels. However in this con guration, a single dip can shift for a temperature variation, but also for a variation of magnetic eld. To decouple the two contributions it is sucient to monitor both jm = +1 > and jm = 1 > spin states at the same s s time, using simultaneous driving of the microwaves in ODMR technique . As it can be seen in the Figure 5, by simultaneously monitoring the initial dips (red curves), it is in principle possible to understand if there are variations in the magnetic eld (the dips move in opposite directions) or in temperature (the dips move in the same direction). In this perspective, an improved technique exploiting the application of an intermediate transverse bias magnetic ? 61 ? eld B has been implemented . Similarly to the case just discussed, the application of B removes the degeneration 7 FIG. 5: Example of magnetic and thermal shifts of the spin resonance, in ODMR spectra. Dips with equal colors correspond to paired resonances. The colors represent the timeline of the dips. The initial dips is red, then green and nally blue. of jm = 1 > and therefore improves the FWHM. The intensity and the transverse direction of that eld creates a quantum superposition of states, which is insensitive to magnetic elds but sensitive to temperature . In this con guration, the expectation value of the spin along any direction is small, implying the degeneracy of the hyper ne structure between the levels jm = 1 > (except for the quadrupole contribution Q , which seperates jm = 0 > I gs I from jm = 1 >). In Figure 6 the corresponding scheme of the spin energy levels are reported (only the N isotope is considered as it is the most common). In this situation, the ODMR spectrum reduces to two dips (instead of 6), providing a substantial improvement in the signal-to-noise ratio. This particular orientation of the magnetic eld ensures the protection of the measurements from the noise of other possible magnetic elds. In fact, the NV spin is non-sensitive to the magnetic eld uctuation, because the contribution of the magnetic component only appears at the second order in the Hamiltonian. 3 ? FIG. 6: NV ground-state A scheme, in presence of intense transverse magnetic eld B . III. BIO-SENSING As mentioned above, NV sensors are particuarly suitable for biological sensing. Before describing the experiments focusing on the NV-based sensor it is necessary to specify the type of biological specimens of interest, the expected magnitude of the electromagnetic eld produced by these specimens and the principal parameters such as sensitivity and spatial/temporal resolution, required for the NV-based sensors. This section, after reviewing some of the devices typically used for biosensing, analyzes in detail neuronal and cardiac cells. Higher sensitivity and resolution of electromagnetic elds is considered necessary to expand the understanding of the fundamental processes regulating the interaction of these cells. A. From the conventional electrophysiological techniques to NV sensors The electrical activity of excitable cells can be investigated by means of the conventional patch-clamp technique or through Micro Electrode Array (MEA) recordings. Single-cell recordings, performed under voltage- or current-clamp con guration, respectively allow to monitor ion currents or the membrane resting potential, postsynaptic responses and action potential ring activity . Besides having an extreme versatility (monitoring overall electrical events from 8 the whole cell, from microdomains of the cell membrane or even from single channel proteins), patch-clamp has a high temporal resolution and high sensitivity: all these features make this electrophysiological approach the gold standard for measuring the electrical activity. Though, patch-clamp is rather invasive, as it damages the cell membrane through the recording electrode: this implies that only one recording is feasible for each cell . On the contrary, MEA is a non-invasive approach, used to measure the membrane potential variations from many cells simultaneously. The MEA probe is structured as an array of sensing electrodes, of variable geometry dimensions and material, which are im- mersed in a glassy (insulating) double layer. Commonly, sensing electrodes are made of titanium or indium tin oxide 2,3 (ITO), and have a diameter that can vary between 10 and 30 m . By means of MEAs, it is possible to monitor the electrical activity of a neuronal network as a whole, and measuring its changes along with its maturation, even though informations on the biophysical properties of ion channel cannot be directly inferred. This speci c measurement need, requiring non-invasive and iterative detection for biological applications has prompted the study and realization of di erent devices. In the following we discuss and compare the most promising ones. Promising devices for the detection of weak magnetic elds, in addition to the NV-based sensors, are the super- 4{6 7 conducting quantum interference device (SQUID) sensors and chip-scale atomic magnetometers (CSAMs) . Until now, the measurement of very weak magnetic elds was the domain of SQUIDs sensors. These sensors have reached 1=2 2 5 sensitivity levels of (0.9{1.4) fT/Hz with a pick-up coil area of the order of 1 cm . However, SQUIDs require cryogenic cooling, which, in addition to implying signi cant cost and maintenance complexity, requires positioning the sensor a few centimeters from the sample. An alternative is o ered by the CSMAs, that are based on microfabricated alkali vapor cells integrated with small optical components such as diode lasers and ber optics. These devices have 1=2 37 reached sensitivities below 5 fT/Hz at sensor volume 8 mm . However, despite the exceptional sensitivity, the minimum working distance between sensor and magnetic source for CSAM or SQUIDS remains at least few mm, that makes them unsuitable for monitoring individual cell signals or small tissues, being the amplitude of the magnetic eld decreasing quadratically with distance. The sensors for the detection of electric elds, emerging in the last few decades are single-electron transistors (SETs) , that are a promising candidate for achieving higher detection sensitivity due to the Coulomb oscillations. However, the existence of a SET- based biosensor has emerged only in recent years , probably because of their diculty of the room-temperature operation. Finally, in recent years there has been a growing interest in the use of temperature sensors capable of operating on a nanometric scale. It has been known that local temperature variations at the intracellular level play a fundamental role in cellular activities related to body temperature homeostasis and energy balance . Particular attention is paid to the possibility of measuring local temperature variations of cell organelles (i.e. nucleus, mitochondria, etc.) or 67,68 ion channels. For example, di erent simulation model shows a hypothetical variation in temperature at the level of the ion channels, due to the ow of the ions from the inside to the outside of the plasma membrane, during the genesis of the action potential. Due to the diculty of this local measurement, no one has ever measured this thermal variation. Interestingly, temperature changes may drastically alter the neuronal ring frequency, as demonstrated by Guatteo et al. . Currently uorescence probes are powerful method used to study intracellular temperature variation thanks their high spatio-temporal resolution. The probes typically used for this measurement are organic or inorganic 70{73 9,74 uorescent probes, such as uorescent proteins, organic dyes , quantum dots (QDs) and many others. Organic proteins are biocompatible probes, rather stable and very easy to chemically target. But there are di erent problems related the use of these probes: these are often auto uorescent and to avoid the phenomenon it is necessary to add speci c quenchers; they cannot be used for a long time, in fact these sensors su er from photobleaching and unstable photoluminescence. In the best case scenario, the probe degradation consists of uorescence suppression, in the worst case scenario it releases an electron that binds to nearby molecules making them toxic. These probes are organic and by their nature they are also subject to even weak pH variations, for this reason it is fundamental a strict control 75,76 of the cell environment . The inorganic probes such as quantum dots (QDs) have the advantage of being stable in uorescence, have a high sensitivity to temperature variations and their nanometric size allows obtaining a spatial resolution useful for cellular measurements. Although the size of these sensors would allow spatial resolution limited by the di raction limit only, their chemical composition is found to be non-biocompatible in most of the cases. Other 10,66 temperature sensors are based on up converting nanoparticles (UCNPs) : nanoscale particles (diameter 1-100 nm) that exhibit photon upconversion, i.e. when stimulated by incident photons they are able to emit uorescence's of shorter wavelength. They are usually composed of rare-earth based lanthanide or actinide-doped transition metals. Their core-shell structure allows sensor compatibility, however, sensitivity is not high. Extremely interesting devices able to realize all these measurements (magnetic, electrical and temperature sensing) eventually at the same time, are one based on the NV center in diamond. The advantages of these sensors are man- ifold: they have stable photoluminescence in the visible and near-infrared range, their chemical composition ensures resistance to photobleaching and diamond is an inert and therefore biocompatible material , so cell/neurons can be 27,28,77 grown directly on its surface or nanodiamonds can be injected inside them, allowing for sub-cellular spatial resolutions with a non-invasive techniques. Finally, NV sensors can operate at room temperature and, in more de- 9 78 79 tail, their dynamical range of temperature sensing extends further 500 K for both bulk and nanoscale diamonds. In the remainder of this section, the modelization of the target neural and cardiac signals for NV-based biosensing will be reviewed. B. NV center as sensor for neuronal signals In the last decades, neuroscience has attracted great interest beyond the scienti c community. Because of the increase in life expectation, cases of neurodegenerative diseases such as Parkinson's, Alzheimer's, Huntington's disease and many others are constantly growing. Currently, these diseases are incurable, even symptoms mitigation is dicult because of late diagnosis when most of the neurons involved have been irreparably damaged. This reason strongly prompts to develop new increasingly precise and sensitive techniques, allowing a deeper understanding of neuronal circuits ranging from functioning of synaptic sites to the behavior of the entire neuronal network. Neurons are the functional units of the nervous system. They communicate via electrical signals, known as action potentials. The action potential (AP) consists in the variation in time of the membrane potential V , where V = e e m m in out is the electrical potential di erence between intra- and extra-cellular side of the cell membrane. The AP characteristic + + 2+ waveform is shown in Figure 7b. The AP pulse is caused by several ionic species (Na , K , Ca ), which ow through the neuronal membrane. a) b) FIG. 7: a) Single neuron simpli ed sketch. In the upper box a zoom of the neuronal membrane is reported, where the ionic current and the corresponding magnetic eld are schematized. In the lower box the axial current and the relative magnetic eld are shown. b) Schematic representation of neuronal action potential (AP). Resting membrane potential (-V ) is -70 mV. When V is driven and exceeds the threshold (following an initial stimulus), a rapid membrane depolarization occurs. In this phase the Na channels open, allowing sodium to enter in the neuron and bringing V to approximately +35 mV. Then the + + repolarization phase begins, caused by Na channels inactivation and opening of K channels. This outward current drives + + the membrane potential close to -93 mV (hyperpolarization). Finally, the Na /K ATP-ase restores the initial conditions. During the depolarization, the in ux of positive charges produces local internal and external longitudinal currents, which are responsible for the AP propagation in the axon adjacent area. The propagation directionality is guaranteed by the AP refractory period: although the local currents propagate in both directions, a new AP cannot be triggered in refractory membrane area. The two electrophysiological techniques mostly used to study cell excitability and synaptic transmission in a neuronal network are the patch-clamp and the MEA. In the last decade scientists have tried to study more and more speci cally the propagation of the electrical signal from the cell body (or soma) to the whole dendritic tree. In other words, the goal would be to create a device that allows scanning the neuron point by point from the soma to the axon and the dendrites, following and characterizing the electrophysiological variations of the electrical signal during its propagation. The technology closer to this ambitious goal is the one of the CMOS-MEA, that allows having a much higher density of electrodes with respect to the traditional MEA technology. Numerous stud- 80? {82 ies have managed to scan the path of the electrical signal in a neuronal network at the level of the single neuron . Bakkum et al. recently have developed a high electrode density CMOS-MEA device capable of stimulating a speci c area and simultaneously scanning the signal along some points from the soma to the axon. Clearly, this technique is much more sensitive than MEA, but given the stochasticity of the cell's placement in space, it requires cells to be marked in order to follow their path. Recently several groups have correlated this technology to the technique of optogenetics. They tagged the genes of interest and activated them following an optical stimulation and 10 83,84 simultaneously followed the signal thanks to the integration of the CMOS-MEA . However, these techniques do not allow following the entire dynamics of the action potential, but to have a scan of a region depending on the position of the electrodes with respect to the neuron with its axon and its dendritic body. NV sensors may therefore have a huge impact on these applications: nanodiamonds can be targeted on the membrane surface or, alternatively, cells can be plated and cultured on a bulk diamond . Indeed, taking advantage of diamond biocompatibility and the exceptional spatial resolution displayed by color centers in diamonds , it will be worth exploiting these properties for a timely reconstruction of the AP dynamics. Furthermore, the possibility of positioning them adjacent to the cell membrane has the advantage of experiencing stronger magnetic elds. However, since neuronal magnetic elds are extremely weak (' pT ), their detection appears to be challenging even for NV-based sensors, at least for mammalian cells, while measurements have been performed on giant neurons of invertebrates . To predict the electromagnetic elds intensity created by the AP, and therefore to understand what sensitivity of the NV sensors is needed to sense it, it is necessary to model how the AP develops and propagates. 11{14 Hodgkin-Huxley model allows estimating the ionic current owing through the neuron membrane (when the ion channels are open). For the human neuron, the total estimated ionic current, sum of the single channels contribution I is: ion I = I ' 2 pA m ? ion and the current pulse typically lasts t ' 1 ms. Each I generates a magnetic eld (see Figure 7a), which can ion ~ ~ be estimated by means of the Biot-Savart law: B  l =  I . However, the resulting amplitude of these elds ion 0 ion depends on the channels density, which largely varies depending on the axon area being considered. Furthermore, we tot note that the B eld, sum of the contributions of the eld produced by the various channels, can be vanishingly small ion 86,87 on average, because of the di erent elds directions. To this purpose, channel clustering may be very signi cant . Assuming a current of 100 pA=m and considering that the NV sensor positioned at an average distance of few 88,89 nanometers (by selective targeting the channel using functionalized NDs ), a magnetic eld of about 0:1 5 nT (or even higher) could probably be sensed. This hypothesis is now under experimental analysis . Current owing through the membrane is not limited to the charge ow through ion channels, as longitudinal currents, but one should also consider the ow along the neuron axis, that is responsible for the AP propagation. These currents also generate a magnetic eld, around the neuron (see Figure 7a). Both the axial current and 90{92 90 the corresponding magnetic eld have been estimated . In particular, Ref. goes beyond the simpli cation of the Hodgkin-Huxley model, introducing the spatial and temporal progression of the AP along the various neuronal compartments, into which they have divided the axon. The theoretical prediction is a maximum eld B ' 3 pT axial on the external membrane near the Ranvier node and a eld B ' 2.3 pT on the myelin sheath external surface axial in those regions where the axon is wrapped by it. The maximum magnetic eld was also calculated by Isakovic et al. in Ref. for the nerve composed of 100 axons, obtaining only B ' 6 pT. This is due to the cancellation of the magnetic eld component, caused by di erent axial axons within the same nerve, bringing opposite directional currents. This estimated magnetic elds, in reality, are compatible with the elds detected by magnetoencephalography (MEG). MEG is able to detect elds if the order of 15 93 10 T because of the distance from the source . Considering these values, a NV sensor positioned on the neuron surface or a few micrometers from it, should have a temporal resolution of about 0.1 ms (in order to be able to trace the time variation), and spatial resolution of about 10 m (which would allows a good reconstruction of the AP propagation, being the axon length ranging from 0.1 m to 1 m). Thus, the NV sensor should have a minimum sensitivity of : p p 1=2 = B t ' 3 pT 0:1 ms ' 30 fT Hz (6) min The NV sensor optimal sensitivity is in principle limited by the quantum projection noise. This fundamental sensitivity limit for spin-based magnetometers is given by : 1 1 = (7) nT Where is the magnetic coupling coecient (Table I), n represents the number of NV centers and T their characteristic dephasing time. It is important to underline that the number of NV centers n refers to the sensing volume. As mentioned, for the single PA detection the sensing volume should be around 10 m , the size of the cell. 77 6 3 In the Ref. , the estimation of the parameters n ' 3 10 cm and T ' 450 ns determines a spin projection noise 1=2 3 value of  ' 30 pT Hz for the sensing volume of 10 m (the experimental sensitivity reached is instead  ' 15 q 11 1=2 6 3 pT Hz for the sensing volume of 5  10 m ). This value is still 1000 times larger than the sensitivity required for the detection of a single AP. However, as will be discussed in section 5, it is possible to optimize both the above mentioned parameters to improve the performances. Once the biomagnetic eld B (x~ ; t) has been measured, to reconstruct the unknown currents generating it, one should solve an inversion problem. In general, its solution is not unique, due to the existence of the so-called "magnetically silent" currents (i.e. the ones producing magnetic elds that almost cancel each others) and due to 96,97 the fact that the magnetic eld can be in uenced by the electric eld . However, in the single axon case, it can be uniquely resolved. On the contrary, in the biological tissue case and in the 3D structures case, that cannot be traced back to standard models (such as a spherically symmetrical conductor or a horizontally layered medium), the solution is not unique. In some cases this is resolved by the knowledge of the electric eld on the conductor surface . C. NV center as sensor for cardiac signals The human (and animal) heart generates the body's most intense electromagnetic eld. In particular, by comparing measurements performed externally to the human body, the electric eld generated by the heart, measured through the electrocardiogram (ECG) is about 60 times stronger than that the one of the brain, recorded by an electroen- cephalogram (EEG). In addition, the heart magnetic eld detected by the magnetocardiogram (MCG) is about 5000 times higher than the neuronal magnetic eld detected by magnetoencephalography (MEG): 0,05 nT (heart) vs 1 fT (neuron). Thus, ODMR based on NV sensors can also nd very signi cant applications in studying cardiac cells and tissues. To achieve a rst qualitative estimation of the magnitude of the magnetic eld in this case, one can start from a very simpli ed model: the spherical heart . Although this model is not physiologically accurate, it allows to extrapolate analytical solutions. In a more recent work , is proposed a further assumption concerning the origin of the currents. There are two currents sources in the heart: the rst consists of intracellular currents, the second is given by the anisotropy of the tissue . Regarding the rst current contribution, the authors consider a spherical shell of cardiac tissue, which covers a blood cavity and is surrounded by an external bath of unlimited electrical conduction. The heart bers propagate in the z direction and a variation of the membrane potential V is assumed following the activation of the action potential (AP), started at  = 90 (see Figure 8a). In this work the electric eld is evaluated using the bidomain model and considering a situation of quasi-stationarity (although V depends on time due to the action potential propagation, it is assumed that, given a certain V (t ), m m 0 one can derive current and magnetic eld in a quasistatic way). 101,102 Thus, the electric potential is obtained, using the continuity equations and the boundary conditions , the current density distribution is obtained using Ohm's law and nally the magnetic eld using Biot-Savart's law. Considering 103 98 the anisotropic electrical conductance data , the V values and typical heart dimensions , it turns out that the magnetic eld is stronger near the internal and external surfaces tissue while it is weaker in the heart wall. The peak value of the magnetic eld is around 14 nT (see Figure 8b). At the heart center, instead, the magnetic eld reduces to B = 2 nT . This is due to the fact that intracellular and extracellular currents are in opposite directions with almost the same magnitudes in the depths of the tissue and, therefore, the corresponding magnetic elds essentially cancel each other. Considering a planar cardiac tissue sample, the spherical shell method is no longer valid. In this last case it has been found that the magnetic eld reaches a peak value B = 1 nT . The heart AP is about t = 300500 ms long, however for some cardiac cells, such as ventricular or rapid response cells, the AP rapid rise occurs in 1 ms, as in the neuronal case. Considering a human heart, a NV sensor positioned on the heart surface should be sensitive to magnetic eld B = 14 nT, with a temporal resolution of about 0.1 ms (in order to be able to trace the time variation t even in the case of the AP rapid rise), and a spatial resolution of about 10 m (which would allow a good reconstruction of the PA propagation, being the heart radius of about 40 mm ). This corresponds to a minimum sensitivity: p p 1=2 = B t ' 14 nT 0:1 ms ' 140 pT Hz (8) min This value can be considered a useful intermediate step for the application of the actual NV-based biosensing tech- nologies, with the aim of reaching sensitivity that allows the detection on neuronal signals. 12 a) b) FIG. 8: a) The model of a spherical heart, taken from the reference . Part of the spherical shell has been cut out to show the heart wall. The black curves indicate the ber orientation. The pink tissue has a transmembrane potential of 20 mV, and the blue tissue has a transmembrane potential of 80 mV. The green curve shows the magnetic eld. The endocardial (inner) surface has radius r , and the epicardial (outer) surface has radius r . b) The magnetic eld over a cross section of the heart. 1 2 The dashed curves indicate the heart inner and outer surfaces. An area 40 mm by 40 mm is shown. IV. METHODS AND BIO-APPLICATIONS To exploit NV-centers for biosensing, it is necessary to set up an optical microscope equipped with an ODMR apparatus featuring a non-resonant laser and a microwave antenna positioned near the diamond sample designed for a microwave source operating in the 2-4 GHz range, as described in section 2.2. The next sections are devoted to the presentation of several biosensing experiments exploiting NV centers in diamonds. Some of them are "proof-of-principle" in vitro tests on a cell culture, while others are experiments carried on living organisms (in vivo experiments). The reviewed experiments will be divided in two categories describing respectively applications exploiting bulk diamond and experiments relying on nanodiamonds. A. Bulk diamond applications As discussed in the previous section, the electromagnetic elds produced by the excitable cells (as neurons, neuronal-like chroman cells, heart cells) even in mammals, are typically extremely weak (pT ). For this reason, several "proof-of-principle" experiments addressed measurements of elds produced by cells with peculiar electro- magnetic properties. 106{109 III II Among the most suitable ones there are magnetotactic bacteria (MTB) containing magnetite ( Fe FeO ) III II or ferrite ( Fe FeS ) bacteria magnetic particles (BMP). The nanometer size of the BMPs is small to generate a single magnetic domain, but sucient to create a permanent magnetic moment ~ . This produces a cell magnetic BMP moment ~ = ~ , given by the sum of the BMP individual dipoles, which is exploited by the MTB to MTB BMP 110,111 orient itself with respect to the earth's magnetic eld . Among the various uses in the biomedical eld, Sage et al. used Magnetospirillum magneticum AMB-1 for bio- magnetic imaging. The bacteria used in this work create magnetic nanoparticles with cubo-octahedral morphology and an average diameter of 50 nm. The experiment was performed both with bacteria dried on the surface of diamond chip implanted with NV centers, as well as with bacteria stored in phosphate-bu ered saline (PBS) and laid on the chip surface (in vitro experiment). The diamond sensor used to perform this experiment was a high-purity single-crystal diamond chip, with a 10 nm 17 3 layer thickness of NV centers. The estimated surface density of nitrogen-vacancy centers was 3  10 cm in the 18 3 case of experiment on bacteria in the liquid medium and 10 cm for dry bacteria. In the case of dry bacteria the objective was to demonstrate the possibility of measuring their static magnetic eld, exploiting ODMR measurements at di erent bias magnetic eld orientations (B = 3.7 mT). bias In the case of live bacteria in the liquid medium, it was shown that it is possible to evaluate the magnetic eld generated by the bacteria dipole ~ along the [111] crystallographic axis of the diamond, when also the bias MTB magnetic eld is oriented along it. Furthermore, cell viability was assessed immediately after magnetic imaging (lasting 4 minutes), using a standard uorescence-based "live-dead" assay obtaining a viability of about 44%. Cells 13 mortality was attributed to the laser heating, since preliminary tests showed that 1 hour exposure to microwaves did not cause substantial cells mortality. Cells vitality was however partially preserved thanks to the strategy used to decouple laser light from the biological sample. Indeed in this set-up the laser impinged on diamond at an angle greater than the critical angle for the diamond{water interface, resulting in its total internal re ection within the diamond. A wide eld optical microscope was used for both MTB samples, with a eld of view of 100  30 m of the sample surface and a resolution of 400 nm. A CMOS camera was used to image the single magnetic nanoparticles inside the MTB. Their magnetic eld was of the order of mT. Thanks to these measurements, the total magnetic moment ~ was determined by numerically tting the modeled eld distribution to the measured ones, with a mean value MTB 17 2 of 5 10 m A. The magnetic eld estimated from the ODMR measurements was compared with a scanning electron microscope (SEM) measurements. The position of nanoparticles revealed by the SEM was used to model the magnetic eld 108,109 they generated . The two measurements were in excellent agreement and their values were compatible with 112,113 the data reported in . This highlights the potential of NV centers, able to perform sub-cellular magnetic eld measurements at room temperature, allowing real-time imaging of magnetic dipole creation, single MTBs chain 111 114,115 dynamics and magnetic particles formation in various organisms . Another "proof-of-principle" test was carried out by Davis et al. with the aim of measuring the magnetic eld generated by iron oxide nanoparticles (IONs) incorporated in murine RAW 264,7 macrophages (a line estab- lished from a tumor induced by Abelson murine leukemia virus and often studied in relation to immune responses). The cells, after having phagocytized the iron ions (about 200 nm in size), were dried on the surface of a bulk diamond. High resolution magnetic imaging was performed exploiting ODMR measurements with a succession of bias magnetic eld (B = 10 mT) for each orientation of NV complex, for a total duration of 2 hours. Projection bias eld maps were combined to form 3 orthogonal eld maps, from which the cellular magnetic moment was obtained. For this experiment, the central dip of the hyper ne transition was used, which allowed to achieve a sensitivity of 17 nT at 1 m in plane resolution, sucient to reveal these magnetic nanoparticles. To extend this technique to diagnostic imaging, Davis et al. performed NV magnetometry on liver specimens from a mouse model of hepatic iron overload, generated through intravenous administration of 900 nm IONs to C57bl/6 mice. To reduce the deposition of optical and thermal energy, the sample was illuminated for only 5 minutes with a duty cycle of 50% and the ODMR technique was carried out with the bias magnetic eld along only one of four NV axes. Furthermore, the laser beam was directed on the sample in total re ection mode. With these choices, time-lapse images of magnetic elds that coalesced within the macrophages after ION internalization were evaluated along that NV axes. This experiment highlights the possibility of study the spatial distribution of iron deposits in the liver and other tissues. This has been a topic of interest in clinical literature as an indicator of disease state and the magnetic resonance imaging is becoming increasingly important in non-invasive quanti cation of tissue iron, overcoming the drawbacks of traditional techniques (liver biopsy) . Barry et al. studied individual neurons of marine worms (Myxicola infundibulum ) and squids (Loligo pealeii ). The marine worm has a long axon , which stretches over its entire length (tens of mm and diameter of about 5 mm). The giant squid neuron (about 0.5 m long) did not extend over the entire length and is isolated following speci c protocols . An initial proof-of-principle test was performed on isolated neurons for both species. The AP is stimulated by means of a current pulse, received by an electrode directly in contact with the neuron. The pulse was generated by a current of about 10 mA, had a duration of about 1 ms and was repeated with a frequency of 0.4 Hz for the worm and 100 Hz for the squid. The AP generation and its propagation was veri ed by micro-electrodes (see Figure 9A). From this axonal AP 119{121 intracellular time trace, the shape of the associated magnetic eld could be modeled (see Figure 9B). This was compared with the experimentally measured magnetic eld, performed by means of the NV-based sensor in contact with the excised single neuron. Traces are shown in Figure 9C and 9D respectively for the worm and the squid neuron. These measurements were performed using the ODMR technique at bias magnetic eld B = 0.7 bias mT, oriented along two diamond axes and perpendicular to the axon axis (being the magnetic eld generated by the AP pulse perpendicular to this last one). In Ref. , Barry et al. carried on also a measurement on a living worm. The worm was directly xed on the diamond and the distance between the neuron and the active NV layer was about 1.2 mm (see Figure 10A). The magnetic eld generated by the propagation of the AP pulse measured by ODMR technique is shown in Figure 10B. It is smaller than the one measured in the excised neuron, but its value is compatible with the increasing sensor distance. The diamond sensor, exploited an electronic grade (N < 5 ppb) single crystal chip, with a NV center layer of 13 17 3 m. This layer had a NV centers density of d = 3 10 cm and a characteristic dephasing time T = 450 ns. The 6 3 11 sensing volume is V = 5 10 cm , consequently the number of potentially stimulated centers was n = 15 10 . 14 FIG. 9: Measured AP voltage and magnetic eld from excised single neurons, taken from the reference . meas A) Measured time trace of intracellular axonal AP voltage  (t) for giant axon from M. infundibulum (worm). in B) Calculated time trace of AP magnetic eld B(t) for M. infundibulum extracted from data in A. C) Measured time trace of AP magnetic eld B(t) for M. infundibulum giant axon with N = 600. avg D) Measured time trace of AP magnetic eld B(t) for L. pealeii (squid) giant axon with N = 375. Gray box indicates avg magnetic artifact from stimulation current. 1=2 Referring to the Equation (7), the fundamental sensitivity limit:  ' 10 fT Hz , while the sensitivity reached 1=2 experimentally was  ' 15 pT Hz , allowing, anyway, a reliable measure of the magnetic elds generated by these animal species (of the order of nT). Further development should be needed for revealing those of human neurons (of the order of pT). FIG. 10: Single-neuron AP magnetic sensing exterior to live intact organism, taken from the reference . A) Overhead view of intact living specimen of M. infundibulum (worm) on top of NV diamond sensor. In con guration shown, animal is stimulated from posterior end by suction electrode, APs propagate toward worm's anterior end, and bipolar electrodes con rm AP stimulation and propagation. (Scale bar 20 mm). B) Measured time trace of AP magnetic eld B(t) from live intact specimen of M. infundibulum for N = 1,650 events. avg B. Nanodiamonds The techniques for the creation of NV centers in diamond are well established also for nanodiamonds (NDs). Nanodiamonds-based sensors exploit colloidal suspensions of single diamond particles of minimum diameter of 4-5 nm, but on average the nanodiamonds typically used in experiment have a size of 50-100 nm. The nanometer size makes nanodiamonds-based sensor of extreme interest for bio-sensing application, as they are potentially usable in vivo experiments. Nonetheless, they have also important drawbacks such as e.g. the increased sensitivity of NV spins to environmental noise. Indeed, while in a bulk diamond the coherence time T is mainly in uenced by the electronic impurities and nuclear spins in the surrounding, for nanodiamonds the coherence time is further reduced due to the surface spin noise. This should be taken into account in the estimation of the sensitivity limit (see Equation (7)). Despite this limitation, nanodiamonds have attracted interest also as a non-toxic alternative to quantum dots for biomedical imaging, as magnetic sensors and nally as drug transporters (thanks to the discovery of the possibility to functionalize the diamonds surface in various ways, exploiting the covalents carbon bonds). The great interest and the exceptional range of applications of NDs is boosting the development of novel fabrication techniques, even if the actual technologies are already able to provide very pure nanodiamonds with controlled surface chemistry at a relative 122,123 low cost . 15 1. Biocompatibility and functionalization studies To understand the perspective in bio-medical application, deep investigation of NDs biocompatibility is required. More speci cally, it is important to understand their di usion in tissues and their acute and long-term biological 124,125 e ect. While bulk diamonds are non-toxic and inert, NDs interaction with cells should be carefully investigated . There is a huge variety of nanodiamond specimens, which di er in dimensions, functionalized surface and potential interaction with the biological sample. 126{128 Various experiments were therefore conducted to evaluate the cell viability, e.g. in HeLa cells (a cell line 43,47 129 deriving from tumoral human cells), in human neurons , in human trachea , in the translucent Caenorhabditis 130 131 elegans worm and intravenous infusions . Brie y, the nanodiamonds of size between 50 and 100 nm have been found to be incorporated by the cells, without producing a signi cant damage. In particular in Guarina et al. an ODMR detection scheme with NV centers in nanodiamonds internalized in hippocampal neurons was performed in suitable conditions (3 mW of excitation power, -20 dBm of continuous- wave MW power). This experiment demonstrated that neuron functionality was not signi cantly a ected by the implementation of the measurement protocol: their spontaneous ring (bursts synchronization) was preserved, as well as the amplitude of spontaneous inhibitory and excitatory events. Even thought some alteration both at the single-cell level and in neuronal networks was observed, this was principally attributed to the e ects of nanoparticles aggregation. The aim of the work was to assess the feasibility of in vitro imaging and targetable drug delivery via nanodiamonds, but the same argument holds for the other sensing applications. Furthermore, if properly functionalized, the NDs can anchor themselves to the surface of the cell sample in the targeted areas. 2. Nanodiamonds applications Once the biocompatibility of nanodiamonds is assessed, it is necessary to understand to which extent the sensing techniques developed for sensor based on NV in bulk diamond can be extended to nanodiamonds based sensors, functionalized and incorporated in the cells of interest. A proof-of-principle demonstration of quantum control techniques to map the intracellular temperature of a neuronal network was performed by Simpson et al. . The NDs were dispersed in cell media in concentration of 6 g/ml, sonicated for few minutes, and then applied to the primary cultures during a routine change of cell media. The 170 nm diameter NDs contained approximately 500 NV centers each. Using ODMR techniques in combination with standard wide- eld microscopy with a eld of view of 80 80 m was possible to observe NV resonance frequency in only 6 seconds. Speci cally, in Ref. the ODMR signal presented two uorescence dips (see section 2.4) because of the strain. This e ect is negligible in bulk diamonds while nanodiamonds crystal lattice su ers strong deformation inducing line splitting. In that paper the two dips, spaced by few MHz, were modeled as a single one with higher spectral broadening. By interpolating the ODMR graph with a Lorentzian function, it was estimated the mean crystal eld splitting D = (2868.590.17) MHz. gs To demonstrate the NV thermo-sensor performance in biological measurement, the temperature of the neuronal solution was reduced by 1:9 C. Repeating the ODMR analysis for a total acquisition time of 12 s, a resonance frequency shift was observed. The respective temperature variation was estimated using the temperature coupling coecient dD=dT ' -74 kHz K (see section 2.5). The distribution reported a mean temperature change of (1:36 0:08) C, consistent with the reduction in environmental temperature. We underline that NDs allowing to create spatial maps of the temperature inside the cells will bring new insight on the understanding of cell activity. There are many biological processes whose knowledge would be enriched by nanoscale thermometry, such as temperature increases following the opening of ion channels , or the correlation of temperature changes and the occurrence of neurological disorders and pathological conditions. Another biological application was reported by Ermakova et al. , using nanodiamonds with NV centers as 132,133 thermo-sensors, exploits optically-induced thermal gradients for thermogenetic neural modulation . This thermal gradient is generated at the transient receptor potential channels (TRP channels): a group of ion channels that are commonly present on the plasma membrane of numerous types of animal cells . A particular specialized form of these ion channels appears to be highly sensitive to temperature changes . Some species of snakes can use TRP channels to detect the thermal build-up caused by infrared IR radiation emitted by nearby prey, allowing them to estimate the direction and distance of the (IR) source . To experimentally recreate this local temperature change and therefore study the TRPs response, Ermakova et al. used IR short pulsed laser. This method, with respect to conventional techniques as environmental heating or TRPs chemical agonists , allows cellular spatial resolution and ultrahigh temporal resolution. The precise temperature control was performed by varying the laser intensity, whose actual thermal impact was monitored by 16 the nitrogen-vacancy complex. This quantum probe (whose dimension was about 300 nm) was integrated on the tip of an optical ber, together with a microwave antenna. The optical ber was positioned near the cell irradiated by the IR laser, allowing a measurement of its temperature by the ODMR technique. In this experiment it was initially evaluated the thermal stimulation via IR laser of two TRP channels of the snake. The TRP channels considered were the Crotalus atrox TRPA1 (caTRPA1) and the Elaphe obsoleta lindheimeri TRPA1 (eolTRPA1). Fluorescent proteins (caTRPA1-IRES-EGFP) had been added to the channels, allowing to monitor the opening and closing of the calcium channels. Thanks to the NV sensors and by slowly changing the cell temperature with properly tuning pulsed laser intensity, it was possible to obtain the threshold temperature, inducing opening of the calcium channels. The threshold temperatures were found to be T =(27.8  0.6) C for caTRPA1 (see Figure 11b) and T =(38.5  0.7) C for eolTRPA1 (see Figure 11c). Once estimated the threshold temperature T , Ermakova et al. proved the technique on other biological samples: mouse neurons and zebra sh larvae, whose thermogenetic activation is induced by TRPA1 channels causing responses. In the case of caTRPA1 channels, the cultured neurons were maintained at a temperature of 27 C lower than the threshold temperature obtained before for this channel; in the case of eolTRPA1-expressing, neurons were kept at basal temperatures of 35.5 C. As expected, they found that the thermal increase induced by the IR laser activates the TRPs channels triggering the generation of the neuronal AP, measured through conventional electrophysiological techniques. When a measurements on live samples is considered, the sample can no longer be kept at the desired temperature, therefore it is necessary to choose the TRP channel suitable for body temperature of the animal species analyzed. As for the zebra sh neurons, whose body temperature is found to be 26 C, the eolTRPA1 channels may be suitable. As for the mammalian brain, the perfect TRP candidate has still to be found. For example, the mouse body temperature is too close to the threshold temperature of eolTRPA1 and it may be desensitized. The results of the application of this technique in living zebra sh showed that it is possible to thermogenetically activate neurons using the IR laser. In particular, the technique demonstrated a spatial resolution of 60 m ( ber size in which the IR laser was focused on the sample), allowing one or few neurons to be stimulated. As for the IR laser intensity, Emarkova et al. observed that 30 mW laser power induced the escape behavior exhibition of 93% of the larvae. NV-based temperature sensors allowed careful monitoring of the temperature reached by the cells with high spatial resolution and temperature sensitivity up to 0.1 C. To preserve cellular integrity and to avoid cell ablation is essential to heat-up the tissues by a few degrees only and for a time interval not exceeding a few minutes. FIG. 11: Activation of snake TRPA1 in cells expressing TRPA1-IRES-EGFP using femtosecond IR laser pulses, taken from the reference . 2+ a) R-GECO1.1 uorescence (black line) re ects Ca dynamics in the cytoplasm with the 20 mW laser beam turned on at t=30 s and o at t=60 s. b,c) With the temperature of HEK293 cells expressing snake TRPA1 increased in a stepwise fashion using properly adjusted IR laser radiation, the activation thresholds of caTRPA1 b) and eolTRPA1 c) were determined. 2+ d) A similar heating of control cells does not induce Ca elevation. The black line is the uorescence response. The red line is the temperature in the medium. Finally, another demonstration of the e ectiveness of NV nanosensors for thermometry comes from Fujiwara et al.'s experiment . A rst test allowed to measure the temperature dynamics inside live C. elegans adults worms 1=2 during environmental temperature changes. The sensitivity reached was 1.4 C Hz . Having obtained this result, Fujiwara et al. successfully determined the temperature increase caused by the worm's thermogenesis under the treatment of mitochondrial uncoupler stimuli. 17 V. TECHNIQUES FOR IMPROVING ODMR SENSITIVITY In this section we discuss some technological solutions to improve the sensitivity of NV-based sensors as well as the precautions to be taken when they are used as a bio-sensor. Equation (7) provides the ultimate sensitivity limit reachable highlighting that the number n of NV centers and their decoherence time T play a key role. To increase n, while maintaining the same spatial resolution, it is necessary to have diamonds with an increased NV centers density. This can be achieved by enhancing the number of nitrogen implanted in the diamond and improving the N-to-NV conversion eciency, minimizing the concentration of residual paramagnetic substitutional nitrogen . In parallel, to increase T , it is also recommended the production of ultra-pure diamonds, with reduced unwanted electronic impurities (e.g. the P1 centers) and nuclear spins impurities 13 139{141 (e.g. the paramagnetic C isotopes, whose natural abundance is about 1.1%) . It is important to note that the NV density increase will necessary worsen the decoherence time of the NVs themselves, because of their mutual interaction. Consequently, an optimal trade-o s between these parameters must be sought. In addition to the NV-density and diamond sample engineering, the sensitivity can be improved by imple- menting speci c experimental techniques, that are based on laser and microwave pulses of particular duration, 28,46,142 synchronized appropriately . For example, if an unknown electromagnetic eld, responsible for the ODMR resonance frequency shift, is constant or 143 38 slowly varying, it is possible to adopt the experimental pulsed ODMR protocols or the Ramsey method instead of the continuous wave (CW ) ODMR . The CW ODMR is the simplest and most widely employed magnetometry method with NV-based sensors, wherein the microwave driving and the optical polarization and readout (laser pumping) occur simultaneously. Although this technique is easy to be implemented, the relative ODMR spectrum dips are a ected by the broadening induced by the continuous exposure of the laser beam and microwave eld on the sample. With pulsed ODMR techniques this broadening e ect is substantially suppressed, allowing to obtain a narrower ODMR spectrum dips and therefore to improve the measurement sensitivity. This protocol uses temporally separated optical laser initializations,  microwave control pulses, and laser readout pulses. The  pulses, whose name derives from the representation of the process on the Bloch sphere, is an oscillating microwave eld that brings the electronic state from the state jm = 0 > to jm = 1 >. Ramsey ODMR spectroscopy, on the other s s hand, consists on in the application of two =2 pulses, separated by a time  . Also the =2 pulse is an oscillating microwave eld that brings the electronic state from the state jm = 0 > to a balanced superposition of jm = +1 > s s and jm = 1 >. By varying the time  , the so-called "Ramsey fringes" are obtained, from which it is possible to extrapolate an estimation of the magnetic elds amplitude. Also this technique allows sensitivity improvement with respect to the CW : the decoupling of the MW from the laser power allows increasing the MW power improving the contrast, without degrading the FWHM. In the case of time-varying electromagnetic elds, there are other even more complex microwave pulse sequences, capable of decoupling the measurement from surrounding spin environment . In this way the decoherence time of the NV centers increases and consequently it becomes possible to interrogate the quantum system for longer times, improving the measurement statistic and therefore the sensitivity. One of these experimental protocols is the 127,146 Hahn Echo sequence , which refocuses the dephasing NVs spin, applying an additional  pulse in the middle of Ramsey sequence. The characteristic time of the spin coherence decay, measured with this protocol, is called T and it is typically one or two orders of magnitude longer than T . Even more complex dynamic decoupling sequences, 147{150 which apply multiple refocusing  pulses further improving the decoherence time T have been devised . Among 147{150 147,151 these, the most famous are the Carr-Purcell-Meiboom-Gill (CPMG) and the XY8 sequences , which di er in the rotation axes (around which the spin rotates): the rst method applies the pulses along the same axis, while the second chooses a di erent one for each  pulses. It is useful to underline that, although these techniques allow to extended the coherence time of the NV centers, they cannot go beyond the spin-lattice relaxation time T , that for an NV spin ensemble in bulk diamond is about 3 ms . The Figure 12 brie y summarizes the above mentioned pulse sequences. It is useful to underline that the sensitivity formula in Equation (7) describes an idealized measurement with a perfect readout mechanism. On the contrary, typically the readout mechanism adds noise in the measurement, that can be described introducing, in the previous equation, the spin-readout delity factor F : 1 1 = (9) nT Keeping the usual optical-readout, but improving the photon collection is expected to increase F (see ref. for di erent methods to improve photon collection). Ancilla-assisted repetitive readout, which is based on mapping the NV spin state to the nuclear spin state, also improves F . Finally, in a more far perspective, quantum methods of 154{157 noise reduction can be applied . 18 CW PULSED RAMSEY HANH ECHO DYNAMICAL DECOULPLING LASER MW ( ) π π/2 π/2 π/2 π π/2 π/2 π π/2 READOUT FIG. 12: Scheme of timing and duration of laser pulses, microwave pulses and reading sequences associated with the most common measurement protocols of the external elds for the NV complex. When the ultimate goal is bio-sensing, some constraints rise limiting the implementation of the above de- scribed pulse sequences. One constraint is the frequency bandwidth. In fact, the dynamic decoupling techniques mentioned above are capable of measuring time-varying external eld only if this time variation is of the order of the time interval separating the  pulses. Furthermore, in order to control the system quantum state, the time between these pulses cannot exceed the coherence of the NV center. Consequently, the frequency of the signal to be measured must be of the order of the coherence time of the NV centers. In the biological case, the electromagnetic elds pulse lasts about 1 ms. This value is very far from T , marking a boundary for the use of these techniques in biological applications. Another constraint is associated to the optical laser power. The higher the laser power the better the sensitivity in measurements with ensembles, since it increases the percentage of excited centers and consequently the uorescence signal. However precautions must be taken to avoid cells and proteins damaging. An ecient solution can be to direct the laser beam towards the diamond sample at an angle allowing total re ection (Brewster angle). In this way only the uorescence emitted by the NV centers travels through the cells, placed on the other diamond 44,77 surface . In this case, however, precise control over sensing volume would be lost, deteriorating spatial resolution. In a standard con guration, where the laser impinges perpendicularly on the sample, it is necessary to limit the optical power reaching the cells to few mW. In this regard, Figure 13 shows a sensitivity curve versus the laser optical power, obtained by adopting the technique described in Moreva et al. . As anticipated in the introduction, FIG. 13: Temperature sensitivity versus the laser excitation power at 532 nm. The inset shows the inverse of the thermal sensitivity versus the excitation laser power. 19 the application of a transverse bias magnetic eld B ' 3 mT , allows to improve the sensitivity of a the NV center bias based thermo-sensor with respect to other standard techniques in CW regime. In Ref. , the temperature sensitivity 1=2 3 reached is  ' 4:8 mK=Hz in a sensing volume of 1m , obtained at a power level (80 mW) that can present biocompatibility problems. However, the sensitivity obtained is even beyond the one required to monitor biological mechanisms, usually requiring sensitivities of the order of 1 C. Figure 13 shows that it is possible to perform the temperature measurement with a lower laser power, nding an ideal compromise between the temperature sensitivity and laser intensity impinging on the cell sample. Indeed, with a power of a few mW it is already possible to discriminate biological processes with a sensitivity of the tenth of a degree. VI. CONCLUSION Sensors based on the NV centers in arti cial diamonds are one of the emerging quantum technologies of huge potential interest in biological applications, thanks to both their practicality and their technical performances. In fact, the capability to initialize and read out optically the spin state at room temperature, makes the use of this quantum sensor convenient and powerful even for biological applications. Furthermore, the levels of sensitivity and spatial resolution achieved are extremely high, which in principle allows potential application towards the detection of very weak electromagnetic elds as the one generated by mammalian, and potentially human, cells. Even if an eventual use of NV sensors for the detection of biological electric elds is more problematic due to its weak coupling constant, regarding the magnetic eld sensing and especially temperature measurements astonishing results have already been achieved. Indeed, the thermal gradients generated by biological phenomena can be reliably observed thanks to the actual sensitivity of the NV-based sensors. This is of the utmost importance because localized intracellular temperature gradients may a ect neuronal functionality (including vesicular dynamics and neurotransmitter release) 158,159 or may provide indirect measurement of mitochondrial activity . Regarding the detection of bio-magnetic elds, the NV-based sensors have already shown good results with peculiar biological cells, presenting either an intrinsic magnetic eld (magnetitotactic bacteria) or a generated magnetic eld in axon of squids or long worms, much larger than the one generated in the human ones. The improvement of these devices suggests the possibility of exploiting NV-based sensors also for the detection of weaker but more fascinating biological magnetic elds. In particular, an estimate of the cardiac magnetic eld that is generated on the heart surface was here reported. This value is in the range of present measurement capability exploiting the NV center properties, exploiting optimized diamond sample engineering and the adoption of pulsed measurement protocols in order to improve the diamond coherence time. Furthermore, we have analyzed the magnetic eld associated to human neuron activity. The weakness of these elds requires further improvements of the measurement technique in the case of the single AP, while measurement of clustered channels is likely a reasonable target for the actual technology. However, the considerable interest in the neuronal eld detection as diagnostic and therapeutic tools for neurodegenerative diseases and aging e ects, together with the recent years progress of these techniques (partially covered by this review), is expected to boost the technological developments and eventually the market success of quantum assisted biosensing based on NVs. Acknowledgements This work has received funding from the European Union's PATHOS EU H2020 FET-OPEN grant no. 828946 and Horizon 2020, from the EMPIR Participating States in the context of the projects EMPIR-17FUN06 "SIQUST" and from the project Piemonte Quantum Enabling Technologies (PiQuET) funded by the Piemonte Region. Richard T Mathias, Ira S Cohen, and Carlos Oliva. Limitations of the whole cell patch clamp technique in the control of intracellular concentrations. Biophysical journal, 58(3):759, 1990. Michael Fejtl, Alfred Stett, Wilfried Nisch, Karl-Heinz Boven, and Andreas M oller. On micro-electrode array revival: its development, sophistication of recording, and stimulation. In Advances in network electrophysiology, pages 24{37. Springer, Micha E Spira and Aviad Hai. Multi-electrode array technologies for neuroscience and cardiology. Nature nanotechnology, 8(2):83, 2013. M Vengalattore, JM Higbie, SR Leslie, J Guzman, LE Sadler, and DM Stamper-Kurn. High-resolution magnetometry with a spinor bose-einstein condensate. Physical review letters, 98(20):200801, 2007. MI Faley, U Poppe, K Urban, DN Paulson, and RL Fagaly. A new generation of the hts multilayer dc-squid magnetometers and gradiometers. In Journal of Physics: Conference Series, volume 43, page 1199. IOP Publishing, 2006. F Baudenbacher, LE Fong, JR Holzer, and M Radparvar. Monolithic low-transition-temperature superconducting magne- tometers for high resolution imaging magnetic elds of room temperature samples. Applied Physics Letters, 82(20):3487{ 3489, 2003. 20 Svenja Knappe, Tilmann H Sander, Olaf Kosch, Frank Wiekhorst, John Kitching, and Lutz Trahms. Cross-validation of microfabricated atomic magnetometers with superconducting quantum interference devices for biomagnetic applications. Applied Physics Letters, 97(13):133703, 2010. Anri Nakajima. Application of single-electron transistor to biomolecule and ion sensors. Applied Sciences, 6(4):94, 2016. Laura Martinez Maestro, Emma Martin Rodriguez, Francisco Sanz Rodriguez, MC Iglesias-de la Cruz, Angeles Juarranz, Ra k Naccache, Fiorenzo Vetrone, Daniel Jaque, John A Capobianco, and Jose Garcia Sole. Cdse quantum dots for two-photon uorescence thermal imaging. Nano letters, 10(12):5109{5115, 2010. CDS Brites, A Mill an, and LD Carlos. Lanthanides in luminescent thermometry. In Handbook on the Physics and Chemistry of Rare Earths, volume 49, pages 339{427. Elsevier, 2016. Michael Hines. Neuron|a program for simulation of nerve equations. In Neural systems: Analysis and modeling, pages 127{136. Springer, 1993. Michael L Hines and Nicholas T Carnevale. The neuron simulation environment. Neural computation, 9(6):1179{1209, F Santamaria and JM Bower. Hodgkin-huxley models. In Encyclopedia of Neuroscience, pages 1173{1180. Elsevier Ltd, Blaise Aguera  y Arcas, Adrienne L Fairhall, and William Bialek. Computation in a single neuron: Hodgkin and huxley revisited. Neural Computation, 15(8):1715{1749, 2003. Single-Channel Recording. edited by b. sakmann and n. neher, 1995. Stephen G Waxman. Demyelinating diseases|new pathological insights, new therapeutic targets, 1998. Eugene A Kiyatkin. Brain temperature: from physiology and pharmacology to neuropathology. In Handbook of clinical neurology, volume 157, pages 483{504. Elsevier, 2018. PC Maurer, M Kubo, MD Lukin, HJ Noh, G Kucsko, NY Yao, PK Lo, and H Park. Nanometer scale quantum thermometry in a living cell. 2013. Tingting Bai and Ning Gu. Micro/nanoscale thermometry for cellular thermal sensing. Small, 12(34):4590{4610, 2016. Yuzhou Wu, Fedor Jelezko, Martin B Plenio, and Tanja Weil. Diamond quantum devices in biology. Angewandte Chemie International Edition, 55(23):6586{6598, 2016. Alexios Beveratos, Rosa Brouri, Thierry Gacoin, Jean-Philippe Poizat, and Philippe Grangier. Nonclassical radiation from diamond nanocrystals. Physical Review A, 64(6):061802, 2001. Alexios Beveratos, Sergei Kuhn,  Rosa Brouri, Thierry Gacoin, J-P Poizat, and Philippe Grangier. Room temperature stable single-photon source. The European Physical Journal D-Atomic, Molecular, Optical and Plasma Physics, 18(2):191{196, Alexander M Zaitsev. Optical properties of diamond: a data handbook. Springer Science & Business Media, 2013. Ettore Bernardi, Ekaterina Moreva, Paolo Traina, Giulia Petrini, Sviatoslav Ditalia Tchernij, Jacopo Forneris, Zelijko Pastuovic, Ivo Pietro Degiovanni, Paolo Olivero, and M Genovese. Biocompatible technique for nanoscale magnetic eld sensing with nitrogen-vacancy centers. arXiv preprint arXiv:2005.13230, 2020. Shu-Jung Yu, Ming-Wei Kang, Huan-Cheng Chang, Kuan-Ming Chen, and Yueh-Chung Yu. Bright uorescent nanodia- monds: no photobleaching and low cytotoxicity. Journal of the American Chemical Society, 127(50):17604{17605, 2005. Frank C Barone, Cezary Marcinkiewicz, Jie Li, Yi Feng, Mark Sternberg, Peter I Lelkes, David Rosenbaum-Halevi, Jonathan A Gerstenhaber, and Giora Z Feuerstein. Long-term biocompatibility of uorescent diamonds-(nv)-z~ 800 nm in rats: survival, morbidity, histopathology, particle distribution and excretion studies (part iv). International journal of nanomedicine, 14:1163, 2019. Christian G Specht, Oliver A Williams, Richard B Jackman, and Ralf Schoepfer. Ordered growth of neurons on diamond. Biomaterials, 25(18):4073{4078, 2004. Romana Schirhagl, Kevin Chang, Michael Loretz, and Christian L Degen. Nitrogen-vacancy centers in diamond: nanoscale sensors for physics and biology. Annual review of physical chemistry, 65:83{105, 2014. Georg Kucsko, Peter C Maurer, Norman Ying Yao, MICHAEL Kubo, Hyun Jong Noh, Po Kam Lo, Hongkun Park, and Mikhail D Lukin. Nanometre-scale thermometry in a living cell. Nature, 500(7460):54{58, 2013. Jan Jeske, Jared H Cole, and Andrew D Greentree. Laser threshold magnetometry. New Journal of Physics, 18(1):013015, Marcus W Doherty, Neil B Manson, Paul Delaney, Fedor Jelezko, J org Wrachtrup, and Lloyd CL Hollenberg. The nitrogen- vacancy colour centre in diamond. Physics Reports, 528(1):1{45, 2013. Thomas Wolf, Philipp Neumann, and J org Wrachtrup. Diamond magnetometer, December 26 2017. US Patent 9,851,418. A Gruber, A Dr abenstedt, C Tietz, L Fleury, J Wrachtrup, and C Von Borczyskowski. Scanning confocal optical microscopy and magnetic resonance on single defect centers. Science, 276(5321):2012{2014, 1997. MW Doherty, F Dolde, H Fedder, Fedor Jelezko, J Wrachtrup, NB Manson, and LCL Hollenberg. Theory of the ground-state spin of the nv- center in diamond. Physical Review B, 85(20):205203, 2012. Gopalakrishnan Balasubramanian, IY Chan, Roman Kolesov, Mohannad Al-Hmoud, Julia Tisler, Chang Shin, Changdong Kim, Aleksander Wojcik, Philip R Hemmer, Anke Krueger, et al. Nanoscale imaging magnetometry with diamond spins under ambient conditions. Nature, 455(7213):648{651, 2008. Linh My Pham, David Le Sage, Paul L Stanwix, Tsun Kwan Yeung, D Glenn, Alexei Trifonov, Paola Cappellaro, Philip R Hemmer, Mikhail D Lukin, Hongkun Park, et al. Magnetic eld imaging with nitrogen-vacancy ensembles. New Journal of Physics, 13(4):045021, 2011. Jennifer M Schloss, John F Barry, Matthew J Turner, and Ronald L Walsworth. Simultaneous broadband vector magne- tometry using solid-state spins. Physical Review Applied, 10(3):034044, 2018. 21 JM Taylor, P Cappellaro, L Childress, L Jiang, D Budker, PR Hemmer, A Yacoby, R Walsworth, and MD Lukin. High- sensitivity diamond magnetometer with nanoscale resolution. Nature Physics, 4(10):810{816, 2008. CL Degen. Scanning magnetic eld microscope with a diamond single-spin sensor. Applied Physics Letters, 92(24):243111, Jeronimo R Maze, Paul L Stanwix, James S Hodges, Seungpyo Hong, Jacob M Taylor, Paola Cappellaro, Liang Jiang, MV Gurudev Dutt, Emre Togan, AS Zibrov, et al. Nanoscale magnetic sensing with an individual electronic spin in diamond. Nature, 455(7213):644{647, 2008. Florian Dolde, Helmut Fedder, Marcus W Doherty, Tobias N obauer, Florian Rempp, Gopalakrishnan Balasubramanian, Thomas Wolf, Friedemann Reinhard, Lloyd CL Hollenberg, Fedor Jelezko, et al. Electric- eld sensing using single diamond spins. Nature Physics, 7(6):459{463, 2011. Victor M Acosta, Erik Bauch, Micah P Ledbetter, Amir Waxman, L-S Bouchard, and Dmitry Budker. Temperature dependence of the nitrogen-vacancy magnetic resonance in diamond. Physical review letters, 104(7):070801, 2010. Laura Guarina, C Calorio, D Gavello, E Moreva, P Traina, A Battiato, S Ditalia Tchernij, J Forneris, M Gai, F Picollo, et al. Nanodiamonds-induced e ects on neuronal ring of mouse hippocampal microcircuits. Scienti c Reports, 8(1):1{14, David Le Sage, Koji Arai, David R Glenn, Stephen J DeVience, Linh M Pham, Lilah Rahn-Lee, Mikhail D Lukin, Amir Yacoby, Arash Komeili, and Ronald L Walsworth. Optical magnetic imaging of living cells. Nature, 496(7446):486{489, Hunter C Davis, Pradeep Ramesh, Aadyot Bhatnagar, Audrey Lee-Gosselin, John F Barry, David R Glenn, Ronald L Walsworth, and Mikhail G Shapiro. Mapping the microscale origins of magnetic resonance image contrast with subcellular diamond magnetometry. Nature communications, 9(1):1{9, 2018. John F Barry, Jennifer M Schloss, Erik Bauch, Matthew J Turner, Connor A Hart, Linh M Pham, and Ronald L Walsworth. Sensitivity optimization for nv-diamond magnetometry. Reviews of Modern Physics, 92(1):015004, 2020. David A Simpson, Emma Morrisroe, Julia M McCoey, Alain H Lombard, Dulini C Mendis, Francois Treussart, Liam T Hall, Steven Petrou, and Lloyd CL Hollenberg. Non-neurotoxic nanodiamond probes for intraneuronal temperature mapping. ACS nano, 11(12):12077{12086, 2017. Yulia G Ermakova, Aleksandr A Lanin, Ilya V Fedotov, Matvey Roshchin, Ilya V Kelmanson, Dmitry Kulik, Yulia A Bogdanova, Arina G Shokhina, Dmitry S Bilan, Dmitry B Staroverov, et al. Thermogenetic neurostimulation with single- cell resolution. Nature communications, 8(1):1{15, 2017. Masazumi Fujiwara, Simo Sun, Alexander Dohms, Yushi Nishimura, Ken Suto, Yuka Takezawa, Keisuke Oshimi, Li Zhao, Nikola Sadzak, Yumi Umehara, Yoshio Teki, Naoki Komatsu, Oliver Benson, Yutaka Shikano, and Eriko Kage-Nakadai. Real-time nanodiamond thermometry probing in-vivo thermogenic responses, 2020. JoHoNo Loubser and JoA van Wyk. Electron spin resonance in the study of diamond. Reports on Progress in Physics, 41(8):1201, 1978. Erik Bauch, Connor A Hart, Jennifer M Schloss, Matthew J Turner, John F Barry, Pauli Kehayias, Swati Singh, and Ronald L Walsworth. Ultralong dephasing times in solid-state spin ensembles via quantum control. Physical Review X, 8(3):031025, 2018. Chang S Shin, Mark C Butler, Hai-Jing Wang, Claudia E Avalos, Scott J Seltzer, Ren-Bao Liu, Alexander Pines, and Vikram S Bajaj. Optically detected nuclear quadrupolar interaction of n 14 in nitrogen-vacancy centers in diamond. Physical Review B, 89(20):205202, 2014. S Felton, AM Edmonds, ME Newton, PM Martineau, D Fisher, DJ Twitchen, and JM Baker. Hyper ne interaction in the ground state of the negatively charged nitrogen vacancy center in diamond. Physical Review B, 79(7):075203, 2009. Xing-Fei He, Neil B Manson, and Peter TH Fisk. Paramagnetic resonance of photoexcited n-v defects in diamond. ii. hyper ne interaction with the n 14 nucleus. Physical Review B, 47(14):8816, 1993. Victor M Acosta, Erik Bauch, Micah P Ledbetter, Charles Santori, K-MC Fu, Paul E Barclay, Raymond G Beausoleil, H elo se Linget, Jean Francois Roch, Francois Treussart, et al. Diamonds with a high density of nitrogen-vacancy centers for magnetometry applications. Physical Review B, 80(11):115202, 2009. Ekaterina Moreva. The biosensing with nv centers in diamond: Related challenges. International Journal of Quantum Information, 18(01):1941023, 2020. Varun KA Sreenivasan, Andrei V Zvyagin, and Ewa M Goldys. Luminescent nanoparticles and their applications in the life sciences. Journal of Physics: Condensed Matter, 25(19):194101, 2013. Eric Van Oort and Max Glasbeek. Electric- eld-induced modulation of spin echoes of nv centers in diamond. Chemical Physics Letters, 168(6):529{532, 1990. X-D Chen, C-H Dong, F-W Sun, C-L Zou, J-M Cui, Z-F Han, and G-C Guo. Temperature dependent energy level shifts of nitrogen-vacancy centers in diamond. Applied Physics Letters, 99(16):161903, 2011. Yan-Kai Tzeng, Pei-Chang Tsai, Hsiou-Yuan Liu, Oliver Y Chen, Hsiang Hsu, Fu-Goul Yee, Ming-Shien Chang, and Huan-Cheng Chang. Time-resolved luminescence nanothermometry with nitrogen-vacancy centers in nanodiamonds. Nano letters, 15(6):3945{3952, 2015. E. Moreva, E. Bernardi, P. Traina, A. Sosso, S. Ditalia Tchernij, J. Forneris, F. Picollo, G. Brida, Z. Pastuovi c, I. P. Degiovanni, P. Olivero, and M. Genovese. Practical applications of quantum sensing: A simple method to enhance the sensitivity of nitrogen-vacancy-based temperature sensors. Phys. Rev. Applied, 13:054057, May 2020. E Moreva, E Bernardi, P Traina, A Sosso, S Ditalia Tchernij, J Forneris, F Picollo, G Brida, Z Pastuovic, IP Degio- vanni, et al. Practical applications of quantum sensing: a simple method to enhance sensitivity of nitrogen-vacancy-based temperature sensors. arXiv preprint arXiv:1912.10887, 2019. 22 Erwin Neher and Bert Sakmann. The patch clamp technique. Scienti c American, 266(3):44{51, 1992. Alfred Stett, Ulrich Egert, Elke Guenther, Frank Hofmann, Thomas Meyer, Wilfried Nisch, and Hugo Haemmerle. Biological application of microelectrode arrays in drug discovery and basic research. Analytical and bioanalytical chemistry, 377(3):486{ 495, 2003. Jubayer Jalil, Yong Zhu, Chandima Ekanayake, and Yong Ruan. Sensing of single electrons using micro and nano tech- nologies: a review. Nanotechnology, 28(14):142002, 2017. Yoshiaki Takei, Satoshi Arai, Atsushi Murata, Masao Takabayashi, Kotaro Oyama, Shin'ichi Ishiwata, Shinji Takeoka, and Madoka Suzuki. A nanoparticle-based ratiometric and self-calibrated uorescent thermometer for single living cells. ACS nano, 8(1):198{206, 2014. Duan P Chen, Robert S Eisenberg, Joseph W Jerome, and Chi-Wang Shu. Hydrodynamic model of temperature change in open ionic channels. Biophysical Journal, 69(6):2304, 1995. Ahmed El Hady and Benjamin B Machta. Mechanical surface waves accompany action potential propagation. Nature communications, 6:6697, 2015. Ezia Guatteo, Kenny KH Chung, Tharushini K Bowala, Giorgio Bernardi, Nicola B Mercuri, and Janusz Lipski. Tempera- ture sensitivity of dopaminergic neurons of the substantia nigra pars compacta: involvement of transient receptor potential channels. Journal of neurophysiology, 94(5):3069{3080, 2005. Jon S Donner, Sebastian A Thompson, Mark P Kreuzer, Guillaume Ba ou, and Romain Quidant. Mapping intracellular temperature using green uorescent protein. Nano letters, 12(4):2107{2111, 2012. Jui-Ming Yang, Haw Yang, and Liwei Lin. Quantum dot nano thermometers reveal heterogeneous local thermogenesis in living cells. ACS nano, 5(6):5067{5071, 2011. Soo Ho Kim, Jermim Noh, Min Ku Jeon, Ki Woong Kim, Luke P Lee, and Seong Ihl Woo. Micro-raman thermome- try for measuring the temperature distribution inside the microchannel of a polymerase chain reaction chip. Journal of Micromechanics and Microengineering, 16(3):526, 2006. Fiorenzo Vetrone, Ra k Naccache, Alicia Zamarron, Angeles Juarranz de la Fuente, Francisco Sanz-Rodriguez, Laura Martinez Maestro, Emma Martin Rodriguez, Daniel Jaque, Jose Garciia Sole, and John A Capobianco. Temperature sensing using uorescent nanothermometers. ACS nano, 4(6):3254{3258, 2010. Laura Martinez Maestro, Qiming Zhang, Xiangping Li, Daniel Jaque, and Min Gu. Quantum-dot based nanothermometry in optical plasmonic recording media. Applied Physics Letters, 105(18):181110, 2014. Ellen C Jensen. Use of uorescent probes: their e ect on cell biology and limitations. The Anatomical Record: Advances in Integrative Anatomy and Evolutionary Biology, 295(12):2031{2036, 2012. Tytus Bernas, Bartlomiej P Rajwa, Elikplimki K Asem, and Joseph Paul Robinson. Loss of image quality in photobleaching during microscopic imaging of uorescent probes bound to chromatin. Journal of biomedical optics, 10(6):064015, 2005. John F Barry, Matthew J Turner, Jennifer M Schloss, David R Glenn, Yuyu Song, Mikhail D Lukin, Hongkun Park, and Ronald L Walsworth. Optical magnetic detection of single-neuron action potentials using quantum defects in diamond. Proceedings of the National Academy of Sciences, 113(49):14133{14138, 2016. DM Toyli, DJ Christle, A Alkauskas, BB Buckley, CG Van de Walle, and DD Awschalom. Measurement and control of single nitrogen-vacancy center spins above 600 k. Physical Review X, 2(3):031001, 2012. Taras Plakhotnik, Marcus W Doherty, Jared H Cole, Robert Chapman, and Neil B Manson. All-optical thermometry and thermal properties of the optically detected spin resonances of the nv{center in nanodiamond. Nano letters, 14(9):4989{4996, Jan Muller,  Marco Ballini, Paolo Livi, Yihui Chen, Milos Radivojevic, Amir Shadmani, Vijay Viswam, Ian L Jones, Michele Fiscella, Roland Diggelmann, et al. High-resolution cmos mea platform to study neurons at subcellular, cellular, and network levels. Lab on a Chip, 15(13):2767{2780, 2015. Gabriel Bertotti, Dmytro Velychko, Norman Dodel, Stefan Keil, Dirk Wolansky, Bernd Tillak, Matthias Schreiter, Andreas Grall, Peter Jesinger, Sebastian R ohler, et al. A cmos-based sensor array for in-vitro neural tissue interfacing with 4225 recording sites and 1024 stimulation sites. In 2014 IEEE Biomedical Circuits and Systems Conference (BioCAS) Proceedings, pages 304{307. IEEE, 2014. Douglas J Bakkum, Urs Frey, Milos Radivojevic, Thomas L Russell, Jan Muller,  Michele Fiscella, Hirokazu Takahashi, and Andreas Hierlemann. Tracking axonal action potential propagation on a high-density microelectrode array across hundreds of sites. Nature communications, 4(1):1{12, 2013. Gy orgy Buzs aki, Eran Stark, Antal Ber enyi, Dion Khodagholy, Daryl R Kipke, Euisik Yoon, and Kensall D Wise. Tools for probing local circuits: high-density silicon probes combined with optogenetics. Neuron, 86(1):92{105, 2015. Milos Radivojevic, David J ackel, Michael Altermatt, Jan Muller,  Vijay Viswam, Andreas Hierlemann, and Douglas J Bakkum. Electrical identi cation and selective microstimulation of neuronal compartments based on features of extracellular action potentials. Scienti c reports, 6:31332, 2016. Giulia Tomagra, Pietro Apr a, Al o Battiato, Cecilia Colla Ruvolo, Alberto Pasquarelli, Andrea Marcantoni, Emilio Car- bone, Valentina Carabelli, Paolo Olivero, and Federico Picollo. Micro graphite-patterned diamond sensors: Towards the simultaneous in vitro detection of molecular release and action potentials generation from excitable cells. Carbon, 152:424{ 433, 2019. Matthew N Rasband and James S Trimmer. Developmental clustering of ion channels at and near the node of ranvier. Developmental biology, 236(1):5{16, 2001. Daisuke Sato, Gonzalo Hern andez-Hern andez, Collin Matsumoto, Sendoa Tajada, Claudia M Moreno, Rose E Dixon, Samantha O'Dwyer, Manuel F Navedo, James S Trimmer, Colleen E Clancy, et al. A stochastic model of ion channel cluster formation in the plasma membrane. Journal of General Physiology, 151(9):1116{1134, 2019. 23 Bokai Zhang, Xi Feng, Hang Yin, Zhenpeng Ge, Yanhuan Wang, Zhiqin Chu, Helena Raabova, Jan Vavra, Petr Cigler, Renbao Liu, et al. Anchored but not internalized: shape dependent endocytosis of nanodiamond. Scienti c reports, 7:46462, Torsten Rendler, Jitka Neburkova, Ondrej Zemek, Jan Kotek, Andrea Zappe, Zhiqin Chu, Petr Cigler, and J org Wrachtrup. Optical imaging of localized chemical events using programmable diamond quantum nanosensors. Nature communications, 8(1):1{9, 2017. Jasmina Isakovic, Ian Dobbs-Dixon, Dipesh Chaudhury, and Dinko Mitrecic. Modeling of inhomogeneous electromagnetic elds in the nervous system: a novel paradigm in understanding cell interactions, disease etiology and therapy. Scienti c reports, 8(1):1{20, 2018. Mursel  Karadas, Adam M Wojciechowski, Alexander Huck, Nils Ole Dalby, Ulrik Lund Andersen, and Axel Thielscher. Feasibility and resolution limits of opto-magnetic imaging of neural network activity in brain slices using color centers in diamond. Scienti c reports, 8(1):1{14, 2018. LT Hall, GCG Beart, EA Thomas, DA Simpson, LP McGuinness, JH Cole, JH Manton, RE Scholten, Fedor Jelezko, J org Wrachtrup, et al. High spatial and temporal resolution wide- eld imaging of neuron activity using quantum nv-diamond. Scienti c reports, 2:401, 2012. Matti H am al ainen, Riitta Hari, Risto J Ilmoniemi, Jukka Knuutila, and Olli V Lounasmaa. Magnetoencephalogra- phy|theory, instrumentation, and applications to noninvasive studies of the working human brain. Reviews of modern Physics, 65(2):413, 1993. Rolf Simon Schoenfeld and Wolfgang Harneit. Real time magnetic eld sensing and imaging using a single spin in diamond. Physical review letters, 106(3):030802, 2011. Dmitry Budker and Michael Romalis. Optical magnetometry. Nature physics, 3(4):227{234, 2007. Jukka Sarvas. Basic mathematical and electromagnetic concepts of the biomagnetic inverse problem. Physics in Medicine & Biology, 32(1):11, 1987. Christoph M Michel, Micah M Murray, G oran Lantz, Sara Gonzalez, Laurent Spinelli, and Rolando Grave de Peralta. Eeg source imaging. Clinical neurophysiology, 115(10):2195{2222, 2004. Natalia A Trayanova, Bradley J Roth, and Lisa J Malden. The response of a spherical heart to a uniform electric eld: a bidomain analysis of cardiac stimulation. IEEE transactions on biomedical engineering, 40(9):899{908, 1993. Dan Xu and Bradley J Roth. The magnetic eld produced by the heart and its in uence on mri. Mathematical Problems in Engineering, 2017, 2017. Ryan A Murdick and BJ Roth. A comparative model of two mechanisms from which a magnetic eld arises in the heart. Journal of applied physics, 95(9):5116{5122, 2004. Bradley J Roth. A comparison of two boundary conditions used with the bidomain model of cardiac tissue. Annals of biomedical engineering, 19(6):669{678, 1991. Wanda Krassowska and John C Neu. E ective boundary conditions for syncytial tissues. IEEE transactions on biomedical engineering, 41(2):143{150, 1994. Bradley J Roth. Electrical conductivity values used with the bidomain model of cardiac tissue. IEEE Transactions on Biomedical Engineering, 44(4):326{328, 1997. Krista Kay McBride, Bradley J Roth, VY Sidorov, John P Wikswo, and Franz J Baudenbacher. Measurements of trans- membrane potential and magnetic eld at the apex of the heart. Biophysical journal, 99(10):3113{3118, 2010. Jenny R Holzer, Luis E Fong, Veniamin Y Sidorov, John P Wikswo Jr, and Franz Baudenbacher. High resolution magnetic images of planar wave fronts reveal bidomain properties of cardiac tissue. Biophysical journal, 87(6):4326{4332, 2004. Roger B Proksch, TE Sch a er, BM Moskowitz, ED Dahlberg, Dennis A Bazylinski, and Richard B Frankel. Magnetic force microscopy of the submicron magnetic assembly in a magnetotactic bacterium. Applied Physics Letters, 66(19):2582{2584, Lisa Qian, Beena Kalisky, Amanda Hamilton, Bo Dwyer, AC Matin, and Kathryn Moler. Magnetic characterization of individual magnetotactic bacteria. In APS Meeting Abstracts, 2011. Karen P Lam, Adam P Hitchcock, Martin Obst, John R Lawrence, George DW Swerhone, Gary G Leppard, Tolek Tyliszczak, Chithra Karunakaran, Jian Wang, Konstantin Kaznatcheev, et al. Characterizing magnetism of individual magnetosomes by x-ray magnetic circular dichroism in a scanning transmission x-ray microscope. Chemical Geology, 270(1- 4):110{116, 2010. Rafal E Dunin-Borkowski, Martha R McCartney, Richard B Frankel, Dennis A Bazylinski, Mihaly Posfai, and Peter R Buseck. Magnetic microstructure of magnetotactic bacteria by electron holography. Science, 282(5395):1868{1870, 1998. Arash Komeili. Molecular mechanisms of compartmentalization and biomineralization in magnetotactic bacteria. FEMS microbiology reviews, 36(1):232{255, 2012. Damien Faivre and Dirk Schuler. Magnetotactic bacteria and magnetosomes. Chemical Reviews, 108(11):4875{4898, 2008. A Krichevsky, MJ Smith, LJ Whitman, MB Johnson, TW Clinton, LL Perry, BM Applegate, K O'Connor, and LN Csonka. Trapping motile magnetotactic bacteria with a magnetic recording head. Journal of applied physics, 101(1):014701, 2007. Bruce M Moskowitz, Richard B Frankel, and Dennis A Bazylinski. Rock magnetic criteria for the detection of biogenic magnetite. Earth and Planetary Science Letters, 120(3-4):283{300, 1993. Mih aly P osfai and Rafal E Dunin-Borkowski. Magnetic nanocrystals in organisms. Elements, 5(4):235{240, 2009. Stephan HK Eder, Herv e Cadiou, Airina Muhamad, Peter A McNaughton, Joseph L Kirschvink, and Michael Winklhofer. Magnetic characterization of isolated candidate vertebrate magnetoreceptor cells. Proceedings of the National Academy of Sciences, 109(30):12022{12027, 2012. Nilesh R Ghugre and John C Wood. Relaxivity-iron calibration in hepatic iron overload: probing underlying biophysical 24 mechanisms using a monte carlo model. Magnetic resonance in medicine, 65(3):837{847, 2011. JAC Nicol. The giant nerve- bres in the central nervous system of myxicola (polychaeta, sabellidae). Quarterly Journal of Microscopical Science, 3(5):1{45, 1948. Yuyu Song and Scott T Brady. Analysis of microtubules in isolated axoplasm from the squid giant axon. In Methods in cell biology, volume 115, pages 125{137. Elsevier, 2013. KR Swinney and JP Wikswo Jr. A calculation of the magnetic eld of a nerve action potential. Biophysical journal, 32(2):719{731, 1980. Bradley J Roth and John P Wikswo Jr. The magnetic eld of a single axon. a comparison of theory and experiment. Biophysical journal, 48(1):93{109, 1985. John P Wikswo Jr and Bradley J Roth. Magnetic determination of the spatial extent of a single cortical current source: A theoretical analysis. Electroencephalography and clinical neurophysiology, 69(3):266{276, 1988. ML Markham, JM Dodson, GA Scarsbrook, DJ Twitchen, G Balasubramanian, F Jelezko, and J Wrachtrup. Cvd diamond for spintronics. Diamond and Related Materials, 20(2):134{139, 2011. Jagdish Narayan and Anagh Bhaumik. Novel synthesis and properties of pure and nv-doped nanodiamonds and other nanostructures. Materials Research Letters, 5(4):242{250, 2017. Vadym N Mochalin, Olga Shenderova, Dean Ho, and Yury Gogotsi. The properties and applications of nanodiamonds. Nature nanotechnology, 7(1):11, 2012. Amanda M Schrand, Houjin Huang, Cataleya Carlson, John J Schlager, Eiji Osawa, Saber M Hussain, and Liming Dai. Are diamond nanoparticles cytotoxic? The journal of physical chemistry B, 111(1):2{7, 2007. Felix Neugart, Andrea Zappe, Fedor Jelezko, C Tietz, Jean Paul Boudou, Anke Krueger, and J org Wrachtrup. Dynamics of diamond nanoparticles in solution and cells. Nano letters, 7(12):3588{3591, 2007. Liam P McGuinness, Yuling Yan, Alastair Stacey, David A Simpson, Liam T Hall, Dougal Maclaurin, Steven Prawer, P Mulvaney, J Wrachtrup, F Caruso, et al. Quantum measurement and orientation tracking of uorescent nanodiamonds inside living cells. Nature nanotechnology, 6(6):358, 2011. Wesley Wei-Wen Hsiao, Yuen Yung Hui, Pei-Chang Tsai, and Huan-Cheng Chang. Fluorescent nanodiamond: a versatile tool for long-term cell tracking, super-resolution imaging, and nanoscale temperature sensing. Accounts of chemical research, 49(3):400{407, 2016. Yuan Yuan, Xiang Wang, Guang Jia, Jia-Hui Liu, Tiancheng Wang, Yiqun Gu, Sheng-Tao Yang, Sen Zhen, Haifang Wang, and Yuanfang Liu. Pulmonary toxicity and translocation of nanodiamonds in mice. Diamond and Related Materials, 19(4):291{299, 2010. Nitin Mohan, Chao-Sheng Chen, Hsiao-Han Hsieh, Yi-Chun Wu, and Huan-Cheng Chang. In vivo imaging and toxicity assessments of uorescent nanodiamonds in caenorhabditis elegans. Nano letters, 10(9):3692{3699, 2010. Edward K Chow, Xue-Qing Zhang, Mark Chen, Robert Lam, Erik Robinson, Houjin Huang, Daniel Scha er, Eiji Osawa, Andrei Goga, and Dean Ho. Nanodiamond therapeutic delivery agents mediate enhanced chemoresistant tumor treatment. Science translational medicine, 3(73):73ra21{73ra21, 2011. Susana Q Lima and Gero Miesenb ock. Remote control of behavior through genetically targeted photostimulation of neurons. Cell, 121(1):141{152, 2005. Edward S Boyden, Feng Zhang, Ernst Bamberg, Georg Nagel, and Karl Deisseroth. Millisecond-timescale, genetically targeted optical control of neural activity. Nature neuroscience, 8(9):1263{1268, 2005. Jacob G Bernstein, Paul A Garrity, and Edward S Boyden. Optogenetics and thermogenetics: technologies for controlling the activity of targeted cells within intact neural circuits. Current opinion in neurobiology, 22(1):61{71, 2012. Daniel E Bath, John R Stowers, Dorothea H ormann, Andreas Poehlmann, Barry J Dickson, and Andrew D Straw. Flymad: rapid thermogenetic control of neuronal activity in freely walking drosophila. Nature methods, 11(7):756, 2014. Eric A Newman and Peter H Hartline. The infrared" vision" of snakes. Scienti c American, 246(3):116{127, 1982. Fumika N Hamada, Mark Rosenzweig, Kyeongjin Kang, Stefan R Pulver, Alfredo Ghezzi, Timothy J Jegla, and Paul A Garrity. An internal thermal sensor controlling temperature preference in drosophila. Nature, 454(7201):217{220, 2008. Shijia Chen, Cindy N Chiu, Kimberly L McArthur, Joseph R Fetcho, and David A Prober. Trp channel mediated neuronal activation and ablation in freely behaving zebra sh. Nature methods, 13(2):147{150, 2016. Gopalakrishnan Balasubramanian, Philipp Neumann, Daniel Twitchen, Matthew Markham, Roman Kolesov, Norikazu Mizuochi, Junichi Isoya, Jocelyn Achard, Johannes Beck, Julia Tissler, et al. Ultralong spin coherence time in isotopically engineered diamond. Nature materials, 8(5):383{387, 2009. Jocelyn Achard, Vincent Jacques, and Alexandre Tallaire. Cvd diamond single crystals with nv centres: a review of material synthesis and technology for quantum sensing applications. Journal of Physics D: Applied Physics, 2020. Ettore Bernardi, Richard Nelz, Selda Sonusen, and Elke Neu. Nanoscale sensing using point defects in single-crystal diamond: recent progress on nitrogen vacancy center-based sensors. Crystals, 7(5):124, 2017. Edlyn V Levine, Matthew J Turner, Pauli Kehayias, Connor A Hart, Nicholas Langellier, Raisa Trubko, David R Glenn, Roger R Fu, and Ronald L Walsworth. Principles and techniques of the quantum diamond microscope. Nanophotonics, 8(11):1945{1973, 2019. A Dr eau, M Lesik, L Rondin, P Spinicelli, O Arcizet, J-F Roch, and V Jacques. Avoiding power broadening in optically detected magnetic resonance of single nv defects for enhanced dc magnetic eld sensitivity. Physical Review B, 84(19):195204, Chang S Shin, Claudia E Avalos, Mark C Butler, David R Trease, Scott J Seltzer, J Peter Mustonen, Daniel J Kennedy, Victor M Acosta, Dmitry Budker, Alexander Pines, et al. Room-temperature operation of a radiofrequency diamond magnetometer near the shot-noise limit. Journal of Applied Physics, 112(12):124519, 2012. 25 Erik Bauch, Swati Singh, Junghyun Lee, Connor A Hart, Jennifer M Schloss, Matthew J Turner, John F Barry, Linh Pham, Nir Bar-Gill, Susanne F Yelin, et al. Decoherence of dipolar spin ensembles in diamond. arXiv preprint arXiv:1904.08763, Michael Sean Grinolds, Sungkun Hong, Patrick Maletinsky, Lan Luan, Mikhail D Lukin, Ronald Lee Walsworth, and Amir Yacoby. Nanoscale magnetic imaging of a single electron spin under ambient conditions. Nature Physics, 9(4):215{219, G De Lange, ZH Wang, D Riste, VV Dobrovitski, and R Hanson. Universal dynamical decoupling of a single solid-state spin from a spin bath. Science, 330(6000):60{63, 2010. Piotr Szank  owski, Guy Ramon, Jan Krzywda, Damian Kwiatkowski, et al. Environmental noise spectroscopy with qubits subjected to dynamical decoupling. Journal of Physics: Condensed Matter, 29(33):333001, 2017. Boris Naydenov, Florian Dolde, Liam T Hall, Chang Shin, Helmut Fedder, Lloyd CL Hollenberg, Fedor Jelezko, and J org Wrachtrup. Dynamical decoupling of a single-electron spin at room temperature. Physical Review B, 83(8):081201, 2011. Helena S Knowles, Dhiren M Kara, and Mete Atatu re. Observing bulk diamond spin coherence in high-purity nanodiamonds. Nature materials, 13(1):21{25, 2014. Q-Y Cao, P-C Yang, M-S Gong, M Yu, A Retzker, Martin B Plenio, C Muller,  N Tomek, B Naydenov, LP McGuinness, et al. Protecting quantum spin coherence of nanodiamonds in living cells. Physical Review Applied, 13(2):024021, 2020. Mariusz Mr ozek, Daniel Rudnicki, Pauli Kehayias, Andrey Jarmola, Dmitry Budker, and Wojciech Gawlik. Longitudinal spin relaxation in nitrogen-vacancy ensembles in diamond. EPJ Quantum Technology, 2(1):22, 2015. Mariusz Radtke, Ettore Bernardi, Abdallah Slablab, Richard Nelz, and Elke Neu. Nanoscale sensing based on nitrogen vacancy centers in single crystal diamond and nanodiamonds: achievements and challenges. Nano Futures, 3(4):042004, Nir Bar-Gill, Linh My Pham, Chinmay Belthangady, David Le Sage, Paola Cappellaro, JR Maze, Mikhail D Lukin, Amir Yacoby, and Ronald Walsworth. Suppression of spin-bath dynamics for improved coherence of multi-spin-qubit systems. Nature communications, 3(1):1{6, 2012. Analia Zwick, Gonzalo A Alvarez, and Gershon Kurizki. Maximizing information on the environment by dynamically controlled qubit probes. Physical Review Applied, 5(1):014007, 2016. Matthias M Muller,  Stefano Gherardini, Nicola Dalla Pozza, and Filippo Caruso. Noise sensing via stochastic quantum zeno. Physics Letters A, 384(13):126244, 2020. Christian Arenz, Daniel Burgarth, Paolo Facchi, Vittorio Giovannetti, Hiromichi Nakazato, Saverio Pascazio, and Kazuya Yuasa. Universal control induced by noise. Physical Review A, 93(6):062308, 2016. Ruslan I Stanika, Idalis Villanueva, Galina Kazanina, S Brian Andrews, and Natalia B Pivovarova. Comparative impact of voltage-gated calcium channels and nmda receptors on mitochondria-mediated neuronal injury. Journal of Neuroscience, 32(19):6642{6650, 2012. Juliette Humeau, Jos e Manuel Bravo-San Pedro, Ilio Vitale, Lucia Nunez, ~ Carlos Villalobos, Guido Kroemer, and Laura Senovilla. Calcium signaling and cell cycle: progression or death. Cell calcium, 70:3{15, 2018. http://www.deepdyve.com/assets/images/DeepDyve-Logo-lg.png Quantum Physics arXiv (Cornell University)

Is a quantum biosensing revolution approaching?

Download PDF
25 pages

Loading next page...
 
/lp/arxiv-cornell-university/is-a-quantum-biosensing-revolution-approaching-ePCm76MZ0F
ISSN
2511-9044
eISSN
ARCH-3342
DOI
10.1002/qute.202000066
Publisher site
See Article on Publisher Site

Abstract

Is biosensing revolution approaching? Review on Biocompatible ODMR Techniques 1;2 2 2 2 3 3 2 2;4 G. Petrini , E. Moreva , E. Bernardi , P. Traina , G. Tomagra , V. Carabelli ,I. P. Degiovanni , M. Genovese . Physics Department, University of Torino, Torino, Italy Istituto Nazionale di Ricerca Metrologica (INRIM), Strada delle cacce 91, Torino, Italy Department of Drug and Science Technology and NIS Inter-departmental Centre Torino, University of Torino, Italy Istituto Nazionale di Fisica Nucleare (INFN) Sez. Torino, Torino, Italy Understanding the human brain is one of the most signi cant challenges of the 21st century. As theoretical studies continue to improve the description of the complex mechanisms that regu- late biological processes, in parallel numerous experiments are conducted to enrich or verify these theoretical predictions also with the aim of extrapolating more accurate models. In the eld of mag- netometry, among the various sensors proposed for biological application, NV centers are emerging as a promising solution due to their perfect biocompatibility and the possibility of being positioned in close proximity to the cell membrane, thus allowing a nanometric spatial resolution down to the nano-scale. Still many issues must be overcome to obtain both spatial resolution and sensitivity capable of revealing the very weak electromagnetic elds generated by neurons, or other excitable cells, during their ring activity. However, over the last few years, signi cant improvements have been achieved in this direction, thanks to the use of innovative techniques. In this review, the new results regarding the application of NV centers will be analyzed and the main challenges that must be a orded for leading to practical applications will be discussed. PACS numbers: I. INTRODUCTION Electromagnetic eld sensing is of the utmost importance for several applications in current scienti c research, fostering the search for novel high sensitivity sensors. Several innovative electromagnetic eld sensors emerged in 1{10 the last years , whose main goal is revealing less and less intense elds with an increased spatial resolution. In particular high sensitivity sensing coupled to high resolution is of the utmost relevance in biological research, especially, for instance, in studies of human brain cell currents, which are typically extremely faint. Localized 11{14 monitoring of neuronal elds would allow not only the investigation of brain currents during cognitive processes in order to improve neurological diagnostic systems, but also identifying the early stages of neurodegenerative disease, 15,16 like Parkinson's, Alzheimers's disease and other forms of dementia . Furthermore, since localized temperature gradients and heat dissipation occur within cellular microdomains, the 17,18 exploitation of sensing probes for multi-task applications would be very fruitful . Among the various sensing 19,20 devices that have emerged over the years , promising sensors for the detection of biological elds are color 21,22 23{30 centers in diamond . Color centers are impurities in the crystalline matrix that, when stimulated, emit uorescence. In particular, the nitrogen-vacancy (NV) complex is by far the most promising due to its level structure . This dependence allows the realization of techniques for optical initialization and spin readout by means 33 34 of the Optically Detected Magnetic Resonance (ODMR) technique . Furthermore, its spin energy levels are 35{41 42 sensitive not only to electromagnetic elds , but also to temperature variations . These exceptional properties, together with their photostability at room-temperature and the non-toxicity of diamond , promotes the NV complex 44{49 as a very promising candidate for biological application . The paper is structured as follows: section 2 reports a brief description of the theory of quantum sensing, section 3 analyzes the magnetic eld generated by mammalian neuronal cells and cardiac tissue, section 4 deals with experiments aimed at the detection of cells elds and nally in section 5 the experimental techniques, used to enhance the sensitivity of NV centers to be used as biosensors, are highlighted. II. THE THEORY OF QUANTUM SENSING WITH NV CENTERS The nitrogen-vacancy (NV) defect is a natural complex of impurities in diamond crystalline matrix. This complex is composed of a substitution nitrogen atom and a vacancy-type defect, located in adjacent reticular sites . This system has a pyramidal symmetry (C ) and it has, as axis of symmetry, the line that connects the nitrogen atom 3v with the vacancy (see Figure 1a). With respect to the tetrahedral structure of the diamond, there are 4 possible arXiv:2006.00831v2 [physics.app-ph] 1 Jul 2020 2 orientations of this defect, all equiprobable in conditions of conventional syntheses. Moreover, there are two charged states in which it is possible to nd the nitrogen-vacancy defects, distinguished by the number of electrons involved. The 3 carbon atoms surrounding the vacancy contribute to sharing 1 electron each to the complex, while nitrogen contributes with 2. If, in total, only these 5 electrons are present in the system, the center is electrically neutral and it is referred to as NV , with total electronic spin S = 1/2. Alternatively, the defect can trap 1 additional electron from the surrounding lattice, creating the NV center. In this case the total electronic spin becomes S = 1, with spin component along symmetry axis fjm = 0 >;jm = +1 >;jm = 1 >g. The most promising con guration for s s s quantum sensing exploits the spin property of the NV complex. Its sp orbitals linearly combine to form 4 molecular orbitals: the lowest energy state of the ground con guration that is the orbital singlet, spin triplet state A , the 3 1 1 electronic excited states that are orbital doublet, spin triplet E, spin singlet orbital singlet E and A . a) b) FIG. 1: a) Diamond crystalline structure with nitrogen-vacancy defect; b) NV radiative state transitions that occur during 3 3 1 laser pumping. Radiative optical transition E ! A with 637 nm zero phonon line (ZPL), and non optical transition E ! 1 3 1 1 A with 1042 nm ZPL. Non radiative intersystem crossing (ISC) transitions subsist between E ! A and between E ! 1 1 A . As we will discuss, due to its electronic levels the NV complex is the most promising con guration for quantum sensing. By irradiating the complex with a 532 nm pump laser (Figure 1b), the electronic state is excited in a non-resonant way and afterwards it relaxes to the fundamental state emitting a photon with wave lenght between 637 nm (zero phonon line) and 800 nm (phonon sideband). While the optical excitation from the jm = 0 > state is spin preserving, the transition from jm = 1 > has a nite branching ratio into the metastable singlet E, with a lifetime of 300 ns. This singlet state relaxes into jm = 0 > through a non-radiative processes and weak infrared emission peaking at 1042 nm, leading to a drop in uorescence output up to 30% for a single NV , or 1-2% for a large NV ensemble, with respect to the situation when the system is initialized in jm = 0 >. This allows the optical readout of the spin state. Since the main focus of the present work is on sensing, from now on we will refer to NV as NV for simplicity. A. NV ground electronic state 3 34,50 The Hamiltonian of A , the ground spin state of the NV system, can be written in the following form : gs ^ ^ ^ ^ ^ ^ = SDS + SAI + IQI (1) ^ ^ ^ ^ ^ ^ ^ ^ where S = (S ; S ; S ) and I = (I ; I ; I ) are the dimensionless electron and nitrogen nuclear spin operators, x y z x y z respectively. The rst term represents the ne structure splitting due to the electronic spin-spin interaction, coupled by the ne structure tensor D. The second term is generated by the hyper ne interaction between NV electrons and 14 15 the nitrogen nucleus (I = 1 for a N nucleus, while I = 1/2 for a N nucleus), with the hyper ne tensor A. Finally, the third term represents the nuclear electric quadrupole interaction, with the electric quadrupole tensor Q. It should be noted that, in this notation, the component z coincides with the NV axis of symmetry. Due to the symmetry of 51,52 the NV center, D, A, and Q are diagonal in the NV coordinate system and, in terms of the natural spin-triplet 3 basis fjm = 0 >;jm = +1 >;jm = 1 >g, the Hamiltonian can be written as: s s s gs 2 2 == ? 2 2 ^ ^ ^ ^ ^ ^ ^ ^ ^ ^ = D [S S =3] + A S I + A [S I + S I ] + Q [I I =3] (2) gs z z x x y y gs z gs gs z h | {z } | {z } | {z } electronic spinspin interaction nuclear spinspin interaction electronnucleus spin interaction == where D ' 2.87 GHz is the zero eld splitting, Q is the nuclear electric quadrupole parameter, A and A are gs gs gs gs 53,54 the axial and non-axial magnetic hyper ne parameters . The parameters values are reported in Table I. Hyper ne parameters Value Zero eld splitting D ' 2.87 GHz gs == Axial hyper ne term A ' -2.14 MHz gs; N == A ' 3.03 MHz gs; N Transverse hyper ne term A ' -2.70 MHz gs; N A ' 3.65 MHz gs; N Nuclear electric quadrupole term Q ' -5 MHz gs TABLE I: Hyper ne parameters for the NV defect determined at room temperature. B. The Optically Detected Magnetic Resonance technique One of the characteristics that makes NV centers so attractive and convenient as key element in various type of sensors is the possibility to discriminate the spin components of the electronic state. This is allowed by the di erent coupling of the jm = 0 > state with a metastable level, compared to the jm = 1 > states and results s s in a variation of the photoluminescence (PL) of the defect under laser non resonant excitation. Optically Detected Magnetic Resonance (ODMR) consists in the application of a microwave eld (MW) on the sample, simultaneously with its exposure to a non-resonant laser at a frequency higher than the resonant one, corresponding to the energy gap between the ground and the E level (e.g. 532 nm) (see Figure 2a). When the frequency of the MW reaches the ground state resonance D of the NVs, with a certain probability (depending on the MW power), those NV centers gs will be initialized in the states jm = 1 > rather than jm = 0 >. As mentioned, this corresponds to a reduction in s s photoluminescence of the NV centers, as it can be observed, e.g., in Figure 2b, where a typical ODMR spectrum is reported with the expected uorescence dip at the zero eld splitting frequency D . gs a) b) FIG. 2: a) NV radiative state transitions that occur during laser pumping and MW excitation. The coupling of the state jm = 1 > with the metastable level generates a statistically lower uorescent emission than the jm = 0 > initialized s s electronic state. b) Fluorescence collected from the NV center as a function of the MW frequency varies. A dip in correspondence of the zero eld splitting D (resonance frequency of the undisturbed NV center, at room temperature) can be observed. gs The coupling terms of a NV center with the electric, magnetic elds and local temperature variations will be analyzed in following. 4 C. Magnetic eld sensing A static magnetic eld produces the well-known Zeeman e ect , that it is described by: == V  g  g gs B gs B gs N N ^ ^ ^ ^~ = S B + (S B + S B ) + IB (3) z z x x y y h h h h | {z } | {z } Zeeman interaction nuclear Zeeman interaction == where  is the Bohr magneton,  is the nuclear magneton, g and g are the components of the ground state B N gs gs electronic g-factor tensor and g is the isotropic nuclear g-factor. In the presence of relatively weak magnetic elds, it is possible to approximate the almost diagonal g-factor tensor in a diagonal form, with constant g = 2.003 . As reported in Table II, the interaction of the magnetic eld with the nucleus is 2000 times smaller and, consequently, it is typically neglected . The presence of external elds eliminates the energy degeneracy of the levels jm = 1 >, B e whose splitting become B , where = (see Figure 3). e z e a) b) 3 14 15 FIG. 3: NV ground-state A scheme. Above: a) N hyper ne states and b) N hyper ne states. Below: schematic ODMR spectra. The spectra are shown considering Zeeman splitting and hyper ne splitting. If, instead of a single NV center, an ensemble of NV centers is considered, up to eight magnetic resonance dips can be observed, due to the four possible orientation of the NV axis in the diamond's crystalline matrix (see Figure 4). For certain directions of the magnetic eld, some resonances can be degenerate. A NV-based magnetometer can be realized, for example, by applying a bias eld along the NV axis, removing the degeneracy, so that changes in the magnetic eld projection along this axis a ect the resonance frequencies almost linearly. Another option is to use all four NV alignments; although the eight ODMR frequencies have more complicated dependence on B, this option yields information about the direction of magnetic eld . 38,39 35,40 The use of NV center as a magnetic eld sensor rstly was proposed in and demonstrated with single NV and NV ensembles in 2008. D. Electric eld sensing The Hamiltonian describing the interaction with the electric eld was derived from molecular orbit theory by Doherty et al. and it can be written in the following form: ^ ^ V S gs == 2 ? 2 2 ? ^ ^ ^ ^ ^ ^ ^ = d (E + F )[S ] + d (E + F )(S S ) + d (E + F )(S S + S S ) (4) z z x x y y x y y x gs z gs y x gs h 3 5 B≠0 B=0 |0> |-1> |+1> |±1> b) MW [GHz] a) MW [GHz] B≠0 c) MW [GHz] FIG. 4: ODMR spectra in the absence a) and in the presence b) of an external bias magnetic eld. The magnetic eld lifts the degeneracy of the jm = 1 > states and results in two separate dips in the ODMR spectrum. c) An example ODMR spectrum (excited at 532 nm) with a magnetic eld in an arbitrary direction for an ensemble NV centers in diamond. Each of the four NV alignments has a di erent magnetic eld projection along its quantization axis, leading to eight ODMR peaks (two for each NV alignment). For each dip a coupling with the nuclear spin of the N atom generates additional three hyper ne levels . == where d and d are respectively the axial and non-axial Stark shift components of the permanent electric dipole gs gs ? 57 ~ ~ moment d in the ground triplet state , E is the electric eld and F is the mechanical strain. gs 41,58 ~ ~ According to Equation (4) the e ect of the electric eld E plays the same role as mechanical strain F . The strain depends on the diamond material: in single-crystal samples, the mechanical strain eld is substantially negligible; while, in polycrystalline ones, a relatively high strain eld is induced by the growth conditions, leading to a splitting of the spin state jm = 1 > even in absence of external elds. The frequency shift caused by the electric eld is much smaller than the shift produced by the presence of a magnetic eld (see Table II). For this reason, in order to reliably measure this second-order e ect caused by the Stark shift, it is necessary to decouple it from the Zeeman shift. Brie y, the ne structure Hamiltonian of the NV ground state, describing the energy levels of the electronic spin ^ ~ ~ ~ states due to the spin (S) interaction with the static magnetic (B), electric (E), and strain (F ) elds, can be written in terms of the natural spin-triplet basis fjm = 0 > , jm = +1 >, jm = 1 >g in the following matrix form: s s s 0 1 B iB B +iB x y x y p p 0  g  g B e B e 2 2 B C B +iB x y ? ^ B C H =  g hD +  g B hd (P iP ) gs B e B e z x y gs @ 2 A B iB x y g hd (P + iP ) hD  g B B e x y B e z gs where it is possible to observe that the natural-spin basis vectors are eigenstates of the Hamiltonian only in the == presence of both the magnetic and electric eld aligned with the NV axis. In this condition, D = D + d P gs gs z describes the frequency shift of the resonance lines resulting from the zero- eld splitting and from the Stark e ect ~ ~ ~ associated with the component of the vector P = E +F . Otherwise, external elds not aligned to NV symmetry axis produce a non-diagonal matrix, and therefore energy levels of unde ned spin. In particular, the presence of additional transverse strain and electric- eld components P modi es the ground-state structure. The Hamiltonian assumes a quasidiagonal form considering a new spin basis fj0 >, j+ >, j >g, obtained by a eld-dependent mixing of the jm = +1 > and jm = 1 > spin states according to the following unitary operator: s s 0 1 1 0 0 B   C ^ i  i 2 2 U = 0 e sin( ) e sin( ) @ A 2 2 i i 2 2 0 e cos( ) e sin( ) 2 2 Fluorescence [a.u] Fluorescence [a.u] Fluorescence [a.u] 6 ? ? where tan() = P =P and tan() = (d P )=( g B ) are the eld-dependent phases de ning the spin state mixing. x y B e z gs The Hamiltonian takes the following form in the fj0 >,j+ >,j >g basis: 0 1 ? ? 0 c  g B c  g B 1 B e 2 B e B C ^0 ^ ^ y  ? H = UH U = gs gs @c  g B hD + W 0 A B e c  g B 0 hD W B e with ? ? 2 2 W = (hd P ) + ( g B ) (5) B e z gs The complex constants c and c represent the phase of the matrix elements and B is the transverse component 1 2 of the magnetic eld with respect to the NV axis. If B  0, the non-diagonal terms can be neglected and the Hamiltonian can be regarded as diagonal in the basis fj0 >, j+ >, j >g. The energy di erence between the j0 > and the j > states is hD  W , corresponding to ODMR resonances separated by 2W/h depending on the strengths of the magnetic, electric, and strain elds, as well as their orientations with respect to the axes of the NV center. Since the states j > are a coherent superposition of the states jm = 1 >, we underline that the ODMR resonance is observed also in this case as a reduction in the uorescence emission at the new MW resonance frequencies. E. Temperature sensing Another interesting feature of the NV complex is the temperature dependence of its spin levels . Indeed, the microscopic origin of D , also called the zero eld splitting (ZFS) parameter, is due to spin-spin interactions in gs the NV's orbital structures, and the value depends on the lattice length, which is strongly correlated to the local temperature. When the local temperature increases the diamond lattice spacing of the NV center increases as well, lowering the spin-spin interaction and reducing the ZFS parameter D . Under ambient conditions D ' 2:87 GHz gs gs and the temperature dependence is dD=dT ' 74kHz=K . In general, the ZFS parameter shows a non-linear dependence, and its value increases when temperature decreases . Property Coupling coecient e 1 Magnetic eld = ' 28 GHz T N N = ' 15 MHz T Electric eld d ' 3.5 mHz V m ==;gs d ' 0.17 Hz V m ?;gs Temperature @D =@T ' -74 kHz K gs TABLE II: Coupling coecient of the NV center with the external elds and temperature. To realize a NV-based temperature sensor, the most obvious solution is exploiting the D temperature dependence. gs ~ ~ ~ This requires that no external eld is present (B; E; F = 0), i.e. jm = 1 > is degenerate. In this case, an increase in temperature leads to a decrease in the resonance frequency, associated with a shift of the degenerate levels jm = 1 > towards the level jm = 0 >. Nevertheless, this is the simplest but not the optimal solution, since, even in the absence of applied elds, the sample may have an internal strain and may be a ected by the Earth's magnetic eld. Unless it is possible to nd a diamond sample with negligible F and to design an experimental set up able to reasonably compensate for the external magnetic eld (e.g. Helmholtz coils), the dips would not be perfectly overlapped because of the non perfect degeneracy of jm = 1 >, thus showing a larger full-width-at-half-maximum (FWHM) and therefore a lower resolution. A better solution is to apply an external magnetic eld in order to signi cantly separate the spin levels. However in this con guration, a single dip can shift for a temperature variation, but also for a variation of magnetic eld. To decouple the two contributions it is sucient to monitor both jm = +1 > and jm = 1 > spin states at the same s s time, using simultaneous driving of the microwaves in ODMR technique . As it can be seen in the Figure 5, by simultaneously monitoring the initial dips (red curves), it is in principle possible to understand if there are variations in the magnetic eld (the dips move in opposite directions) or in temperature (the dips move in the same direction). In this perspective, an improved technique exploiting the application of an intermediate transverse bias magnetic ? 61 ? eld B has been implemented . Similarly to the case just discussed, the application of B removes the degeneration 7 FIG. 5: Example of magnetic and thermal shifts of the spin resonance, in ODMR spectra. Dips with equal colors correspond to paired resonances. The colors represent the timeline of the dips. The initial dips is red, then green and nally blue. of jm = 1 > and therefore improves the FWHM. The intensity and the transverse direction of that eld creates a quantum superposition of states, which is insensitive to magnetic elds but sensitive to temperature . In this con guration, the expectation value of the spin along any direction is small, implying the degeneracy of the hyper ne structure between the levels jm = 1 > (except for the quadrupole contribution Q , which seperates jm = 0 > I gs I from jm = 1 >). In Figure 6 the corresponding scheme of the spin energy levels are reported (only the N isotope is considered as it is the most common). In this situation, the ODMR spectrum reduces to two dips (instead of 6), providing a substantial improvement in the signal-to-noise ratio. This particular orientation of the magnetic eld ensures the protection of the measurements from the noise of other possible magnetic elds. In fact, the NV spin is non-sensitive to the magnetic eld uctuation, because the contribution of the magnetic component only appears at the second order in the Hamiltonian. 3 ? FIG. 6: NV ground-state A scheme, in presence of intense transverse magnetic eld B . III. BIO-SENSING As mentioned above, NV sensors are particuarly suitable for biological sensing. Before describing the experiments focusing on the NV-based sensor it is necessary to specify the type of biological specimens of interest, the expected magnitude of the electromagnetic eld produced by these specimens and the principal parameters such as sensitivity and spatial/temporal resolution, required for the NV-based sensors. This section, after reviewing some of the devices typically used for biosensing, analyzes in detail neuronal and cardiac cells. Higher sensitivity and resolution of electromagnetic elds is considered necessary to expand the understanding of the fundamental processes regulating the interaction of these cells. A. From the conventional electrophysiological techniques to NV sensors The electrical activity of excitable cells can be investigated by means of the conventional patch-clamp technique or through Micro Electrode Array (MEA) recordings. Single-cell recordings, performed under voltage- or current-clamp con guration, respectively allow to monitor ion currents or the membrane resting potential, postsynaptic responses and action potential ring activity . Besides having an extreme versatility (monitoring overall electrical events from 8 the whole cell, from microdomains of the cell membrane or even from single channel proteins), patch-clamp has a high temporal resolution and high sensitivity: all these features make this electrophysiological approach the gold standard for measuring the electrical activity. Though, patch-clamp is rather invasive, as it damages the cell membrane through the recording electrode: this implies that only one recording is feasible for each cell . On the contrary, MEA is a non-invasive approach, used to measure the membrane potential variations from many cells simultaneously. The MEA probe is structured as an array of sensing electrodes, of variable geometry dimensions and material, which are im- mersed in a glassy (insulating) double layer. Commonly, sensing electrodes are made of titanium or indium tin oxide 2,3 (ITO), and have a diameter that can vary between 10 and 30 m . By means of MEAs, it is possible to monitor the electrical activity of a neuronal network as a whole, and measuring its changes along with its maturation, even though informations on the biophysical properties of ion channel cannot be directly inferred. This speci c measurement need, requiring non-invasive and iterative detection for biological applications has prompted the study and realization of di erent devices. In the following we discuss and compare the most promising ones. Promising devices for the detection of weak magnetic elds, in addition to the NV-based sensors, are the super- 4{6 7 conducting quantum interference device (SQUID) sensors and chip-scale atomic magnetometers (CSAMs) . Until now, the measurement of very weak magnetic elds was the domain of SQUIDs sensors. These sensors have reached 1=2 2 5 sensitivity levels of (0.9{1.4) fT/Hz with a pick-up coil area of the order of 1 cm . However, SQUIDs require cryogenic cooling, which, in addition to implying signi cant cost and maintenance complexity, requires positioning the sensor a few centimeters from the sample. An alternative is o ered by the CSMAs, that are based on microfabricated alkali vapor cells integrated with small optical components such as diode lasers and ber optics. These devices have 1=2 37 reached sensitivities below 5 fT/Hz at sensor volume 8 mm . However, despite the exceptional sensitivity, the minimum working distance between sensor and magnetic source for CSAM or SQUIDS remains at least few mm, that makes them unsuitable for monitoring individual cell signals or small tissues, being the amplitude of the magnetic eld decreasing quadratically with distance. The sensors for the detection of electric elds, emerging in the last few decades are single-electron transistors (SETs) , that are a promising candidate for achieving higher detection sensitivity due to the Coulomb oscillations. However, the existence of a SET- based biosensor has emerged only in recent years , probably because of their diculty of the room-temperature operation. Finally, in recent years there has been a growing interest in the use of temperature sensors capable of operating on a nanometric scale. It has been known that local temperature variations at the intracellular level play a fundamental role in cellular activities related to body temperature homeostasis and energy balance . Particular attention is paid to the possibility of measuring local temperature variations of cell organelles (i.e. nucleus, mitochondria, etc.) or 67,68 ion channels. For example, di erent simulation model shows a hypothetical variation in temperature at the level of the ion channels, due to the ow of the ions from the inside to the outside of the plasma membrane, during the genesis of the action potential. Due to the diculty of this local measurement, no one has ever measured this thermal variation. Interestingly, temperature changes may drastically alter the neuronal ring frequency, as demonstrated by Guatteo et al. . Currently uorescence probes are powerful method used to study intracellular temperature variation thanks their high spatio-temporal resolution. The probes typically used for this measurement are organic or inorganic 70{73 9,74 uorescent probes, such as uorescent proteins, organic dyes , quantum dots (QDs) and many others. Organic proteins are biocompatible probes, rather stable and very easy to chemically target. But there are di erent problems related the use of these probes: these are often auto uorescent and to avoid the phenomenon it is necessary to add speci c quenchers; they cannot be used for a long time, in fact these sensors su er from photobleaching and unstable photoluminescence. In the best case scenario, the probe degradation consists of uorescence suppression, in the worst case scenario it releases an electron that binds to nearby molecules making them toxic. These probes are organic and by their nature they are also subject to even weak pH variations, for this reason it is fundamental a strict control 75,76 of the cell environment . The inorganic probes such as quantum dots (QDs) have the advantage of being stable in uorescence, have a high sensitivity to temperature variations and their nanometric size allows obtaining a spatial resolution useful for cellular measurements. Although the size of these sensors would allow spatial resolution limited by the di raction limit only, their chemical composition is found to be non-biocompatible in most of the cases. Other 10,66 temperature sensors are based on up converting nanoparticles (UCNPs) : nanoscale particles (diameter 1-100 nm) that exhibit photon upconversion, i.e. when stimulated by incident photons they are able to emit uorescence's of shorter wavelength. They are usually composed of rare-earth based lanthanide or actinide-doped transition metals. Their core-shell structure allows sensor compatibility, however, sensitivity is not high. Extremely interesting devices able to realize all these measurements (magnetic, electrical and temperature sensing) eventually at the same time, are one based on the NV center in diamond. The advantages of these sensors are man- ifold: they have stable photoluminescence in the visible and near-infrared range, their chemical composition ensures resistance to photobleaching and diamond is an inert and therefore biocompatible material , so cell/neurons can be 27,28,77 grown directly on its surface or nanodiamonds can be injected inside them, allowing for sub-cellular spatial resolutions with a non-invasive techniques. Finally, NV sensors can operate at room temperature and, in more de- 9 78 79 tail, their dynamical range of temperature sensing extends further 500 K for both bulk and nanoscale diamonds. In the remainder of this section, the modelization of the target neural and cardiac signals for NV-based biosensing will be reviewed. B. NV center as sensor for neuronal signals In the last decades, neuroscience has attracted great interest beyond the scienti c community. Because of the increase in life expectation, cases of neurodegenerative diseases such as Parkinson's, Alzheimer's, Huntington's disease and many others are constantly growing. Currently, these diseases are incurable, even symptoms mitigation is dicult because of late diagnosis when most of the neurons involved have been irreparably damaged. This reason strongly prompts to develop new increasingly precise and sensitive techniques, allowing a deeper understanding of neuronal circuits ranging from functioning of synaptic sites to the behavior of the entire neuronal network. Neurons are the functional units of the nervous system. They communicate via electrical signals, known as action potentials. The action potential (AP) consists in the variation in time of the membrane potential V , where V = e e m m in out is the electrical potential di erence between intra- and extra-cellular side of the cell membrane. The AP characteristic + + 2+ waveform is shown in Figure 7b. The AP pulse is caused by several ionic species (Na , K , Ca ), which ow through the neuronal membrane. a) b) FIG. 7: a) Single neuron simpli ed sketch. In the upper box a zoom of the neuronal membrane is reported, where the ionic current and the corresponding magnetic eld are schematized. In the lower box the axial current and the relative magnetic eld are shown. b) Schematic representation of neuronal action potential (AP). Resting membrane potential (-V ) is -70 mV. When V is driven and exceeds the threshold (following an initial stimulus), a rapid membrane depolarization occurs. In this phase the Na channels open, allowing sodium to enter in the neuron and bringing V to approximately +35 mV. Then the + + repolarization phase begins, caused by Na channels inactivation and opening of K channels. This outward current drives + + the membrane potential close to -93 mV (hyperpolarization). Finally, the Na /K ATP-ase restores the initial conditions. During the depolarization, the in ux of positive charges produces local internal and external longitudinal currents, which are responsible for the AP propagation in the axon adjacent area. The propagation directionality is guaranteed by the AP refractory period: although the local currents propagate in both directions, a new AP cannot be triggered in refractory membrane area. The two electrophysiological techniques mostly used to study cell excitability and synaptic transmission in a neuronal network are the patch-clamp and the MEA. In the last decade scientists have tried to study more and more speci cally the propagation of the electrical signal from the cell body (or soma) to the whole dendritic tree. In other words, the goal would be to create a device that allows scanning the neuron point by point from the soma to the axon and the dendrites, following and characterizing the electrophysiological variations of the electrical signal during its propagation. The technology closer to this ambitious goal is the one of the CMOS-MEA, that allows having a much higher density of electrodes with respect to the traditional MEA technology. Numerous stud- 80? {82 ies have managed to scan the path of the electrical signal in a neuronal network at the level of the single neuron . Bakkum et al. recently have developed a high electrode density CMOS-MEA device capable of stimulating a speci c area and simultaneously scanning the signal along some points from the soma to the axon. Clearly, this technique is much more sensitive than MEA, but given the stochasticity of the cell's placement in space, it requires cells to be marked in order to follow their path. Recently several groups have correlated this technology to the technique of optogenetics. They tagged the genes of interest and activated them following an optical stimulation and 10 83,84 simultaneously followed the signal thanks to the integration of the CMOS-MEA . However, these techniques do not allow following the entire dynamics of the action potential, but to have a scan of a region depending on the position of the electrodes with respect to the neuron with its axon and its dendritic body. NV sensors may therefore have a huge impact on these applications: nanodiamonds can be targeted on the membrane surface or, alternatively, cells can be plated and cultured on a bulk diamond . Indeed, taking advantage of diamond biocompatibility and the exceptional spatial resolution displayed by color centers in diamonds , it will be worth exploiting these properties for a timely reconstruction of the AP dynamics. Furthermore, the possibility of positioning them adjacent to the cell membrane has the advantage of experiencing stronger magnetic elds. However, since neuronal magnetic elds are extremely weak (' pT ), their detection appears to be challenging even for NV-based sensors, at least for mammalian cells, while measurements have been performed on giant neurons of invertebrates . To predict the electromagnetic elds intensity created by the AP, and therefore to understand what sensitivity of the NV sensors is needed to sense it, it is necessary to model how the AP develops and propagates. 11{14 Hodgkin-Huxley model allows estimating the ionic current owing through the neuron membrane (when the ion channels are open). For the human neuron, the total estimated ionic current, sum of the single channels contribution I is: ion I = I ' 2 pA m ? ion and the current pulse typically lasts t ' 1 ms. Each I generates a magnetic eld (see Figure 7a), which can ion ~ ~ be estimated by means of the Biot-Savart law: B  l =  I . However, the resulting amplitude of these elds ion 0 ion depends on the channels density, which largely varies depending on the axon area being considered. Furthermore, we tot note that the B eld, sum of the contributions of the eld produced by the various channels, can be vanishingly small ion 86,87 on average, because of the di erent elds directions. To this purpose, channel clustering may be very signi cant . Assuming a current of 100 pA=m and considering that the NV sensor positioned at an average distance of few 88,89 nanometers (by selective targeting the channel using functionalized NDs ), a magnetic eld of about 0:1 5 nT (or even higher) could probably be sensed. This hypothesis is now under experimental analysis . Current owing through the membrane is not limited to the charge ow through ion channels, as longitudinal currents, but one should also consider the ow along the neuron axis, that is responsible for the AP propagation. These currents also generate a magnetic eld, around the neuron (see Figure 7a). Both the axial current and 90{92 90 the corresponding magnetic eld have been estimated . In particular, Ref. goes beyond the simpli cation of the Hodgkin-Huxley model, introducing the spatial and temporal progression of the AP along the various neuronal compartments, into which they have divided the axon. The theoretical prediction is a maximum eld B ' 3 pT axial on the external membrane near the Ranvier node and a eld B ' 2.3 pT on the myelin sheath external surface axial in those regions where the axon is wrapped by it. The maximum magnetic eld was also calculated by Isakovic et al. in Ref. for the nerve composed of 100 axons, obtaining only B ' 6 pT. This is due to the cancellation of the magnetic eld component, caused by di erent axial axons within the same nerve, bringing opposite directional currents. This estimated magnetic elds, in reality, are compatible with the elds detected by magnetoencephalography (MEG). MEG is able to detect elds if the order of 15 93 10 T because of the distance from the source . Considering these values, a NV sensor positioned on the neuron surface or a few micrometers from it, should have a temporal resolution of about 0.1 ms (in order to be able to trace the time variation), and spatial resolution of about 10 m (which would allows a good reconstruction of the AP propagation, being the axon length ranging from 0.1 m to 1 m). Thus, the NV sensor should have a minimum sensitivity of : p p 1=2 = B t ' 3 pT 0:1 ms ' 30 fT Hz (6) min The NV sensor optimal sensitivity is in principle limited by the quantum projection noise. This fundamental sensitivity limit for spin-based magnetometers is given by : 1 1 = (7) nT Where is the magnetic coupling coecient (Table I), n represents the number of NV centers and T their characteristic dephasing time. It is important to underline that the number of NV centers n refers to the sensing volume. As mentioned, for the single PA detection the sensing volume should be around 10 m , the size of the cell. 77 6 3 In the Ref. , the estimation of the parameters n ' 3 10 cm and T ' 450 ns determines a spin projection noise 1=2 3 value of  ' 30 pT Hz for the sensing volume of 10 m (the experimental sensitivity reached is instead  ' 15 q 11 1=2 6 3 pT Hz for the sensing volume of 5  10 m ). This value is still 1000 times larger than the sensitivity required for the detection of a single AP. However, as will be discussed in section 5, it is possible to optimize both the above mentioned parameters to improve the performances. Once the biomagnetic eld B (x~ ; t) has been measured, to reconstruct the unknown currents generating it, one should solve an inversion problem. In general, its solution is not unique, due to the existence of the so-called "magnetically silent" currents (i.e. the ones producing magnetic elds that almost cancel each others) and due to 96,97 the fact that the magnetic eld can be in uenced by the electric eld . However, in the single axon case, it can be uniquely resolved. On the contrary, in the biological tissue case and in the 3D structures case, that cannot be traced back to standard models (such as a spherically symmetrical conductor or a horizontally layered medium), the solution is not unique. In some cases this is resolved by the knowledge of the electric eld on the conductor surface . C. NV center as sensor for cardiac signals The human (and animal) heart generates the body's most intense electromagnetic eld. In particular, by comparing measurements performed externally to the human body, the electric eld generated by the heart, measured through the electrocardiogram (ECG) is about 60 times stronger than that the one of the brain, recorded by an electroen- cephalogram (EEG). In addition, the heart magnetic eld detected by the magnetocardiogram (MCG) is about 5000 times higher than the neuronal magnetic eld detected by magnetoencephalography (MEG): 0,05 nT (heart) vs 1 fT (neuron). Thus, ODMR based on NV sensors can also nd very signi cant applications in studying cardiac cells and tissues. To achieve a rst qualitative estimation of the magnitude of the magnetic eld in this case, one can start from a very simpli ed model: the spherical heart . Although this model is not physiologically accurate, it allows to extrapolate analytical solutions. In a more recent work , is proposed a further assumption concerning the origin of the currents. There are two currents sources in the heart: the rst consists of intracellular currents, the second is given by the anisotropy of the tissue . Regarding the rst current contribution, the authors consider a spherical shell of cardiac tissue, which covers a blood cavity and is surrounded by an external bath of unlimited electrical conduction. The heart bers propagate in the z direction and a variation of the membrane potential V is assumed following the activation of the action potential (AP), started at  = 90 (see Figure 8a). In this work the electric eld is evaluated using the bidomain model and considering a situation of quasi-stationarity (although V depends on time due to the action potential propagation, it is assumed that, given a certain V (t ), m m 0 one can derive current and magnetic eld in a quasistatic way). 101,102 Thus, the electric potential is obtained, using the continuity equations and the boundary conditions , the current density distribution is obtained using Ohm's law and nally the magnetic eld using Biot-Savart's law. Considering 103 98 the anisotropic electrical conductance data , the V values and typical heart dimensions , it turns out that the magnetic eld is stronger near the internal and external surfaces tissue while it is weaker in the heart wall. The peak value of the magnetic eld is around 14 nT (see Figure 8b). At the heart center, instead, the magnetic eld reduces to B = 2 nT . This is due to the fact that intracellular and extracellular currents are in opposite directions with almost the same magnitudes in the depths of the tissue and, therefore, the corresponding magnetic elds essentially cancel each other. Considering a planar cardiac tissue sample, the spherical shell method is no longer valid. In this last case it has been found that the magnetic eld reaches a peak value B = 1 nT . The heart AP is about t = 300500 ms long, however for some cardiac cells, such as ventricular or rapid response cells, the AP rapid rise occurs in 1 ms, as in the neuronal case. Considering a human heart, a NV sensor positioned on the heart surface should be sensitive to magnetic eld B = 14 nT, with a temporal resolution of about 0.1 ms (in order to be able to trace the time variation t even in the case of the AP rapid rise), and a spatial resolution of about 10 m (which would allow a good reconstruction of the PA propagation, being the heart radius of about 40 mm ). This corresponds to a minimum sensitivity: p p 1=2 = B t ' 14 nT 0:1 ms ' 140 pT Hz (8) min This value can be considered a useful intermediate step for the application of the actual NV-based biosensing tech- nologies, with the aim of reaching sensitivity that allows the detection on neuronal signals. 12 a) b) FIG. 8: a) The model of a spherical heart, taken from the reference . Part of the spherical shell has been cut out to show the heart wall. The black curves indicate the ber orientation. The pink tissue has a transmembrane potential of 20 mV, and the blue tissue has a transmembrane potential of 80 mV. The green curve shows the magnetic eld. The endocardial (inner) surface has radius r , and the epicardial (outer) surface has radius r . b) The magnetic eld over a cross section of the heart. 1 2 The dashed curves indicate the heart inner and outer surfaces. An area 40 mm by 40 mm is shown. IV. METHODS AND BIO-APPLICATIONS To exploit NV-centers for biosensing, it is necessary to set up an optical microscope equipped with an ODMR apparatus featuring a non-resonant laser and a microwave antenna positioned near the diamond sample designed for a microwave source operating in the 2-4 GHz range, as described in section 2.2. The next sections are devoted to the presentation of several biosensing experiments exploiting NV centers in diamonds. Some of them are "proof-of-principle" in vitro tests on a cell culture, while others are experiments carried on living organisms (in vivo experiments). The reviewed experiments will be divided in two categories describing respectively applications exploiting bulk diamond and experiments relying on nanodiamonds. A. Bulk diamond applications As discussed in the previous section, the electromagnetic elds produced by the excitable cells (as neurons, neuronal-like chroman cells, heart cells) even in mammals, are typically extremely weak (pT ). For this reason, several "proof-of-principle" experiments addressed measurements of elds produced by cells with peculiar electro- magnetic properties. 106{109 III II Among the most suitable ones there are magnetotactic bacteria (MTB) containing magnetite ( Fe FeO ) III II or ferrite ( Fe FeS ) bacteria magnetic particles (BMP). The nanometer size of the BMPs is small to generate a single magnetic domain, but sucient to create a permanent magnetic moment ~ . This produces a cell magnetic BMP moment ~ = ~ , given by the sum of the BMP individual dipoles, which is exploited by the MTB to MTB BMP 110,111 orient itself with respect to the earth's magnetic eld . Among the various uses in the biomedical eld, Sage et al. used Magnetospirillum magneticum AMB-1 for bio- magnetic imaging. The bacteria used in this work create magnetic nanoparticles with cubo-octahedral morphology and an average diameter of 50 nm. The experiment was performed both with bacteria dried on the surface of diamond chip implanted with NV centers, as well as with bacteria stored in phosphate-bu ered saline (PBS) and laid on the chip surface (in vitro experiment). The diamond sensor used to perform this experiment was a high-purity single-crystal diamond chip, with a 10 nm 17 3 layer thickness of NV centers. The estimated surface density of nitrogen-vacancy centers was 3  10 cm in the 18 3 case of experiment on bacteria in the liquid medium and 10 cm for dry bacteria. In the case of dry bacteria the objective was to demonstrate the possibility of measuring their static magnetic eld, exploiting ODMR measurements at di erent bias magnetic eld orientations (B = 3.7 mT). bias In the case of live bacteria in the liquid medium, it was shown that it is possible to evaluate the magnetic eld generated by the bacteria dipole ~ along the [111] crystallographic axis of the diamond, when also the bias MTB magnetic eld is oriented along it. Furthermore, cell viability was assessed immediately after magnetic imaging (lasting 4 minutes), using a standard uorescence-based "live-dead" assay obtaining a viability of about 44%. Cells 13 mortality was attributed to the laser heating, since preliminary tests showed that 1 hour exposure to microwaves did not cause substantial cells mortality. Cells vitality was however partially preserved thanks to the strategy used to decouple laser light from the biological sample. Indeed in this set-up the laser impinged on diamond at an angle greater than the critical angle for the diamond{water interface, resulting in its total internal re ection within the diamond. A wide eld optical microscope was used for both MTB samples, with a eld of view of 100  30 m of the sample surface and a resolution of 400 nm. A CMOS camera was used to image the single magnetic nanoparticles inside the MTB. Their magnetic eld was of the order of mT. Thanks to these measurements, the total magnetic moment ~ was determined by numerically tting the modeled eld distribution to the measured ones, with a mean value MTB 17 2 of 5 10 m A. The magnetic eld estimated from the ODMR measurements was compared with a scanning electron microscope (SEM) measurements. The position of nanoparticles revealed by the SEM was used to model the magnetic eld 108,109 they generated . The two measurements were in excellent agreement and their values were compatible with 112,113 the data reported in . This highlights the potential of NV centers, able to perform sub-cellular magnetic eld measurements at room temperature, allowing real-time imaging of magnetic dipole creation, single MTBs chain 111 114,115 dynamics and magnetic particles formation in various organisms . Another "proof-of-principle" test was carried out by Davis et al. with the aim of measuring the magnetic eld generated by iron oxide nanoparticles (IONs) incorporated in murine RAW 264,7 macrophages (a line estab- lished from a tumor induced by Abelson murine leukemia virus and often studied in relation to immune responses). The cells, after having phagocytized the iron ions (about 200 nm in size), were dried on the surface of a bulk diamond. High resolution magnetic imaging was performed exploiting ODMR measurements with a succession of bias magnetic eld (B = 10 mT) for each orientation of NV complex, for a total duration of 2 hours. Projection bias eld maps were combined to form 3 orthogonal eld maps, from which the cellular magnetic moment was obtained. For this experiment, the central dip of the hyper ne transition was used, which allowed to achieve a sensitivity of 17 nT at 1 m in plane resolution, sucient to reveal these magnetic nanoparticles. To extend this technique to diagnostic imaging, Davis et al. performed NV magnetometry on liver specimens from a mouse model of hepatic iron overload, generated through intravenous administration of 900 nm IONs to C57bl/6 mice. To reduce the deposition of optical and thermal energy, the sample was illuminated for only 5 minutes with a duty cycle of 50% and the ODMR technique was carried out with the bias magnetic eld along only one of four NV axes. Furthermore, the laser beam was directed on the sample in total re ection mode. With these choices, time-lapse images of magnetic elds that coalesced within the macrophages after ION internalization were evaluated along that NV axes. This experiment highlights the possibility of study the spatial distribution of iron deposits in the liver and other tissues. This has been a topic of interest in clinical literature as an indicator of disease state and the magnetic resonance imaging is becoming increasingly important in non-invasive quanti cation of tissue iron, overcoming the drawbacks of traditional techniques (liver biopsy) . Barry et al. studied individual neurons of marine worms (Myxicola infundibulum ) and squids (Loligo pealeii ). The marine worm has a long axon , which stretches over its entire length (tens of mm and diameter of about 5 mm). The giant squid neuron (about 0.5 m long) did not extend over the entire length and is isolated following speci c protocols . An initial proof-of-principle test was performed on isolated neurons for both species. The AP is stimulated by means of a current pulse, received by an electrode directly in contact with the neuron. The pulse was generated by a current of about 10 mA, had a duration of about 1 ms and was repeated with a frequency of 0.4 Hz for the worm and 100 Hz for the squid. The AP generation and its propagation was veri ed by micro-electrodes (see Figure 9A). From this axonal AP 119{121 intracellular time trace, the shape of the associated magnetic eld could be modeled (see Figure 9B). This was compared with the experimentally measured magnetic eld, performed by means of the NV-based sensor in contact with the excised single neuron. Traces are shown in Figure 9C and 9D respectively for the worm and the squid neuron. These measurements were performed using the ODMR technique at bias magnetic eld B = 0.7 bias mT, oriented along two diamond axes and perpendicular to the axon axis (being the magnetic eld generated by the AP pulse perpendicular to this last one). In Ref. , Barry et al. carried on also a measurement on a living worm. The worm was directly xed on the diamond and the distance between the neuron and the active NV layer was about 1.2 mm (see Figure 10A). The magnetic eld generated by the propagation of the AP pulse measured by ODMR technique is shown in Figure 10B. It is smaller than the one measured in the excised neuron, but its value is compatible with the increasing sensor distance. The diamond sensor, exploited an electronic grade (N < 5 ppb) single crystal chip, with a NV center layer of 13 17 3 m. This layer had a NV centers density of d = 3 10 cm and a characteristic dephasing time T = 450 ns. The 6 3 11 sensing volume is V = 5 10 cm , consequently the number of potentially stimulated centers was n = 15 10 . 14 FIG. 9: Measured AP voltage and magnetic eld from excised single neurons, taken from the reference . meas A) Measured time trace of intracellular axonal AP voltage  (t) for giant axon from M. infundibulum (worm). in B) Calculated time trace of AP magnetic eld B(t) for M. infundibulum extracted from data in A. C) Measured time trace of AP magnetic eld B(t) for M. infundibulum giant axon with N = 600. avg D) Measured time trace of AP magnetic eld B(t) for L. pealeii (squid) giant axon with N = 375. Gray box indicates avg magnetic artifact from stimulation current. 1=2 Referring to the Equation (7), the fundamental sensitivity limit:  ' 10 fT Hz , while the sensitivity reached 1=2 experimentally was  ' 15 pT Hz , allowing, anyway, a reliable measure of the magnetic elds generated by these animal species (of the order of nT). Further development should be needed for revealing those of human neurons (of the order of pT). FIG. 10: Single-neuron AP magnetic sensing exterior to live intact organism, taken from the reference . A) Overhead view of intact living specimen of M. infundibulum (worm) on top of NV diamond sensor. In con guration shown, animal is stimulated from posterior end by suction electrode, APs propagate toward worm's anterior end, and bipolar electrodes con rm AP stimulation and propagation. (Scale bar 20 mm). B) Measured time trace of AP magnetic eld B(t) from live intact specimen of M. infundibulum for N = 1,650 events. avg B. Nanodiamonds The techniques for the creation of NV centers in diamond are well established also for nanodiamonds (NDs). Nanodiamonds-based sensors exploit colloidal suspensions of single diamond particles of minimum diameter of 4-5 nm, but on average the nanodiamonds typically used in experiment have a size of 50-100 nm. The nanometer size makes nanodiamonds-based sensor of extreme interest for bio-sensing application, as they are potentially usable in vivo experiments. Nonetheless, they have also important drawbacks such as e.g. the increased sensitivity of NV spins to environmental noise. Indeed, while in a bulk diamond the coherence time T is mainly in uenced by the electronic impurities and nuclear spins in the surrounding, for nanodiamonds the coherence time is further reduced due to the surface spin noise. This should be taken into account in the estimation of the sensitivity limit (see Equation (7)). Despite this limitation, nanodiamonds have attracted interest also as a non-toxic alternative to quantum dots for biomedical imaging, as magnetic sensors and nally as drug transporters (thanks to the discovery of the possibility to functionalize the diamonds surface in various ways, exploiting the covalents carbon bonds). The great interest and the exceptional range of applications of NDs is boosting the development of novel fabrication techniques, even if the actual technologies are already able to provide very pure nanodiamonds with controlled surface chemistry at a relative 122,123 low cost . 15 1. Biocompatibility and functionalization studies To understand the perspective in bio-medical application, deep investigation of NDs biocompatibility is required. More speci cally, it is important to understand their di usion in tissues and their acute and long-term biological 124,125 e ect. While bulk diamonds are non-toxic and inert, NDs interaction with cells should be carefully investigated . There is a huge variety of nanodiamond specimens, which di er in dimensions, functionalized surface and potential interaction with the biological sample. 126{128 Various experiments were therefore conducted to evaluate the cell viability, e.g. in HeLa cells (a cell line 43,47 129 deriving from tumoral human cells), in human neurons , in human trachea , in the translucent Caenorhabditis 130 131 elegans worm and intravenous infusions . Brie y, the nanodiamonds of size between 50 and 100 nm have been found to be incorporated by the cells, without producing a signi cant damage. In particular in Guarina et al. an ODMR detection scheme with NV centers in nanodiamonds internalized in hippocampal neurons was performed in suitable conditions (3 mW of excitation power, -20 dBm of continuous- wave MW power). This experiment demonstrated that neuron functionality was not signi cantly a ected by the implementation of the measurement protocol: their spontaneous ring (bursts synchronization) was preserved, as well as the amplitude of spontaneous inhibitory and excitatory events. Even thought some alteration both at the single-cell level and in neuronal networks was observed, this was principally attributed to the e ects of nanoparticles aggregation. The aim of the work was to assess the feasibility of in vitro imaging and targetable drug delivery via nanodiamonds, but the same argument holds for the other sensing applications. Furthermore, if properly functionalized, the NDs can anchor themselves to the surface of the cell sample in the targeted areas. 2. Nanodiamonds applications Once the biocompatibility of nanodiamonds is assessed, it is necessary to understand to which extent the sensing techniques developed for sensor based on NV in bulk diamond can be extended to nanodiamonds based sensors, functionalized and incorporated in the cells of interest. A proof-of-principle demonstration of quantum control techniques to map the intracellular temperature of a neuronal network was performed by Simpson et al. . The NDs were dispersed in cell media in concentration of 6 g/ml, sonicated for few minutes, and then applied to the primary cultures during a routine change of cell media. The 170 nm diameter NDs contained approximately 500 NV centers each. Using ODMR techniques in combination with standard wide- eld microscopy with a eld of view of 80 80 m was possible to observe NV resonance frequency in only 6 seconds. Speci cally, in Ref. the ODMR signal presented two uorescence dips (see section 2.4) because of the strain. This e ect is negligible in bulk diamonds while nanodiamonds crystal lattice su ers strong deformation inducing line splitting. In that paper the two dips, spaced by few MHz, were modeled as a single one with higher spectral broadening. By interpolating the ODMR graph with a Lorentzian function, it was estimated the mean crystal eld splitting D = (2868.590.17) MHz. gs To demonstrate the NV thermo-sensor performance in biological measurement, the temperature of the neuronal solution was reduced by 1:9 C. Repeating the ODMR analysis for a total acquisition time of 12 s, a resonance frequency shift was observed. The respective temperature variation was estimated using the temperature coupling coecient dD=dT ' -74 kHz K (see section 2.5). The distribution reported a mean temperature change of (1:36 0:08) C, consistent with the reduction in environmental temperature. We underline that NDs allowing to create spatial maps of the temperature inside the cells will bring new insight on the understanding of cell activity. There are many biological processes whose knowledge would be enriched by nanoscale thermometry, such as temperature increases following the opening of ion channels , or the correlation of temperature changes and the occurrence of neurological disorders and pathological conditions. Another biological application was reported by Ermakova et al. , using nanodiamonds with NV centers as 132,133 thermo-sensors, exploits optically-induced thermal gradients for thermogenetic neural modulation . This thermal gradient is generated at the transient receptor potential channels (TRP channels): a group of ion channels that are commonly present on the plasma membrane of numerous types of animal cells . A particular specialized form of these ion channels appears to be highly sensitive to temperature changes . Some species of snakes can use TRP channels to detect the thermal build-up caused by infrared IR radiation emitted by nearby prey, allowing them to estimate the direction and distance of the (IR) source . To experimentally recreate this local temperature change and therefore study the TRPs response, Ermakova et al. used IR short pulsed laser. This method, with respect to conventional techniques as environmental heating or TRPs chemical agonists , allows cellular spatial resolution and ultrahigh temporal resolution. The precise temperature control was performed by varying the laser intensity, whose actual thermal impact was monitored by 16 the nitrogen-vacancy complex. This quantum probe (whose dimension was about 300 nm) was integrated on the tip of an optical ber, together with a microwave antenna. The optical ber was positioned near the cell irradiated by the IR laser, allowing a measurement of its temperature by the ODMR technique. In this experiment it was initially evaluated the thermal stimulation via IR laser of two TRP channels of the snake. The TRP channels considered were the Crotalus atrox TRPA1 (caTRPA1) and the Elaphe obsoleta lindheimeri TRPA1 (eolTRPA1). Fluorescent proteins (caTRPA1-IRES-EGFP) had been added to the channels, allowing to monitor the opening and closing of the calcium channels. Thanks to the NV sensors and by slowly changing the cell temperature with properly tuning pulsed laser intensity, it was possible to obtain the threshold temperature, inducing opening of the calcium channels. The threshold temperatures were found to be T =(27.8  0.6) C for caTRPA1 (see Figure 11b) and T =(38.5  0.7) C for eolTRPA1 (see Figure 11c). Once estimated the threshold temperature T , Ermakova et al. proved the technique on other biological samples: mouse neurons and zebra sh larvae, whose thermogenetic activation is induced by TRPA1 channels causing responses. In the case of caTRPA1 channels, the cultured neurons were maintained at a temperature of 27 C lower than the threshold temperature obtained before for this channel; in the case of eolTRPA1-expressing, neurons were kept at basal temperatures of 35.5 C. As expected, they found that the thermal increase induced by the IR laser activates the TRPs channels triggering the generation of the neuronal AP, measured through conventional electrophysiological techniques. When a measurements on live samples is considered, the sample can no longer be kept at the desired temperature, therefore it is necessary to choose the TRP channel suitable for body temperature of the animal species analyzed. As for the zebra sh neurons, whose body temperature is found to be 26 C, the eolTRPA1 channels may be suitable. As for the mammalian brain, the perfect TRP candidate has still to be found. For example, the mouse body temperature is too close to the threshold temperature of eolTRPA1 and it may be desensitized. The results of the application of this technique in living zebra sh showed that it is possible to thermogenetically activate neurons using the IR laser. In particular, the technique demonstrated a spatial resolution of 60 m ( ber size in which the IR laser was focused on the sample), allowing one or few neurons to be stimulated. As for the IR laser intensity, Emarkova et al. observed that 30 mW laser power induced the escape behavior exhibition of 93% of the larvae. NV-based temperature sensors allowed careful monitoring of the temperature reached by the cells with high spatial resolution and temperature sensitivity up to 0.1 C. To preserve cellular integrity and to avoid cell ablation is essential to heat-up the tissues by a few degrees only and for a time interval not exceeding a few minutes. FIG. 11: Activation of snake TRPA1 in cells expressing TRPA1-IRES-EGFP using femtosecond IR laser pulses, taken from the reference . 2+ a) R-GECO1.1 uorescence (black line) re ects Ca dynamics in the cytoplasm with the 20 mW laser beam turned on at t=30 s and o at t=60 s. b,c) With the temperature of HEK293 cells expressing snake TRPA1 increased in a stepwise fashion using properly adjusted IR laser radiation, the activation thresholds of caTRPA1 b) and eolTRPA1 c) were determined. 2+ d) A similar heating of control cells does not induce Ca elevation. The black line is the uorescence response. The red line is the temperature in the medium. Finally, another demonstration of the e ectiveness of NV nanosensors for thermometry comes from Fujiwara et al.'s experiment . A rst test allowed to measure the temperature dynamics inside live C. elegans adults worms 1=2 during environmental temperature changes. The sensitivity reached was 1.4 C Hz . Having obtained this result, Fujiwara et al. successfully determined the temperature increase caused by the worm's thermogenesis under the treatment of mitochondrial uncoupler stimuli. 17 V. TECHNIQUES FOR IMPROVING ODMR SENSITIVITY In this section we discuss some technological solutions to improve the sensitivity of NV-based sensors as well as the precautions to be taken when they are used as a bio-sensor. Equation (7) provides the ultimate sensitivity limit reachable highlighting that the number n of NV centers and their decoherence time T play a key role. To increase n, while maintaining the same spatial resolution, it is necessary to have diamonds with an increased NV centers density. This can be achieved by enhancing the number of nitrogen implanted in the diamond and improving the N-to-NV conversion eciency, minimizing the concentration of residual paramagnetic substitutional nitrogen . In parallel, to increase T , it is also recommended the production of ultra-pure diamonds, with reduced unwanted electronic impurities (e.g. the P1 centers) and nuclear spins impurities 13 139{141 (e.g. the paramagnetic C isotopes, whose natural abundance is about 1.1%) . It is important to note that the NV density increase will necessary worsen the decoherence time of the NVs themselves, because of their mutual interaction. Consequently, an optimal trade-o s between these parameters must be sought. In addition to the NV-density and diamond sample engineering, the sensitivity can be improved by imple- menting speci c experimental techniques, that are based on laser and microwave pulses of particular duration, 28,46,142 synchronized appropriately . For example, if an unknown electromagnetic eld, responsible for the ODMR resonance frequency shift, is constant or 143 38 slowly varying, it is possible to adopt the experimental pulsed ODMR protocols or the Ramsey method instead of the continuous wave (CW ) ODMR . The CW ODMR is the simplest and most widely employed magnetometry method with NV-based sensors, wherein the microwave driving and the optical polarization and readout (laser pumping) occur simultaneously. Although this technique is easy to be implemented, the relative ODMR spectrum dips are a ected by the broadening induced by the continuous exposure of the laser beam and microwave eld on the sample. With pulsed ODMR techniques this broadening e ect is substantially suppressed, allowing to obtain a narrower ODMR spectrum dips and therefore to improve the measurement sensitivity. This protocol uses temporally separated optical laser initializations,  microwave control pulses, and laser readout pulses. The  pulses, whose name derives from the representation of the process on the Bloch sphere, is an oscillating microwave eld that brings the electronic state from the state jm = 0 > to jm = 1 >. Ramsey ODMR spectroscopy, on the other s s hand, consists on in the application of two =2 pulses, separated by a time  . Also the =2 pulse is an oscillating microwave eld that brings the electronic state from the state jm = 0 > to a balanced superposition of jm = +1 > s s and jm = 1 >. By varying the time  , the so-called "Ramsey fringes" are obtained, from which it is possible to extrapolate an estimation of the magnetic elds amplitude. Also this technique allows sensitivity improvement with respect to the CW : the decoupling of the MW from the laser power allows increasing the MW power improving the contrast, without degrading the FWHM. In the case of time-varying electromagnetic elds, there are other even more complex microwave pulse sequences, capable of decoupling the measurement from surrounding spin environment . In this way the decoherence time of the NV centers increases and consequently it becomes possible to interrogate the quantum system for longer times, improving the measurement statistic and therefore the sensitivity. One of these experimental protocols is the 127,146 Hahn Echo sequence , which refocuses the dephasing NVs spin, applying an additional  pulse in the middle of Ramsey sequence. The characteristic time of the spin coherence decay, measured with this protocol, is called T and it is typically one or two orders of magnitude longer than T . Even more complex dynamic decoupling sequences, 147{150 which apply multiple refocusing  pulses further improving the decoherence time T have been devised . Among 147{150 147,151 these, the most famous are the Carr-Purcell-Meiboom-Gill (CPMG) and the XY8 sequences , which di er in the rotation axes (around which the spin rotates): the rst method applies the pulses along the same axis, while the second chooses a di erent one for each  pulses. It is useful to underline that, although these techniques allow to extended the coherence time of the NV centers, they cannot go beyond the spin-lattice relaxation time T , that for an NV spin ensemble in bulk diamond is about 3 ms . The Figure 12 brie y summarizes the above mentioned pulse sequences. It is useful to underline that the sensitivity formula in Equation (7) describes an idealized measurement with a perfect readout mechanism. On the contrary, typically the readout mechanism adds noise in the measurement, that can be described introducing, in the previous equation, the spin-readout delity factor F : 1 1 = (9) nT Keeping the usual optical-readout, but improving the photon collection is expected to increase F (see ref. for di erent methods to improve photon collection). Ancilla-assisted repetitive readout, which is based on mapping the NV spin state to the nuclear spin state, also improves F . Finally, in a more far perspective, quantum methods of 154{157 noise reduction can be applied . 18 CW PULSED RAMSEY HANH ECHO DYNAMICAL DECOULPLING LASER MW ( ) π π/2 π/2 π/2 π π/2 π/2 π π/2 READOUT FIG. 12: Scheme of timing and duration of laser pulses, microwave pulses and reading sequences associated with the most common measurement protocols of the external elds for the NV complex. When the ultimate goal is bio-sensing, some constraints rise limiting the implementation of the above de- scribed pulse sequences. One constraint is the frequency bandwidth. In fact, the dynamic decoupling techniques mentioned above are capable of measuring time-varying external eld only if this time variation is of the order of the time interval separating the  pulses. Furthermore, in order to control the system quantum state, the time between these pulses cannot exceed the coherence of the NV center. Consequently, the frequency of the signal to be measured must be of the order of the coherence time of the NV centers. In the biological case, the electromagnetic elds pulse lasts about 1 ms. This value is very far from T , marking a boundary for the use of these techniques in biological applications. Another constraint is associated to the optical laser power. The higher the laser power the better the sensitivity in measurements with ensembles, since it increases the percentage of excited centers and consequently the uorescence signal. However precautions must be taken to avoid cells and proteins damaging. An ecient solution can be to direct the laser beam towards the diamond sample at an angle allowing total re ection (Brewster angle). In this way only the uorescence emitted by the NV centers travels through the cells, placed on the other diamond 44,77 surface . In this case, however, precise control over sensing volume would be lost, deteriorating spatial resolution. In a standard con guration, where the laser impinges perpendicularly on the sample, it is necessary to limit the optical power reaching the cells to few mW. In this regard, Figure 13 shows a sensitivity curve versus the laser optical power, obtained by adopting the technique described in Moreva et al. . As anticipated in the introduction, FIG. 13: Temperature sensitivity versus the laser excitation power at 532 nm. The inset shows the inverse of the thermal sensitivity versus the excitation laser power. 19 the application of a transverse bias magnetic eld B ' 3 mT , allows to improve the sensitivity of a the NV center bias based thermo-sensor with respect to other standard techniques in CW regime. In Ref. , the temperature sensitivity 1=2 3 reached is  ' 4:8 mK=Hz in a sensing volume of 1m , obtained at a power level (80 mW) that can present biocompatibility problems. However, the sensitivity obtained is even beyond the one required to monitor biological mechanisms, usually requiring sensitivities of the order of 1 C. Figure 13 shows that it is possible to perform the temperature measurement with a lower laser power, nding an ideal compromise between the temperature sensitivity and laser intensity impinging on the cell sample. Indeed, with a power of a few mW it is already possible to discriminate biological processes with a sensitivity of the tenth of a degree. VI. CONCLUSION Sensors based on the NV centers in arti cial diamonds are one of the emerging quantum technologies of huge potential interest in biological applications, thanks to both their practicality and their technical performances. In fact, the capability to initialize and read out optically the spin state at room temperature, makes the use of this quantum sensor convenient and powerful even for biological applications. Furthermore, the levels of sensitivity and spatial resolution achieved are extremely high, which in principle allows potential application towards the detection of very weak electromagnetic elds as the one generated by mammalian, and potentially human, cells. Even if an eventual use of NV sensors for the detection of biological electric elds is more problematic due to its weak coupling constant, regarding the magnetic eld sensing and especially temperature measurements astonishing results have already been achieved. Indeed, the thermal gradients generated by biological phenomena can be reliably observed thanks to the actual sensitivity of the NV-based sensors. This is of the utmost importance because localized intracellular temperature gradients may a ect neuronal functionality (including vesicular dynamics and neurotransmitter release) 158,159 or may provide indirect measurement of mitochondrial activity . Regarding the detection of bio-magnetic elds, the NV-based sensors have already shown good results with peculiar biological cells, presenting either an intrinsic magnetic eld (magnetitotactic bacteria) or a generated magnetic eld in axon of squids or long worms, much larger than the one generated in the human ones. The improvement of these devices suggests the possibility of exploiting NV-based sensors also for the detection of weaker but more fascinating biological magnetic elds. In particular, an estimate of the cardiac magnetic eld that is generated on the heart surface was here reported. This value is in the range of present measurement capability exploiting the NV center properties, exploiting optimized diamond sample engineering and the adoption of pulsed measurement protocols in order to improve the diamond coherence time. Furthermore, we have analyzed the magnetic eld associated to human neuron activity. The weakness of these elds requires further improvements of the measurement technique in the case of the single AP, while measurement of clustered channels is likely a reasonable target for the actual technology. However, the considerable interest in the neuronal eld detection as diagnostic and therapeutic tools for neurodegenerative diseases and aging e ects, together with the recent years progress of these techniques (partially covered by this review), is expected to boost the technological developments and eventually the market success of quantum assisted biosensing based on NVs. Acknowledgements This work has received funding from the European Union's PATHOS EU H2020 FET-OPEN grant no. 828946 and Horizon 2020, from the EMPIR Participating States in the context of the projects EMPIR-17FUN06 "SIQUST" and from the project Piemonte Quantum Enabling Technologies (PiQuET) funded by the Piemonte Region. Richard T Mathias, Ira S Cohen, and Carlos Oliva. Limitations of the whole cell patch clamp technique in the control of intracellular concentrations. Biophysical journal, 58(3):759, 1990. Michael Fejtl, Alfred Stett, Wilfried Nisch, Karl-Heinz Boven, and Andreas M oller. On micro-electrode array revival: its development, sophistication of recording, and stimulation. In Advances in network electrophysiology, pages 24{37. Springer, Micha E Spira and Aviad Hai. Multi-electrode array technologies for neuroscience and cardiology. Nature nanotechnology, 8(2):83, 2013. M Vengalattore, JM Higbie, SR Leslie, J Guzman, LE Sadler, and DM Stamper-Kurn. High-resolution magnetometry with a spinor bose-einstein condensate. Physical review letters, 98(20):200801, 2007. MI Faley, U Poppe, K Urban, DN Paulson, and RL Fagaly. A new generation of the hts multilayer dc-squid magnetometers and gradiometers. In Journal of Physics: Conference Series, volume 43, page 1199. IOP Publishing, 2006. F Baudenbacher, LE Fong, JR Holzer, and M Radparvar. Monolithic low-transition-temperature superconducting magne- tometers for high resolution imaging magnetic elds of room temperature samples. Applied Physics Letters, 82(20):3487{ 3489, 2003. 20 Svenja Knappe, Tilmann H Sander, Olaf Kosch, Frank Wiekhorst, John Kitching, and Lutz Trahms. Cross-validation of microfabricated atomic magnetometers with superconducting quantum interference devices for biomagnetic applications. Applied Physics Letters, 97(13):133703, 2010. Anri Nakajima. Application of single-electron transistor to biomolecule and ion sensors. Applied Sciences, 6(4):94, 2016. Laura Martinez Maestro, Emma Martin Rodriguez, Francisco Sanz Rodriguez, MC Iglesias-de la Cruz, Angeles Juarranz, Ra k Naccache, Fiorenzo Vetrone, Daniel Jaque, John A Capobianco, and Jose Garcia Sole. Cdse quantum dots for two-photon uorescence thermal imaging. Nano letters, 10(12):5109{5115, 2010. CDS Brites, A Mill an, and LD Carlos. Lanthanides in luminescent thermometry. In Handbook on the Physics and Chemistry of Rare Earths, volume 49, pages 339{427. Elsevier, 2016. Michael Hines. Neuron|a program for simulation of nerve equations. In Neural systems: Analysis and modeling, pages 127{136. Springer, 1993. Michael L Hines and Nicholas T Carnevale. The neuron simulation environment. Neural computation, 9(6):1179{1209, F Santamaria and JM Bower. Hodgkin-huxley models. In Encyclopedia of Neuroscience, pages 1173{1180. Elsevier Ltd, Blaise Aguera  y Arcas, Adrienne L Fairhall, and William Bialek. Computation in a single neuron: Hodgkin and huxley revisited. Neural Computation, 15(8):1715{1749, 2003. Single-Channel Recording. edited by b. sakmann and n. neher, 1995. Stephen G Waxman. Demyelinating diseases|new pathological insights, new therapeutic targets, 1998. Eugene A Kiyatkin. Brain temperature: from physiology and pharmacology to neuropathology. In Handbook of clinical neurology, volume 157, pages 483{504. Elsevier, 2018. PC Maurer, M Kubo, MD Lukin, HJ Noh, G Kucsko, NY Yao, PK Lo, and H Park. Nanometer scale quantum thermometry in a living cell. 2013. Tingting Bai and Ning Gu. Micro/nanoscale thermometry for cellular thermal sensing. Small, 12(34):4590{4610, 2016. Yuzhou Wu, Fedor Jelezko, Martin B Plenio, and Tanja Weil. Diamond quantum devices in biology. Angewandte Chemie International Edition, 55(23):6586{6598, 2016. Alexios Beveratos, Rosa Brouri, Thierry Gacoin, Jean-Philippe Poizat, and Philippe Grangier. Nonclassical radiation from diamond nanocrystals. Physical Review A, 64(6):061802, 2001. Alexios Beveratos, Sergei Kuhn,  Rosa Brouri, Thierry Gacoin, J-P Poizat, and Philippe Grangier. Room temperature stable single-photon source. The European Physical Journal D-Atomic, Molecular, Optical and Plasma Physics, 18(2):191{196, Alexander M Zaitsev. Optical properties of diamond: a data handbook. Springer Science & Business Media, 2013. Ettore Bernardi, Ekaterina Moreva, Paolo Traina, Giulia Petrini, Sviatoslav Ditalia Tchernij, Jacopo Forneris, Zelijko Pastuovic, Ivo Pietro Degiovanni, Paolo Olivero, and M Genovese. Biocompatible technique for nanoscale magnetic eld sensing with nitrogen-vacancy centers. arXiv preprint arXiv:2005.13230, 2020. Shu-Jung Yu, Ming-Wei Kang, Huan-Cheng Chang, Kuan-Ming Chen, and Yueh-Chung Yu. Bright uorescent nanodia- monds: no photobleaching and low cytotoxicity. Journal of the American Chemical Society, 127(50):17604{17605, 2005. Frank C Barone, Cezary Marcinkiewicz, Jie Li, Yi Feng, Mark Sternberg, Peter I Lelkes, David Rosenbaum-Halevi, Jonathan A Gerstenhaber, and Giora Z Feuerstein. Long-term biocompatibility of uorescent diamonds-(nv)-z~ 800 nm in rats: survival, morbidity, histopathology, particle distribution and excretion studies (part iv). International journal of nanomedicine, 14:1163, 2019. Christian G Specht, Oliver A Williams, Richard B Jackman, and Ralf Schoepfer. Ordered growth of neurons on diamond. Biomaterials, 25(18):4073{4078, 2004. Romana Schirhagl, Kevin Chang, Michael Loretz, and Christian L Degen. Nitrogen-vacancy centers in diamond: nanoscale sensors for physics and biology. Annual review of physical chemistry, 65:83{105, 2014. Georg Kucsko, Peter C Maurer, Norman Ying Yao, MICHAEL Kubo, Hyun Jong Noh, Po Kam Lo, Hongkun Park, and Mikhail D Lukin. Nanometre-scale thermometry in a living cell. Nature, 500(7460):54{58, 2013. Jan Jeske, Jared H Cole, and Andrew D Greentree. Laser threshold magnetometry. New Journal of Physics, 18(1):013015, Marcus W Doherty, Neil B Manson, Paul Delaney, Fedor Jelezko, J org Wrachtrup, and Lloyd CL Hollenberg. The nitrogen- vacancy colour centre in diamond. Physics Reports, 528(1):1{45, 2013. Thomas Wolf, Philipp Neumann, and J org Wrachtrup. Diamond magnetometer, December 26 2017. US Patent 9,851,418. A Gruber, A Dr abenstedt, C Tietz, L Fleury, J Wrachtrup, and C Von Borczyskowski. Scanning confocal optical microscopy and magnetic resonance on single defect centers. Science, 276(5321):2012{2014, 1997. MW Doherty, F Dolde, H Fedder, Fedor Jelezko, J Wrachtrup, NB Manson, and LCL Hollenberg. Theory of the ground-state spin of the nv- center in diamond. Physical Review B, 85(20):205203, 2012. Gopalakrishnan Balasubramanian, IY Chan, Roman Kolesov, Mohannad Al-Hmoud, Julia Tisler, Chang Shin, Changdong Kim, Aleksander Wojcik, Philip R Hemmer, Anke Krueger, et al. Nanoscale imaging magnetometry with diamond spins under ambient conditions. Nature, 455(7213):648{651, 2008. Linh My Pham, David Le Sage, Paul L Stanwix, Tsun Kwan Yeung, D Glenn, Alexei Trifonov, Paola Cappellaro, Philip R Hemmer, Mikhail D Lukin, Hongkun Park, et al. Magnetic eld imaging with nitrogen-vacancy ensembles. New Journal of Physics, 13(4):045021, 2011. Jennifer M Schloss, John F Barry, Matthew J Turner, and Ronald L Walsworth. Simultaneous broadband vector magne- tometry using solid-state spins. Physical Review Applied, 10(3):034044, 2018. 21 JM Taylor, P Cappellaro, L Childress, L Jiang, D Budker, PR Hemmer, A Yacoby, R Walsworth, and MD Lukin. High- sensitivity diamond magnetometer with nanoscale resolution. Nature Physics, 4(10):810{816, 2008. CL Degen. Scanning magnetic eld microscope with a diamond single-spin sensor. Applied Physics Letters, 92(24):243111, Jeronimo R Maze, Paul L Stanwix, James S Hodges, Seungpyo Hong, Jacob M Taylor, Paola Cappellaro, Liang Jiang, MV Gurudev Dutt, Emre Togan, AS Zibrov, et al. Nanoscale magnetic sensing with an individual electronic spin in diamond. Nature, 455(7213):644{647, 2008. Florian Dolde, Helmut Fedder, Marcus W Doherty, Tobias N obauer, Florian Rempp, Gopalakrishnan Balasubramanian, Thomas Wolf, Friedemann Reinhard, Lloyd CL Hollenberg, Fedor Jelezko, et al. Electric- eld sensing using single diamond spins. Nature Physics, 7(6):459{463, 2011. Victor M Acosta, Erik Bauch, Micah P Ledbetter, Amir Waxman, L-S Bouchard, and Dmitry Budker. Temperature dependence of the nitrogen-vacancy magnetic resonance in diamond. Physical review letters, 104(7):070801, 2010. Laura Guarina, C Calorio, D Gavello, E Moreva, P Traina, A Battiato, S Ditalia Tchernij, J Forneris, M Gai, F Picollo, et al. Nanodiamonds-induced e ects on neuronal ring of mouse hippocampal microcircuits. Scienti c Reports, 8(1):1{14, David Le Sage, Koji Arai, David R Glenn, Stephen J DeVience, Linh M Pham, Lilah Rahn-Lee, Mikhail D Lukin, Amir Yacoby, Arash Komeili, and Ronald L Walsworth. Optical magnetic imaging of living cells. Nature, 496(7446):486{489, Hunter C Davis, Pradeep Ramesh, Aadyot Bhatnagar, Audrey Lee-Gosselin, John F Barry, David R Glenn, Ronald L Walsworth, and Mikhail G Shapiro. Mapping the microscale origins of magnetic resonance image contrast with subcellular diamond magnetometry. Nature communications, 9(1):1{9, 2018. John F Barry, Jennifer M Schloss, Erik Bauch, Matthew J Turner, Connor A Hart, Linh M Pham, and Ronald L Walsworth. Sensitivity optimization for nv-diamond magnetometry. Reviews of Modern Physics, 92(1):015004, 2020. David A Simpson, Emma Morrisroe, Julia M McCoey, Alain H Lombard, Dulini C Mendis, Francois Treussart, Liam T Hall, Steven Petrou, and Lloyd CL Hollenberg. Non-neurotoxic nanodiamond probes for intraneuronal temperature mapping. ACS nano, 11(12):12077{12086, 2017. Yulia G Ermakova, Aleksandr A Lanin, Ilya V Fedotov, Matvey Roshchin, Ilya V Kelmanson, Dmitry Kulik, Yulia A Bogdanova, Arina G Shokhina, Dmitry S Bilan, Dmitry B Staroverov, et al. Thermogenetic neurostimulation with single- cell resolution. Nature communications, 8(1):1{15, 2017. Masazumi Fujiwara, Simo Sun, Alexander Dohms, Yushi Nishimura, Ken Suto, Yuka Takezawa, Keisuke Oshimi, Li Zhao, Nikola Sadzak, Yumi Umehara, Yoshio Teki, Naoki Komatsu, Oliver Benson, Yutaka Shikano, and Eriko Kage-Nakadai. Real-time nanodiamond thermometry probing in-vivo thermogenic responses, 2020. JoHoNo Loubser and JoA van Wyk. Electron spin resonance in the study of diamond. Reports on Progress in Physics, 41(8):1201, 1978. Erik Bauch, Connor A Hart, Jennifer M Schloss, Matthew J Turner, John F Barry, Pauli Kehayias, Swati Singh, and Ronald L Walsworth. Ultralong dephasing times in solid-state spin ensembles via quantum control. Physical Review X, 8(3):031025, 2018. Chang S Shin, Mark C Butler, Hai-Jing Wang, Claudia E Avalos, Scott J Seltzer, Ren-Bao Liu, Alexander Pines, and Vikram S Bajaj. Optically detected nuclear quadrupolar interaction of n 14 in nitrogen-vacancy centers in diamond. Physical Review B, 89(20):205202, 2014. S Felton, AM Edmonds, ME Newton, PM Martineau, D Fisher, DJ Twitchen, and JM Baker. Hyper ne interaction in the ground state of the negatively charged nitrogen vacancy center in diamond. Physical Review B, 79(7):075203, 2009. Xing-Fei He, Neil B Manson, and Peter TH Fisk. Paramagnetic resonance of photoexcited n-v defects in diamond. ii. hyper ne interaction with the n 14 nucleus. Physical Review B, 47(14):8816, 1993. Victor M Acosta, Erik Bauch, Micah P Ledbetter, Charles Santori, K-MC Fu, Paul E Barclay, Raymond G Beausoleil, H elo se Linget, Jean Francois Roch, Francois Treussart, et al. Diamonds with a high density of nitrogen-vacancy centers for magnetometry applications. Physical Review B, 80(11):115202, 2009. Ekaterina Moreva. The biosensing with nv centers in diamond: Related challenges. International Journal of Quantum Information, 18(01):1941023, 2020. Varun KA Sreenivasan, Andrei V Zvyagin, and Ewa M Goldys. Luminescent nanoparticles and their applications in the life sciences. Journal of Physics: Condensed Matter, 25(19):194101, 2013. Eric Van Oort and Max Glasbeek. Electric- eld-induced modulation of spin echoes of nv centers in diamond. Chemical Physics Letters, 168(6):529{532, 1990. X-D Chen, C-H Dong, F-W Sun, C-L Zou, J-M Cui, Z-F Han, and G-C Guo. Temperature dependent energy level shifts of nitrogen-vacancy centers in diamond. Applied Physics Letters, 99(16):161903, 2011. Yan-Kai Tzeng, Pei-Chang Tsai, Hsiou-Yuan Liu, Oliver Y Chen, Hsiang Hsu, Fu-Goul Yee, Ming-Shien Chang, and Huan-Cheng Chang. Time-resolved luminescence nanothermometry with nitrogen-vacancy centers in nanodiamonds. Nano letters, 15(6):3945{3952, 2015. E. Moreva, E. Bernardi, P. Traina, A. Sosso, S. Ditalia Tchernij, J. Forneris, F. Picollo, G. Brida, Z. Pastuovi c, I. P. Degiovanni, P. Olivero, and M. Genovese. Practical applications of quantum sensing: A simple method to enhance the sensitivity of nitrogen-vacancy-based temperature sensors. Phys. Rev. Applied, 13:054057, May 2020. E Moreva, E Bernardi, P Traina, A Sosso, S Ditalia Tchernij, J Forneris, F Picollo, G Brida, Z Pastuovic, IP Degio- vanni, et al. Practical applications of quantum sensing: a simple method to enhance sensitivity of nitrogen-vacancy-based temperature sensors. arXiv preprint arXiv:1912.10887, 2019. 22 Erwin Neher and Bert Sakmann. The patch clamp technique. Scienti c American, 266(3):44{51, 1992. Alfred Stett, Ulrich Egert, Elke Guenther, Frank Hofmann, Thomas Meyer, Wilfried Nisch, and Hugo Haemmerle. Biological application of microelectrode arrays in drug discovery and basic research. Analytical and bioanalytical chemistry, 377(3):486{ 495, 2003. Jubayer Jalil, Yong Zhu, Chandima Ekanayake, and Yong Ruan. Sensing of single electrons using micro and nano tech- nologies: a review. Nanotechnology, 28(14):142002, 2017. Yoshiaki Takei, Satoshi Arai, Atsushi Murata, Masao Takabayashi, Kotaro Oyama, Shin'ichi Ishiwata, Shinji Takeoka, and Madoka Suzuki. A nanoparticle-based ratiometric and self-calibrated uorescent thermometer for single living cells. ACS nano, 8(1):198{206, 2014. Duan P Chen, Robert S Eisenberg, Joseph W Jerome, and Chi-Wang Shu. Hydrodynamic model of temperature change in open ionic channels. Biophysical Journal, 69(6):2304, 1995. Ahmed El Hady and Benjamin B Machta. Mechanical surface waves accompany action potential propagation. Nature communications, 6:6697, 2015. Ezia Guatteo, Kenny KH Chung, Tharushini K Bowala, Giorgio Bernardi, Nicola B Mercuri, and Janusz Lipski. Tempera- ture sensitivity of dopaminergic neurons of the substantia nigra pars compacta: involvement of transient receptor potential channels. Journal of neurophysiology, 94(5):3069{3080, 2005. Jon S Donner, Sebastian A Thompson, Mark P Kreuzer, Guillaume Ba ou, and Romain Quidant. Mapping intracellular temperature using green uorescent protein. Nano letters, 12(4):2107{2111, 2012. Jui-Ming Yang, Haw Yang, and Liwei Lin. Quantum dot nano thermometers reveal heterogeneous local thermogenesis in living cells. ACS nano, 5(6):5067{5071, 2011. Soo Ho Kim, Jermim Noh, Min Ku Jeon, Ki Woong Kim, Luke P Lee, and Seong Ihl Woo. Micro-raman thermome- try for measuring the temperature distribution inside the microchannel of a polymerase chain reaction chip. Journal of Micromechanics and Microengineering, 16(3):526, 2006. Fiorenzo Vetrone, Ra k Naccache, Alicia Zamarron, Angeles Juarranz de la Fuente, Francisco Sanz-Rodriguez, Laura Martinez Maestro, Emma Martin Rodriguez, Daniel Jaque, Jose Garciia Sole, and John A Capobianco. Temperature sensing using uorescent nanothermometers. ACS nano, 4(6):3254{3258, 2010. Laura Martinez Maestro, Qiming Zhang, Xiangping Li, Daniel Jaque, and Min Gu. Quantum-dot based nanothermometry in optical plasmonic recording media. Applied Physics Letters, 105(18):181110, 2014. Ellen C Jensen. Use of uorescent probes: their e ect on cell biology and limitations. The Anatomical Record: Advances in Integrative Anatomy and Evolutionary Biology, 295(12):2031{2036, 2012. Tytus Bernas, Bartlomiej P Rajwa, Elikplimki K Asem, and Joseph Paul Robinson. Loss of image quality in photobleaching during microscopic imaging of uorescent probes bound to chromatin. Journal of biomedical optics, 10(6):064015, 2005. John F Barry, Matthew J Turner, Jennifer M Schloss, David R Glenn, Yuyu Song, Mikhail D Lukin, Hongkun Park, and Ronald L Walsworth. Optical magnetic detection of single-neuron action potentials using quantum defects in diamond. Proceedings of the National Academy of Sciences, 113(49):14133{14138, 2016. DM Toyli, DJ Christle, A Alkauskas, BB Buckley, CG Van de Walle, and DD Awschalom. Measurement and control of single nitrogen-vacancy center spins above 600 k. Physical Review X, 2(3):031001, 2012. Taras Plakhotnik, Marcus W Doherty, Jared H Cole, Robert Chapman, and Neil B Manson. All-optical thermometry and thermal properties of the optically detected spin resonances of the nv{center in nanodiamond. Nano letters, 14(9):4989{4996, Jan Muller,  Marco Ballini, Paolo Livi, Yihui Chen, Milos Radivojevic, Amir Shadmani, Vijay Viswam, Ian L Jones, Michele Fiscella, Roland Diggelmann, et al. High-resolution cmos mea platform to study neurons at subcellular, cellular, and network levels. Lab on a Chip, 15(13):2767{2780, 2015. Gabriel Bertotti, Dmytro Velychko, Norman Dodel, Stefan Keil, Dirk Wolansky, Bernd Tillak, Matthias Schreiter, Andreas Grall, Peter Jesinger, Sebastian R ohler, et al. A cmos-based sensor array for in-vitro neural tissue interfacing with 4225 recording sites and 1024 stimulation sites. In 2014 IEEE Biomedical Circuits and Systems Conference (BioCAS) Proceedings, pages 304{307. IEEE, 2014. Douglas J Bakkum, Urs Frey, Milos Radivojevic, Thomas L Russell, Jan Muller,  Michele Fiscella, Hirokazu Takahashi, and Andreas Hierlemann. Tracking axonal action potential propagation on a high-density microelectrode array across hundreds of sites. Nature communications, 4(1):1{12, 2013. Gy orgy Buzs aki, Eran Stark, Antal Ber enyi, Dion Khodagholy, Daryl R Kipke, Euisik Yoon, and Kensall D Wise. Tools for probing local circuits: high-density silicon probes combined with optogenetics. Neuron, 86(1):92{105, 2015. Milos Radivojevic, David J ackel, Michael Altermatt, Jan Muller,  Vijay Viswam, Andreas Hierlemann, and Douglas J Bakkum. Electrical identi cation and selective microstimulation of neuronal compartments based on features of extracellular action potentials. Scienti c reports, 6:31332, 2016. Giulia Tomagra, Pietro Apr a, Al o Battiato, Cecilia Colla Ruvolo, Alberto Pasquarelli, Andrea Marcantoni, Emilio Car- bone, Valentina Carabelli, Paolo Olivero, and Federico Picollo. Micro graphite-patterned diamond sensors: Towards the simultaneous in vitro detection of molecular release and action potentials generation from excitable cells. Carbon, 152:424{ 433, 2019. Matthew N Rasband and James S Trimmer. Developmental clustering of ion channels at and near the node of ranvier. Developmental biology, 236(1):5{16, 2001. Daisuke Sato, Gonzalo Hern andez-Hern andez, Collin Matsumoto, Sendoa Tajada, Claudia M Moreno, Rose E Dixon, Samantha O'Dwyer, Manuel F Navedo, James S Trimmer, Colleen E Clancy, et al. A stochastic model of ion channel cluster formation in the plasma membrane. Journal of General Physiology, 151(9):1116{1134, 2019. 23 Bokai Zhang, Xi Feng, Hang Yin, Zhenpeng Ge, Yanhuan Wang, Zhiqin Chu, Helena Raabova, Jan Vavra, Petr Cigler, Renbao Liu, et al. Anchored but not internalized: shape dependent endocytosis of nanodiamond. Scienti c reports, 7:46462, Torsten Rendler, Jitka Neburkova, Ondrej Zemek, Jan Kotek, Andrea Zappe, Zhiqin Chu, Petr Cigler, and J org Wrachtrup. Optical imaging of localized chemical events using programmable diamond quantum nanosensors. Nature communications, 8(1):1{9, 2017. Jasmina Isakovic, Ian Dobbs-Dixon, Dipesh Chaudhury, and Dinko Mitrecic. Modeling of inhomogeneous electromagnetic elds in the nervous system: a novel paradigm in understanding cell interactions, disease etiology and therapy. Scienti c reports, 8(1):1{20, 2018. Mursel  Karadas, Adam M Wojciechowski, Alexander Huck, Nils Ole Dalby, Ulrik Lund Andersen, and Axel Thielscher. Feasibility and resolution limits of opto-magnetic imaging of neural network activity in brain slices using color centers in diamond. Scienti c reports, 8(1):1{14, 2018. LT Hall, GCG Beart, EA Thomas, DA Simpson, LP McGuinness, JH Cole, JH Manton, RE Scholten, Fedor Jelezko, J org Wrachtrup, et al. High spatial and temporal resolution wide- eld imaging of neuron activity using quantum nv-diamond. Scienti c reports, 2:401, 2012. Matti H am al ainen, Riitta Hari, Risto J Ilmoniemi, Jukka Knuutila, and Olli V Lounasmaa. Magnetoencephalogra- phy|theory, instrumentation, and applications to noninvasive studies of the working human brain. Reviews of modern Physics, 65(2):413, 1993. Rolf Simon Schoenfeld and Wolfgang Harneit. Real time magnetic eld sensing and imaging using a single spin in diamond. Physical review letters, 106(3):030802, 2011. Dmitry Budker and Michael Romalis. Optical magnetometry. Nature physics, 3(4):227{234, 2007. Jukka Sarvas. Basic mathematical and electromagnetic concepts of the biomagnetic inverse problem. Physics in Medicine & Biology, 32(1):11, 1987. Christoph M Michel, Micah M Murray, G oran Lantz, Sara Gonzalez, Laurent Spinelli, and Rolando Grave de Peralta. Eeg source imaging. Clinical neurophysiology, 115(10):2195{2222, 2004. Natalia A Trayanova, Bradley J Roth, and Lisa J Malden. The response of a spherical heart to a uniform electric eld: a bidomain analysis of cardiac stimulation. IEEE transactions on biomedical engineering, 40(9):899{908, 1993. Dan Xu and Bradley J Roth. The magnetic eld produced by the heart and its in uence on mri. Mathematical Problems in Engineering, 2017, 2017. Ryan A Murdick and BJ Roth. A comparative model of two mechanisms from which a magnetic eld arises in the heart. Journal of applied physics, 95(9):5116{5122, 2004. Bradley J Roth. A comparison of two boundary conditions used with the bidomain model of cardiac tissue. Annals of biomedical engineering, 19(6):669{678, 1991. Wanda Krassowska and John C Neu. E ective boundary conditions for syncytial tissues. IEEE transactions on biomedical engineering, 41(2):143{150, 1994. Bradley J Roth. Electrical conductivity values used with the bidomain model of cardiac tissue. IEEE Transactions on Biomedical Engineering, 44(4):326{328, 1997. Krista Kay McBride, Bradley J Roth, VY Sidorov, John P Wikswo, and Franz J Baudenbacher. Measurements of trans- membrane potential and magnetic eld at the apex of the heart. Biophysical journal, 99(10):3113{3118, 2010. Jenny R Holzer, Luis E Fong, Veniamin Y Sidorov, John P Wikswo Jr, and Franz Baudenbacher. High resolution magnetic images of planar wave fronts reveal bidomain properties of cardiac tissue. Biophysical journal, 87(6):4326{4332, 2004. Roger B Proksch, TE Sch a er, BM Moskowitz, ED Dahlberg, Dennis A Bazylinski, and Richard B Frankel. Magnetic force microscopy of the submicron magnetic assembly in a magnetotactic bacterium. Applied Physics Letters, 66(19):2582{2584, Lisa Qian, Beena Kalisky, Amanda Hamilton, Bo Dwyer, AC Matin, and Kathryn Moler. Magnetic characterization of individual magnetotactic bacteria. In APS Meeting Abstracts, 2011. Karen P Lam, Adam P Hitchcock, Martin Obst, John R Lawrence, George DW Swerhone, Gary G Leppard, Tolek Tyliszczak, Chithra Karunakaran, Jian Wang, Konstantin Kaznatcheev, et al. Characterizing magnetism of individual magnetosomes by x-ray magnetic circular dichroism in a scanning transmission x-ray microscope. Chemical Geology, 270(1- 4):110{116, 2010. Rafal E Dunin-Borkowski, Martha R McCartney, Richard B Frankel, Dennis A Bazylinski, Mihaly Posfai, and Peter R Buseck. Magnetic microstructure of magnetotactic bacteria by electron holography. Science, 282(5395):1868{1870, 1998. Arash Komeili. Molecular mechanisms of compartmentalization and biomineralization in magnetotactic bacteria. FEMS microbiology reviews, 36(1):232{255, 2012. Damien Faivre and Dirk Schuler. Magnetotactic bacteria and magnetosomes. Chemical Reviews, 108(11):4875{4898, 2008. A Krichevsky, MJ Smith, LJ Whitman, MB Johnson, TW Clinton, LL Perry, BM Applegate, K O'Connor, and LN Csonka. Trapping motile magnetotactic bacteria with a magnetic recording head. Journal of applied physics, 101(1):014701, 2007. Bruce M Moskowitz, Richard B Frankel, and Dennis A Bazylinski. Rock magnetic criteria for the detection of biogenic magnetite. Earth and Planetary Science Letters, 120(3-4):283{300, 1993. Mih aly P osfai and Rafal E Dunin-Borkowski. Magnetic nanocrystals in organisms. Elements, 5(4):235{240, 2009. Stephan HK Eder, Herv e Cadiou, Airina Muhamad, Peter A McNaughton, Joseph L Kirschvink, and Michael Winklhofer. Magnetic characterization of isolated candidate vertebrate magnetoreceptor cells. Proceedings of the National Academy of Sciences, 109(30):12022{12027, 2012. Nilesh R Ghugre and John C Wood. Relaxivity-iron calibration in hepatic iron overload: probing underlying biophysical 24 mechanisms using a monte carlo model. Magnetic resonance in medicine, 65(3):837{847, 2011. JAC Nicol. The giant nerve- bres in the central nervous system of myxicola (polychaeta, sabellidae). Quarterly Journal of Microscopical Science, 3(5):1{45, 1948. Yuyu Song and Scott T Brady. Analysis of microtubules in isolated axoplasm from the squid giant axon. In Methods in cell biology, volume 115, pages 125{137. Elsevier, 2013. KR Swinney and JP Wikswo Jr. A calculation of the magnetic eld of a nerve action potential. Biophysical journal, 32(2):719{731, 1980. Bradley J Roth and John P Wikswo Jr. The magnetic eld of a single axon. a comparison of theory and experiment. Biophysical journal, 48(1):93{109, 1985. John P Wikswo Jr and Bradley J Roth. Magnetic determination of the spatial extent of a single cortical current source: A theoretical analysis. Electroencephalography and clinical neurophysiology, 69(3):266{276, 1988. ML Markham, JM Dodson, GA Scarsbrook, DJ Twitchen, G Balasubramanian, F Jelezko, and J Wrachtrup. Cvd diamond for spintronics. Diamond and Related Materials, 20(2):134{139, 2011. Jagdish Narayan and Anagh Bhaumik. Novel synthesis and properties of pure and nv-doped nanodiamonds and other nanostructures. Materials Research Letters, 5(4):242{250, 2017. Vadym N Mochalin, Olga Shenderova, Dean Ho, and Yury Gogotsi. The properties and applications of nanodiamonds. Nature nanotechnology, 7(1):11, 2012. Amanda M Schrand, Houjin Huang, Cataleya Carlson, John J Schlager, Eiji Osawa, Saber M Hussain, and Liming Dai. Are diamond nanoparticles cytotoxic? The journal of physical chemistry B, 111(1):2{7, 2007. Felix Neugart, Andrea Zappe, Fedor Jelezko, C Tietz, Jean Paul Boudou, Anke Krueger, and J org Wrachtrup. Dynamics of diamond nanoparticles in solution and cells. Nano letters, 7(12):3588{3591, 2007. Liam P McGuinness, Yuling Yan, Alastair Stacey, David A Simpson, Liam T Hall, Dougal Maclaurin, Steven Prawer, P Mulvaney, J Wrachtrup, F Caruso, et al. Quantum measurement and orientation tracking of uorescent nanodiamonds inside living cells. Nature nanotechnology, 6(6):358, 2011. Wesley Wei-Wen Hsiao, Yuen Yung Hui, Pei-Chang Tsai, and Huan-Cheng Chang. Fluorescent nanodiamond: a versatile tool for long-term cell tracking, super-resolution imaging, and nanoscale temperature sensing. Accounts of chemical research, 49(3):400{407, 2016. Yuan Yuan, Xiang Wang, Guang Jia, Jia-Hui Liu, Tiancheng Wang, Yiqun Gu, Sheng-Tao Yang, Sen Zhen, Haifang Wang, and Yuanfang Liu. Pulmonary toxicity and translocation of nanodiamonds in mice. Diamond and Related Materials, 19(4):291{299, 2010. Nitin Mohan, Chao-Sheng Chen, Hsiao-Han Hsieh, Yi-Chun Wu, and Huan-Cheng Chang. In vivo imaging and toxicity assessments of uorescent nanodiamonds in caenorhabditis elegans. Nano letters, 10(9):3692{3699, 2010. Edward K Chow, Xue-Qing Zhang, Mark Chen, Robert Lam, Erik Robinson, Houjin Huang, Daniel Scha er, Eiji Osawa, Andrei Goga, and Dean Ho. Nanodiamond therapeutic delivery agents mediate enhanced chemoresistant tumor treatment. Science translational medicine, 3(73):73ra21{73ra21, 2011. Susana Q Lima and Gero Miesenb ock. Remote control of behavior through genetically targeted photostimulation of neurons. Cell, 121(1):141{152, 2005. Edward S Boyden, Feng Zhang, Ernst Bamberg, Georg Nagel, and Karl Deisseroth. Millisecond-timescale, genetically targeted optical control of neural activity. Nature neuroscience, 8(9):1263{1268, 2005. Jacob G Bernstein, Paul A Garrity, and Edward S Boyden. Optogenetics and thermogenetics: technologies for controlling the activity of targeted cells within intact neural circuits. Current opinion in neurobiology, 22(1):61{71, 2012. Daniel E Bath, John R Stowers, Dorothea H ormann, Andreas Poehlmann, Barry J Dickson, and Andrew D Straw. Flymad: rapid thermogenetic control of neuronal activity in freely walking drosophila. Nature methods, 11(7):756, 2014. Eric A Newman and Peter H Hartline. The infrared" vision" of snakes. Scienti c American, 246(3):116{127, 1982. Fumika N Hamada, Mark Rosenzweig, Kyeongjin Kang, Stefan R Pulver, Alfredo Ghezzi, Timothy J Jegla, and Paul A Garrity. An internal thermal sensor controlling temperature preference in drosophila. Nature, 454(7201):217{220, 2008. Shijia Chen, Cindy N Chiu, Kimberly L McArthur, Joseph R Fetcho, and David A Prober. Trp channel mediated neuronal activation and ablation in freely behaving zebra sh. Nature methods, 13(2):147{150, 2016. Gopalakrishnan Balasubramanian, Philipp Neumann, Daniel Twitchen, Matthew Markham, Roman Kolesov, Norikazu Mizuochi, Junichi Isoya, Jocelyn Achard, Johannes Beck, Julia Tissler, et al. Ultralong spin coherence time in isotopically engineered diamond. Nature materials, 8(5):383{387, 2009. Jocelyn Achard, Vincent Jacques, and Alexandre Tallaire. Cvd diamond single crystals with nv centres: a review of material synthesis and technology for quantum sensing applications. Journal of Physics D: Applied Physics, 2020. Ettore Bernardi, Richard Nelz, Selda Sonusen, and Elke Neu. Nanoscale sensing using point defects in single-crystal diamond: recent progress on nitrogen vacancy center-based sensors. Crystals, 7(5):124, 2017. Edlyn V Levine, Matthew J Turner, Pauli Kehayias, Connor A Hart, Nicholas Langellier, Raisa Trubko, David R Glenn, Roger R Fu, and Ronald L Walsworth. Principles and techniques of the quantum diamond microscope. Nanophotonics, 8(11):1945{1973, 2019. A Dr eau, M Lesik, L Rondin, P Spinicelli, O Arcizet, J-F Roch, and V Jacques. Avoiding power broadening in optically detected magnetic resonance of single nv defects for enhanced dc magnetic eld sensitivity. Physical Review B, 84(19):195204, Chang S Shin, Claudia E Avalos, Mark C Butler, David R Trease, Scott J Seltzer, J Peter Mustonen, Daniel J Kennedy, Victor M Acosta, Dmitry Budker, Alexander Pines, et al. Room-temperature operation of a radiofrequency diamond magnetometer near the shot-noise limit. Journal of Applied Physics, 112(12):124519, 2012. 25 Erik Bauch, Swati Singh, Junghyun Lee, Connor A Hart, Jennifer M Schloss, Matthew J Turner, John F Barry, Linh Pham, Nir Bar-Gill, Susanne F Yelin, et al. Decoherence of dipolar spin ensembles in diamond. arXiv preprint arXiv:1904.08763, Michael Sean Grinolds, Sungkun Hong, Patrick Maletinsky, Lan Luan, Mikhail D Lukin, Ronald Lee Walsworth, and Amir Yacoby. Nanoscale magnetic imaging of a single electron spin under ambient conditions. Nature Physics, 9(4):215{219, G De Lange, ZH Wang, D Riste, VV Dobrovitski, and R Hanson. Universal dynamical decoupling of a single solid-state spin from a spin bath. Science, 330(6000):60{63, 2010. Piotr Szank  owski, Guy Ramon, Jan Krzywda, Damian Kwiatkowski, et al. Environmental noise spectroscopy with qubits subjected to dynamical decoupling. Journal of Physics: Condensed Matter, 29(33):333001, 2017. Boris Naydenov, Florian Dolde, Liam T Hall, Chang Shin, Helmut Fedder, Lloyd CL Hollenberg, Fedor Jelezko, and J org Wrachtrup. Dynamical decoupling of a single-electron spin at room temperature. Physical Review B, 83(8):081201, 2011. Helena S Knowles, Dhiren M Kara, and Mete Atatu re. Observing bulk diamond spin coherence in high-purity nanodiamonds. Nature materials, 13(1):21{25, 2014. Q-Y Cao, P-C Yang, M-S Gong, M Yu, A Retzker, Martin B Plenio, C Muller,  N Tomek, B Naydenov, LP McGuinness, et al. Protecting quantum spin coherence of nanodiamonds in living cells. Physical Review Applied, 13(2):024021, 2020. Mariusz Mr ozek, Daniel Rudnicki, Pauli Kehayias, Andrey Jarmola, Dmitry Budker, and Wojciech Gawlik. Longitudinal spin relaxation in nitrogen-vacancy ensembles in diamond. EPJ Quantum Technology, 2(1):22, 2015. Mariusz Radtke, Ettore Bernardi, Abdallah Slablab, Richard Nelz, and Elke Neu. Nanoscale sensing based on nitrogen vacancy centers in single crystal diamond and nanodiamonds: achievements and challenges. Nano Futures, 3(4):042004, Nir Bar-Gill, Linh My Pham, Chinmay Belthangady, David Le Sage, Paola Cappellaro, JR Maze, Mikhail D Lukin, Amir Yacoby, and Ronald Walsworth. Suppression of spin-bath dynamics for improved coherence of multi-spin-qubit systems. Nature communications, 3(1):1{6, 2012. Analia Zwick, Gonzalo A Alvarez, and Gershon Kurizki. Maximizing information on the environment by dynamically controlled qubit probes. Physical Review Applied, 5(1):014007, 2016. Matthias M Muller,  Stefano Gherardini, Nicola Dalla Pozza, and Filippo Caruso. Noise sensing via stochastic quantum zeno. Physics Letters A, 384(13):126244, 2020. Christian Arenz, Daniel Burgarth, Paolo Facchi, Vittorio Giovannetti, Hiromichi Nakazato, Saverio Pascazio, and Kazuya Yuasa. Universal control induced by noise. Physical Review A, 93(6):062308, 2016. Ruslan I Stanika, Idalis Villanueva, Galina Kazanina, S Brian Andrews, and Natalia B Pivovarova. Comparative impact of voltage-gated calcium channels and nmda receptors on mitochondria-mediated neuronal injury. Journal of Neuroscience, 32(19):6642{6650, 2012. Juliette Humeau, Jos e Manuel Bravo-San Pedro, Ilio Vitale, Lucia Nunez, ~ Carlos Villalobos, Guido Kroemer, and Laura Senovilla. Calcium signaling and cell cycle: progression or death. Cell calcium, 70:3{15, 2018.

Journal

Quantum Physics – arXiv (Cornell University)

Published: Jun 1, 2020

There are no references for this article.