Temperature estimation of a pair of trapped ions.

We apply estimation theory to a system formed by two interacting trapped ions. By using the Fisher matrix formalism, we introduce a simple scheme for estimation of the temperature of the longitudinal vibrational modes of the ions. We use the ions interaction to effectively infer the temperature of the individual ions, by optimising the interaction time evolution and by measuring only over one of the ions. We also investigate the effect of a non-thermal reservoir over the inference approach. The non-classicality of one of the ions vibrational modes, introduced due to a squeezed thermal reservoir, does not directly affect the inference of the individual temperatures, although allowing the modes to be entangled. To check actual experimental conditions, we analyze the temperature inference under heating due to surface-electrode noise.

Introduction

In the quantum mechanical context, the temperature of a system is a nonlinear function of the density operator, so it cannot directly correspond to a quantum observable—one has to indirectly estimate its value by measuring another observable. This indirect procedure for temperature estimation implies in an additional uncertainty for the measured value, which should be handled appropriately. Therefore, any strategy aimed to determine the temperature of a quantum system reduces to a parameter estimation problem[1-3]. Quantum theory of estimation (QET) provides a formal framework in order to optimise the inference procedure by minimising the additional uncertainty[4,5].

Techniques of quantum parameter estimation have been devoted to estimate the temperature in the context of quantum thermodynamics[6-9] and also for technological applications in many different branches of Science—ranging from Material Sciences to Biology and Medicine[10,11]. A central purpose is to employ a controlled quantum system (with view to applications in low temperature measurement[12]) to explore the thermodynamics in the regime of small-scale physics[13-16], where quantum effects become predominant[17-19]. Several quantum systems have been used to estimate very low temperature such as Bose-Einstein condensates[20], ultracold lattice gases[21], trapped ions[22-24], single-qubit[6], to mention a few.

In this work, we address the issue of the vibrational degrees of freedom temperature estimation in a system of two interacting trapped ions. Trapped ion systems constitute an outstanding platform for observing fundamental physical phenomena and has surfaced as one of the central potential architectures for implementing scalable quantum computation[25,26] and quantum simulation[27-29]. In those systems, the ions are trapped by oscillating electromagnetic fields in linear or surface-electrode Paul Traps and are held for manipulations at cryogenic temperatures. Usually, fluctuations in the electromagnetic fields introduce perturbations in the equilibrium state of each distinct ion, causing inhomogeneity in the ions temperatures and spurious heating. Therefore, to access the individual ions temperatures is fundamental for the establishment of optimal conditions for the ions manipulation. Our interest here is to investigate how precisely we are allowed to estimate the local temperatures of the vibrational degree of freedom of the two ions simultaneously. For this purpose, we employ the multiparameter Fisher informationto investigate the parameter estimation accuracy[30-32], where denotes the outcome of a measurement, is the conditional probability of measuring if the value of the parameter under consideration is . From the practice standpoint, increasing the Fisher information of the system tends to increase the maximum precision that can be obtained by an estimation scheme. Mathematically, this relationship is quantified by the Cramér-Rao lower bound[33,34],where is the estimator of the unknown parameter . The Cramér-Rao lower bound is saturated asymptotically by the optimization of the elements of the Fisher information matrix via a suitable choice of all its dependent parameters. Quantum parameter estimation purpose in quantum thermometry is to use a measurement over a quantum probe system, to infer the temperature of the reservoir it is immersed. For that purpose, an optimisation over the measurement operation is usually employed in order to find the maximal overall Fisher information (and consequently the measurement with the minimal dispersion). However, the procedures employed, are sometimes cumbersome and quite generally do not bring any information on how and which observable is to be measured[12]. Such a difficulty is even worse in the situation when systems of continuous variables are involved[16,35,36]. Here we take a more pragmatical approach, instead of optimising the measurement operator we employ measurements accessible in actual ion experiments[37,38] for the calculation of the Fisher information. We employ the detection of one of the ions vibrational mode phonon number, which is accessed experimentally, through the observation of the asymmetry between the red and blue motional sidebands of hyperfine Raman transitions[39-42]. The phonon number detection on one of the ions (say ion 1) is then employed for simultaneous inference of the temperature of both ions.

In what follows, firstly we introduce the theoretical model and we show how to describe the ionic vibrational modes. We propose a joint estimation scheme, where a single probe state is used to estimate the temperatures of the two ions ( and ) with a single projective measurement corresponding to the number of excitations in the first ion vibrational degree of freedom. In particular, we calculate the Fisher information as a function of the temperature of the ions and the interaction parameter g. We analyze the Fisher information of both ions in two scenarios—First, we consider the ions innitially in equilibrium with two individual thermal baths with distinct temperatures. Secondly, we investigate how a change in the vibrational state of one of the ion affects the performance of our approach. For that we consider a non-thermal bath in equilibrium with one of the ions, and also include the effect of spurious heating originating from the surface-electrode noise.

Two interacting charges

Ion trapping has a long and rich history, from its origins[43,44], to the usage as a potential platform for quantum computation[25-29,45]. The ions can be confined in free space with oscillating electromagnetic fields supplied by nearby electrodes[45], known as Paul Traps[44]. Originally constituted by hyperbolic electrodes, with time they were substituted, for more practical purposes, by surface-electrode Paul Traps[46]. Those are microfabricated planar linear traps, where a pseudo-potential allow ions to be held at distances of ~100 µm above the electrode surface. A typical configuration represented in Fig. 1, following the design in Ref.[37], shows alternating DC and AC electrodes, which depending on the experimental purposes allow strong confinement in the transversal direction and controllable confinement along the longitudinal direction (x in this paper). This system was previously employed to trap and investigate the interaction of two ions[37,38]. Given the trapping conditions, after cooling, the motional degree of freedom of the ions can be seen as coupled quantized mechanical oscillators held in separate locations. The direct coupling, through the mutual Coulomb interaction of two ions, was used and probed. Controllable coupling between quantized mechanical oscillators is quite useful for investigation and simulation of several properties observed in Nature. In particular, in Ref.[37], Be ions are held 40 µm apart by trapping potentials (See Fig. 1), and are both Doppler-cooled to K (thermal average number of ). After side-band cooling one of the ions to K (thermal average number of ), its population is probed, observing asymmetry between the red and blue motional sidebands of hyperfine Raman transitions. This ion population shows a signature of the population exchange, and therefore the interaction between ions.

Figure 1

Scheme of two interacting trapped ions with the same charge q. Each ion is confined in a potential well along of coordinate x and separated by a distance d. Schematic diagram based on Ref.[37], where the trap electrodes are placed in a single plane of a chip-based trap[46]. The RF and DC electrodes are placed to form a pseudopotential minimum above the surface, with the longitudinal confinement and multiple regions of confinement guaranteed by DC end caps, alternating with RF voltages, as depicted by the Control electrodes.

Here we focus on the detailed derivation and investigation of temperature inference for a similar situation to the one discussed in[37,38]. Let us consider two interacting trapped ions having the same charge q and mass m. Each ion is confined in a potential well along of coordinate x and separated by a distance d, as shown in Fig. 1, accordingly with refs.[37,38].

Considering only the coordinates along the longitudinal separation of the ions, the Hamiltonian that describes this system reads[37,38]where is the frequency of oscillation of each ion due to the trapping potential, and and are the momentum and position of the ion , respectively. The term is the ions electrostatic interaction energy given by[47]where . By expanding (4) in powers of , and considering small relative displacements we can rewrite the Coulomb interaction energy approximately asIntroducing the operatorsandthe total system Hamiltonian becomeswhere are the ions trap frequencies deviated by a small frequency shift due the ions interaction, and . In the interaction picture obtained by the unitary transformation , where corresponds to the trapped ions free evolution, the Hamiltonian becomesand the terms oscillating at frequencies and are averaged to zero for typical scales of time and can be disregard. The ions interaction Hamiltonian reduces to the simple formwith the coupling constant . It is important to note that the interaction (10) can be effectively controlled by detunig the individual ions motional frequencies, as discussed in Ref.[37]. Moreover, this coupling Hamiltonian is similar in form to that one for two radiation modes interacting through a beam-splitter. It is well known that a beam-splitter can only entangle two light modes if at least one of them is non-classical[48-50]. Therefore, it is interesting to investigate whether the non-classicality of one of the states affects the process of temperature inference. For that, in the next section, we develop a formalism that allows to tackle the dynamics of the two ions vibrational modes independently of the initial equilibrium conditions. Since we focus on the specific experimental conditions reported in Refs.[37,38], it is interesting to fix some quantities used along the remainder of this paper, for the reader reference, as presented in Table 1. Moreover, squeezing in the vibrational mode of trapped ions was reported in several occasions, since[41,42] (See Ref.[51] for an excellent account on squeezing in the motional states of ions.). Particularly in Ref.[52] it was reported a squeezing parameter as high as , which is well above the range of , employed here.

Table 1

Parameter used along this article, following experimental values employed in Ref.[37].

Ω	d	g	T  (n¯)	dn¯_i/dt=Γ,i=1,2
2π×4 MHz	40 µm	π×3.1 kHz	1.42-5.32×10^-4K(0.35-2.3)	1.885 quanta/ms

From left to right: is the ion vibrational angular frequency, d is the mean distance between ions, g is the ions coupling constant, T is the temperature with respective thermal average number inside parenthesis, and is the heating rate due to surface-electrode noise.

Bipartite Gaussian states

All states easily accessed experimentally for our purposes, are Gaussian. Therefore, here we develop the formalism more appropriate to investigate the properties of the ionic vibrational degrees of freedom. A two-mode bipartite quantum state is Gaussian if its symmetric characteristic function[53] is given by , where is a displacement operator in the four-vector -space: , , andwhere () and () are annihilation (creation) operators for party 1 and 2, respectively. is a covariance matrix with elements , which can be decomposed in four block matrices,where and are Hermitian matrices containing only local elements, while is a matrix representing the correlation between the two parties and explicitly are written asPositivity and separability for bipartite Gaussian quantum states have been largely investigated[54-56]. Besides the requirement of the uncertainty principlethere is a necessary and sufficient condition, which must be satisfied for separable Gaussian states[56]under a partial phase space mirror reflection, , withThe physical positivity criterion (14) applies only if[56]respectively, with , , and . Similarly the separability condition (15) writes explicitly into (17) andThe ions vibrational modes are assumed as initially uncoupled, so that in (12) , and they are prepared in special local states to be discussed later. However, the effect of the interaction (10) is to correlate the two modes. This can be seen as the following Bogoliubov operation B:where is the density operator for the joint output state. The output symmetric characteristic function is given byNow with the help of Eqs. (11–13), , with , since . Thuswhere , and analogously to (12), can be block decomposed withFor the specific unitary transformation induced by the time-independent two-modes coupling (10) , and , where t is the time variable. Typically, the separability criterion (15) is recast in terms of the local symplectic invariants, belonging to the group, , , , and . Using the , , , and . A state is not entangled if, and only if[55]We will employ the violation of the bound (28) as an entanglement indicator, after we define, the initial motional states.

As mentioned before, we shall consider the case where the ions are initially uncorrelated so the total density operator factorizes aswhere and are the reduced density operators of ions 1 and 2, respectively. We assume that the ion 1 is initially in equilibrium with a thermal reservoir at temperature , while ion 2 can be in equilibrium with a thermal or non-thermal (for future purposes) reservoir at temperature . For that, we introduce a controllable squeezing parameter so that the ion 2 can be prepared in a thermal squeezed state. Thus, the reduced density operator for the ions 1 and 2 in a coherent state representation are written asandrespectively, whereis the one-mode squeezing operator, , and the average thermal photon numbers are given bywhere and correspond respectively to the equilibrium temperature of the ion-1 and ion-2.

The time evolution of the density operator with Hamiltonian (10) occurs with the following expressionwhere , andTo understand the correlation induced by (34), we plot in Fig. 2 the bound (28) as a function of and the squeezing parameter r for the ion 2, for a fixed temperature of K (), and K (). Only for , is that the vibrational mode 2 is non-classical, being a necessary (but not sufficient[56]) condition for the evolution (34) to entangle the vibrational modes of the ions. As it is expected, the maximal entanglement occurs at (or ), which corresponds to a 50:50 beam-splitter. Only for is that the system can be entangled.

Figure 2

Entanglement of the ions logitudinal vibrational modes as a function of the squeezing parameter of ion 2, and the interaction . Only for , is that the vibrational mode 2 is non-classical, allowing that the evolution (34) entangle the ions (represented by the negative values of the vertical axis). Maximal entanglement occurs at (or ), corresponding to a 50:50 beam-splitter.

Fisher information

The trapping conditions require several steps of cooling and electromagnetic confinement of ions. Typically, the ions coupling parameter g is very small ( kHz, here), and therefore the temperature of the system must reach cryostatic regimes of a few hundred microkelvins. Therefore, an inhomogeneity is expected and the individual ions equilibrium temperatures may differ. Further experimental limitations constrains the direct access to all ions in the trap[37], as was suggested by the inference approach[24]. Our purpose here is to use the phonon number detection (over one of the ions only) as a way to infer simultaneously the temperature differences of both ions. In particular, the estimation of the temperatures of the ions can be implemented by choosing the projective measurement corresponding to the number of phonons in the longitudinal vibrational mode of the ion 1, as described in Fig. 3. It shows our estimation scheme, where the dynamics of two trapped ions is governed by the unitary operation (34).

Figure 3

Scheme for the indirect estimation of the temperatures of two interacting trapped ions. After some interaction time between the two ions, a energy measurement corresponding to the number of phonons is performed on the 1-ion, followed by post-processing for the estimation of parameters and .

In the temperature estimation procedure, it is essential to have access to the ion 1 longitudinal vibrational mean number. This is a standard experimental detection procedure, well described in the literature[39,51]. An individual ion vibrational mode phonon number is accessed experimentally, through the observation of the asymmetry between the red and blue motional sidebands of hyperfine Raman transitions[39]. For Be ions, for example, during the Raman cooling, the internal and hyperfine states are accessed by two counter-propagating Raman beams. In this process, transitions from to (where is the vibrational state of the ion) are induced cyclically, together with dissipative repumping, resulting in the cooling of the ion 1. It also allows to access the probabilities of measuring , and , after excitation with red or blue sideband, respectively[39,51]. The average motional occupation is given by , where . In that fashion the population of the vibrational state of ion 1 is measured. We can calculate the Fisher information corresponding to the parameters and via the number of excitations in the longitudinal vibrational degree of freedom of ion 1.

For the number of phonon measurements, the Fisher information for ion 1 and ion 2 is computed aswhere is the probability distribution of the occupation number obtained by measuring the energy of ion 1.

The probability distribution of k phonons is given bywherewithThe probability distribution given by Eq. (36) involves a Gaussian integral that can be evaluated exactly so that the probability distribution is given as followswhereandThe probability distribution does not depend on the squeezing phase since it willl not appear in the term .

Now using the series expansionthe derivatives in relation to A are easily evaluated, so the probability distribution reduces toThe series in Eq. (44) can be identified with the Hypergeometric[58] function so we can writeIn the case that , with both modes in thermal states, the probability distribution reduces to the simple expressionand the Fisher information elements of Eq. (35) are obtained in a closed form given bywhere is the Kronecker delta. To obtain the values of the Fisher information for non null squeezing () we set the sums up to in Eqs. (35) in order to guarantee convergence in the numerical calculations.

In Fig. 4, the Fisher information , and evolution with , is represented for two thermal reservoirs (for ) at , and , corresponding respectively to the experimental values of and in Ref.[37], and for a non-thermal reservoir (for ). We can see by the corresponding curves of that the temperature is best inferred at the instants where , for . This occurs because at those instants the initial equilibrium state of the ion 1 has recurred[59]. In contrast, as can be seen in , at , for the state of ions 1 and 2 has been swapped, by the characteristic of evolution (34), and therefore it is not surprising that those are the best instants for inference of temperature of the vibrational mode of ion 2—it occurs because when the states are interchanged[59], the detection of the ion 1 vibrational mode is, in fact, the detection over the state of ion 2. Simultaneous inference of both ions temperatures occurs when , for , as depicted by . Even when a non-thermal state is involved, the approach is successful for the inference of the ion 2 temperature, as we can see in the same figures by varying r, the same characteristics hold. The effect of squeezing, however, is to improve the entanglement, and therefore it turns out that the distribution of points where the inference is optimal is spread, meaning that there is a slight advantage in using entanglement as a way to improve the accessibility at distinct instants of time for inference.

Figure 4

Fisher information , and as a function of parameter , for (), (), and MHz[37,57].

In Fig. 5, the Fisher information , and are depicted as functions of the temperature for (similar results are observed for odd multiples of ), and different values of r, for a fixed temperature of ion 1 at . This result indicates a precision loss for the estimation of temperature of the ion 1, when we increase the r (this is reproduced in the graphs of and in Fig. 5. The same cannot be said for estimating only the temperature , which is improved with an increase in r (see graphs of in Fig. 5) under the same conditions, working as a super probe in the non-classical regime , for this maximally entangled situation for . In fact, the maximal precision is still observed at , as depicted in Fig. 6, in contrast to and , which are null for any we see that reaches optimal values, which however is disturbed by the non-thermal reservoir feature.

Figure 5

Fisher information , and as a function of parameter for , (), and MHz[37,57].

Figure 6

Fisher information as a function of parameter for , and MHz[37,57]. At , all Fisher information becomes independent of since the ions have exchanged their states, and particularly, .

Lastly, we remark on the observed flat response regions observed for all , in Figs. 5 and 6, for very low temperatures. To understand the reason behind that, we expand the exact expression for for , from Eq. (47) for at . Thus, at very low the Fisher information elements behave aswhere and , with , are constants dependent on the other parameters in(47). Since the exponentials rapidly decrease as goes to zero, it becomes clear the plateau in the Fisher information elements in Figs. 5 and 6 for small values of . The physical reason is that at small , the vibrational energy of ion 1 is not sufficient to disturb the inference information as given by the . This is because the quantum state of the vibrational degree of motion of ion 1 is not considerably disturbed, and therefore in the dispersion in the measurement of the average number of phonons of ion 1 is dependent only on the quantum statistical properties of its state. This is reinforced by the fact that and remain close to zero at the same interval that the plateau occurs in , meaning that no information about the vibrational energy of ion 2 can be obtained.

Heating due to surface-electrode noise

Fluctuations in the trap parameters may introduce perturbations in the equilibrium state of each distinct ion as we discussed in the previous section. In addition, a thermal heating rate due noisy electrical potentials on the trap electrode surfaces is typically present (For an excellent account, please see Ref.[46]). The residual fluctuation in the electric fields couple to the motional degrees of freedom of the ions, inducing transitions between vibrational state. In[37], e.g., a heating of 1.885 quanta per milisecond was observed, while in[38] it was of 1.3 quanta per milisecond. The heating effect has been successfully treated through a master equation formalism[46], or equivalently through the Heisenberg-Langevin equation formalism, as in Ref.[60]. In any case, the individual ions can be considered in thermal equilibrium with their individual reservoirs, whose temperatures are steadily increasing (in Thermodynamics it is formally stated that the system follows an equilibrium path, being sequentially in contact and in thermal equilibrium with an infinite number of reservoirs, each one with an infinitesimally different and increasing temperature.) and then suffering the individual heating at the same rate. This was the procedure successfully considered in[37,38] to describe the heating rate observed experimentally. The Langevin approach turns out to effectively disturb the individual equilibrium bosonic operators for the vibrational degrees of freedom in the following form,where are noise operators satisfyingbeing the decay constant induced by the effective memoryless reservoirs, with mean occupation number , modeling the effects of surface-electrode noise. We have assumed the same decay constant for both the ions, in order to simplify the analysis.

The solution of the system of coupled equations (49) can be obtained in terms of the noise, and modes operators at aswhich, in turn, can be used to evaluate all the elements of the covariance matrix in Eq. (29). For the ion 1, using the properties of the noise operators, from Eqs. (50)–(53), we haveNow, Eq. (55) clearly shows that, with time, the exchange of excitations between the vibrational modes of the ions ceases to exist and , reaching equilibrium with the effective reservoir. However, this timescale is unrealistic and too far from actual experimental conditions[46], since the effective reservoir emulates the effect of spurious surface-electrode noise. Here we are interested in timescales consistent with the experimental conditions of Ref.[37], where the damping in the excitations exchange is negligible. In addition, the damping mechanism introduced also takes into account the decoherence, through the decay of the other correlations appearing in Eqs. (25–27). In this timescale, where the decay is very slow compared to the excitation exchange timescale between ions (), we can approximate the exponential decay to its first relevant order in each term in Eq. (55) to getand the effective effect of noise is to introduce a linearly time dependent heating at rate in the system, as experimentally observed in Ref.[37]. The same consideration, is obviously valid for mode 2. At the same timescale considered, the correlations are not affected by the heating noise, and therefore the motional state is not affected by decoherence. Other decoherence mechanisms, such as phase noise, may appear at a near timescale of a few miliseconds, particularly in experiments involving coupling between the motional and internal electronic states. In Refs.[61,62], e.g., the observed decoherence time due to dephasing is around 5 ms. However, this is about 6.3 times larger than the timescale of , corresponding to ms, of the actual experimental condition of Ref.[37]. Therefore, at the timescale we are considering, decoherence does not take place, although the entanglement between the two vibrational modes is affected, since the global state becomes more classical with the heating, as we will discuss latter on. Again, at larger timescales the heating would be the dominant effect, and would already compromise the present strategy. In any case, when the internal electronic states are considered, one should employ a combined strategy for investigation of decoherence and heating, as in ref[61], for example.

Therefore, Eq. (38), now readsFigure 7 shows the population of mode 1, Eq. (57), as function of time and some values of squeezing parameter r including a heating rate of 1.885 quanta per milisecond estimated in[37] and is in perfect accordance with the experimental and simulated results obtained in Fig. 3 of that reference for . We see an exchange of excitations between the two motional populations with a steady linear increase due to the heating. For , there is a pronounced rise in the population, particularly at , when the two states are swapped, which is due to the squeezing in mode 2. The increase in population is more evident at , as the initial population of mode 2 increases with .

Figure 7

Population of mode 1, as function of for (), (), MHz, kHz, and a heating rate quanta per milisecond, such that [37,57].

In this situation, the Fisher information must be calculated following the same procedures of the previous section, but using Eq. (57) instead. For , now the Fisher information elements of Eqs. (33) readIn Fig. 8, similarly to Fig. 4, we plot the Fisher information , and evolution with , for the thermal and non-thermal reservoirs. We can see that the heating reduces the Fisher information with time, meaning that it introduces noise (dispersion) in the temperatures inference, as the individual temperatures increase. However, the noise effect is not substantial to hinder the inference of the ions temperatures at the heating rate of 1.885 quanta per milisecond. In fact, we observed that even for a heating rate of ten times the experimentally observed rates, the procedure is still useful to infer the temperature of mode 2 by measuring the excitations of mode 1. The same conclusions related to the spreading of the distributions and the advantage of the inference procedure without noise in “Fisher information” section, remain valid with noise. The only pronounced effect is to diminish the previously available periodicity for inference. Since the heating effect is more pronounced as the time goes on, the amplitude of the Fisher information peaks decrease with time, as from Eq. (58), which means that the heating effect introduces a dispersion, increasing with time in the inference of the temperatures. Also, we should remark, that the profiles of the Fisher information as the temperature of ion 2 is changed (Fig. 5) is not significantly changed, as it evaluates the best regimes for the inferred temperatures. This further confirms what we explained above.

Figure 8

Fisher information , and as a function of parameter for (), (), MHz, kHz, and a heating rate quanta per milisecond ( )[37,57].

The only remaining question is how the entanglement between the vibrational modes of the two ions is affected by the heating effect. In Fig. 9 we see that it is severely affected by the heating effect. This occurs, despite the fact that only the populations of the modes, and not the correlations, are affected by the linear heating at the investigated timescale. The heating increases the modes populations linearly with time, while the other correlations are not affected. Therefore, the global system state becomes more classical. For the state to become entangled, it would be necessary to have higher values of r at disposition initially. Since the temperature estimation here presented, is not strictly dependent on the entanglement between the vibrational modes (available due to the squeezing), it is clear why, the Fisher information in Fig. 8 is not so severely affected. As a last remark we should mention that if other additional decoherence mechanism were present at the considered timescale, such as the dephasing we discussed before, the coherences would decay, and the global state would become even more classical, as if the squeezing parameter were smaller. This however only affects the entanglement and do not affect the strategy at the relevant timescale.

Figure 9

Entanglement of the ions logitudinal vibrational as a function of the squeezing parameter of ion 2, r, and the interaction , at a heating rate quanta per milisecond (). The entanglement is severely affected by the heating as the global state becomes more classical.

Concluding remarks

In summary, we have explored the interaction between trapped ions for simultaneous inference of their vibrational temperatures. For entanglement analysis and temperature estimation, we take ion 1 with a thermal initial state, and ion 2 with the initial state to be thermal and non-thermal (due to quadrature squeezing). We use the interaction to effectively infer the temperature of the individual ions, by optimizing the interaction time evolution. It was observed that while the estimation of the temperature of ion 2 is slightly improved when the initial state is non-classical, the entanglement is not an essential ingredient, since what really matters is the ions state exchange, resulting from their interaction. Simultaneous measurements are often difficult to implement, in particular when associated with motional degrees of ions. The most practical way is to estimate the temperature by measuring each ion. Here, we chose to measure the ion that does not contain a possible squeezing, making it as a classical channel probe, contrary to what it has in the literature, where they use detectors with quantum properties for estimating parameters of classical systems[7,63-65]. To appropriately consider actual experimental conditions we have included a heating effect which is observed to occur due to surface-electrode noise. For that we developed a Langevin equation, considering two uncorrelated noise sources. The noise induces a linear (with time) heating effect, in accordance with experimental observations. Remarkably, although the noise, in the regime employed, decreases the entanglement of the joint vibrational degrees o freedom, it does not disturb considerable the Fisher information. This means that the method employed for temperature inference is useful, being robust against noise. Finally, we should remark that in trapped ion systems, thermal, displaced, squeezed and other arbitrary states were already generated, and therefore can be identified[41,42]. Although we considered only thermal, and squeezed thermal reservoirs, the same procedure would work for any arbitrary state. Moreover, the present strategy can be extended to remotely access other statistical properties of the ions.