Research Article | Open Access

J. M. Ávila-Ruiz, A. Moscoso-Mártir, E. Durán-Valdeiglesias, L. Moreno-Pozas, J. de-Oliva-Rubio, I. Molina-Fernández, "Six-Port-Based Architecture for Phase Noise Measurement in the UWB Band", *Journal of Electrical and Computer Engineering*, vol. 2014, Article ID 646738, 10 pages, 2014. https://doi.org/10.1155/2014/646738

# Six-Port-Based Architecture for Phase Noise Measurement in the UWB Band

**Academic Editor:**Mohamed Helaoui

#### Abstract

A six-port-based frequency discriminator for phase noise measurement is proposed. This circuit makes use of a delay line discriminator configuration, thus not requiring reference oscillator. Furthermore, the use of a six-port network allows an extremely simple and completely passive solution well suited for low power oscillator measurement. A detailed study of the architecture is performed including the system noise sources. Besides, a prototype of the proposed six-port based delay line frequency discriminator is evaluated. Phase noise measurements of a commercial RF VCO are performed and compared with the results obtained with commercial metrology equipment with good agreement.

#### 1. Introduction

The high data rate transmission requirements of current radiofrequency standards demand efficient and robust modulation techniques. One of these modulation techniques is Orthogonal Frequency Division Multiplexing (OFDM), which is widely used in wireless communications due to its good performance against noise and multipath effects. On the other hand, OFDM receivers are very sensitive to local oscillator phase noise, which produces adjacent channel crosstalk due to reciprocal mixing [1, 2]. This adjacent channel crosstalk increases bit error rate (BER) [3, 4]. Phase noise effects in radiofrequency communications are the reason why feedback systems, such as phase locked loops (PLL), have been used for phase noise reduction for 80 years [5, 6]. Moreover, a great amount of new on chip low noise oscillators designs are carried out every year [7, 8].

An outstanding issue, as important as the design of low phase noise oscillators, is the accurate measurement of their phase noise. Modern oscillators with low phase noise level require metrology systems with even more demanding specifications as well as measurement methods to compensate their limitations [2, 9]. Several methods have been defined for phase noise measurement, depending on the characteristics of the system to be measured. As a main distinction, these methods can be first classified depending on whether they require a reference oscillator or not. This differentiation is important, since reference oscillator phase noise will increase the noise floor of the measurement system. In order to increase the instrument sensitivity, ultralow phase noise reference oscillators are needed. Therefore, high performance systems that rely on reference oscillators are more expensive than other instrument systems.

The group of phase noise measurement techniques that make use of a reference oscillator comprises the direct spectrum measurement and the PLL techniques. The direct spectrum measurement is the most straightforward method (Figure 1(a)). The spectrum of the oscillator under test, from now on the device under test (DUT), can be directly measured with a spectrum analyzer or previously down-converted to an intermediate frequency (IF) using a reference oscillator and a mixer. The PLL method makes use of a mixer as a phase detector (Figure 1(b)). Then, the output of the phase detector is low-pass-filtered to eliminate the higher frequency terms, retaining only the baseband component. The phase noise of the DUT is extracted from the spectrum of the low-pass filtered output. In this method, a PLL is needed to lock the reference oscillator to the oscillator under test.

**(a)**

**(b)**

The need of a reference oscillator is removed by the second group of phase noise measurement methods, which make use of a frequency discriminator. The most used one is the delay line frequency discriminator. A typical phase noise measurement system based on a delay line discriminator is shown in Figure 2. This method makes use of a mixer-based phase detector, as in the aforementioned PLL method. However, in this architecture, the RF signal from the oscillator under test is compared in the phase detector with a delayed version of itself. A phase shifter must be included in order to keep the quadrature condition between both signals at the phase detector inputs [10]. Even though this method gets rid of the reference oscillator, it has an important drawback: the complex circuitry needed to control the phase shifter and keep the system into the quadrature condition.

Recently, a new architecture for phase noise measurement based on delay line discriminator has been proposed [11]. In this circuit, a dual phase detector architecture and a 90-degree hybrid allow the detection of the phase noise from the in-phase and quadrature signal components without quadrature adjustment. Even though this architecture needs neither reference oscillator nor quadrature adjustment circuitry, it has an important disadvantage: its input-referred phase noise floor considerably increases at lower input power levels [11]. A good solution to improve the phase noise measurement of this new architecture is the utilization of a six-port demodulator (see Figure 3). Six-port demodulators can work with low input power levels and their calibration has been widely studied [12]. The six-port architecture has been previously implemented in [13], for W-band phase noise measurements, but using a reference oscillator, obtaining good accuracy in phase noise measurement of one-port and two-port systems.

The architecture proposed in this work makes use of the idea of dual phase noise detection based on a delay line discriminator [11] and the six-port architecture used in [13] in order to implement a compact phase noise measurement system. This architecture uses neither reference oscillator nor a phase shifter for quadrature adjustments. Moreover, the existence of powerful calibration strategies for six-port demodulators [12] and their good performance for low input power levels make the use of this structure advantageous against traditional mixers.

The rest of this paper is organized as follows. In Section 2, the six-port-based delay line discriminator architecture is described. In Section 3, the noise contributions and their effects are detailed. In Section 4, measurement setup and results are discussed. Finally, Section 5 is devoted to conclusions.

#### 2. The Six-Port-Based Delay Line Discriminator

The proposed phase noise measurement architecture unifies the delay line method with the use of a six-port demodulator as a phase detector to measure the phase noise from the in-phase and quadrature signal components. The block diagram of the six-port-based delay line discriminator is shown in Figure 4. This figure includes the noise sources, , referred to the detector output, which will only be used in the noise analysis carried out in Section 3.

In order to analyze the system response in the phase noise detection, a RF input signal is considered. The time-domain signal is where is the phase noise, is the amplitude noise, is the rms voltage, and is the microwave carrier frequency. Using phasor nomenclature, the rms complex envelope from (1) is

The relation between the time-domain signal and the complex envelope is

Assuming that and are ergodic noise signals, and knowing that is very small compared with , the power of the RF signal is

Now that the input signal has been defined, the six-port based delay line discriminator will be analyzed. In Figure 4, it can be seen that the input RF signal is equally split in two halves. One half travels through a delay line before being fed in one of the inputs of the six-port receiver and the other one is directly sent to the second input of the receiver. If an ideal 3 dB power splitter and a lossless delay line are assumed, the time-domain complex envelopes at the six-port inputs are where is the delay of the line. For the sake of simplicity, the product is defined as . This phase shift depends on the carrier frequency and is constant as long as keeps constant. As the carrier frequency changes, this phase shift will change with a period of Hertz.

For the six-port block, an ideal passive network with four 90-degree hybrids and a 90-degree phase shift is used [14]. Therefore, insertion loss of 6 dB per path is considered and all are balanced in amplitude and phase. With these assumptions, the signals at the six-port outputs are

Using the concept of centers of the six-port architecture, [15], all expressions, (6) to (9), can be written together as where the centres are , , and .

The instantaneous voltage of each output is obtained from a diode power detector. Assuming square law detectors, the voltage at each power detector output is where is the power detector sensitivity and is the instantaneous power at the input of the detector . It is supposed that the four detectors are identical; hence the same is assumed for all of them. Equation (11) can be rewritten as where is a complex variable that relates the two input signals of the six-port, analogous to the received symbol in classical six-port based receivers [16]. This variable has the phase noise information in its phase as it can be seen in the following definition of : where is the instantaneous phase of . For very small , is proportional to the derivative of the input phase noise:

If (12) is operated, knowing that for all , the expression of the four outputs of the power divider is obtained:

These outputs must be combined in order to obtain scaled versions of the in-phase and quadrature components of : where is a complex variable proportional to depicted in the complex plane in Figure 5. As the phase noise is measured from the phase shift of and the amplitude noise only affects the module of this variable, as seen in (13), the phase noise measurement is independent of the AM noise of the DUT.

The phase in both complex variables, and , is the same, as shown in (18). Thus, it is straightforward to obtain the phase noise information from the phase of :

If the DC component of (19), , is eliminated, a time variant phase signal is measured. This phase signal is proportional to the frequency fluctuations at the system input and can be represented in the frequency domain as a power spectral density:

From (14), a relation between the power spectral density of and the power spectral density of the phase noise, , is obtained: where is the angular frequency deviation from the RF carrier frequency and is the power spectral density of the phase fluctuations, . This relation can be simplified for . In this case, the exponential in (21) can be approximated by the first Taylor series expansion term; that is, , becoming

As can be seen in (22), the RF frequency variations in the system input can be detected in the six-port delay line discriminator output. In fact, the time-domain variation of the detected phase signal () is proportional to the time-domain frequency variations of the RF oscillator under test. Thanks to the relation between phase and frequency, it is easy to obtain the phase fluctuations at the input from these frequency variations. The previous analysis has been carried out assuming lossless blocks. However, following [16], it can be shown that, assuming well balanced couplers and dividers, element losses do not affect the phase of . A more detailed analysis must be done if the imbalances of the different elements in the discriminator architecture are considered. Moreover, as it will be shown in Section 3, the losses of the system increase the noise floor of the system, partially degrading the sensitivity of this measurement method.

#### 3. Noise Analysis

Although in the previous section the effects of the internal noise sources of the measurement setup have not been considered, this noise is an important factor to take into account. The internal noise limits the measurement system noise floor. The measurement system is passive, with the exception of the amplifiers used after the power detection and the digital processing circuitry. Therefore, two different noise sources are included in the analysis, depending on the noise origin: (i) the noise from the power detector (diode plus detection circuitry) and (ii) the noise generated in the data acquisition system. The rest of the noise sources (i.e., thermal noise due to losses in couplers and dividers) are negligible because at these components signal level is well over its thermal noise floor. All the noise sources are referred to the power detector output, as seen in Figure 4, and their spectral power density is measured in . Since both internal noises are uncorrelated, the total internal noise at each detector output is where is the noise from each data acquisition channel and is the noise from each detection block. So, if the internal noise is included in the analysis, the signals at the detectors outputs are with for , respectively. It is assumed that the noise signals at every detector output are identical and uncorrelated; that is, for all . The in-phase and quadrature components of defined in (16) and (17) are rewritten as

The noisy complex variable is depicted in Figure 6. It can be seen how the aforementioned internal voltage noise sources, , generate a phase fluctuation in the phase of . With the assumption of small angle, the arc of the angle in the circumference of radius can be approximated as . Finally, the phase shift introduced by the internal noise is obtained calculating the angle of the arc as

Equation (26) can be expressed in power spectral densities as follows: where is the power spectral density of the phase fluctuations of due to the internal noise and is the power spectral density of the internal voltage noise referred to the power detector output.

It can be seen from (27) that the noise contribution of the system depends on the input power and the sensitivity of the power detectors. Therefore, for higher input power level, the influence of the internal noise decreases.

In order to know the limitations of the measurement method, it is necessary to obtain the minimum discernible value of phase noise at the system input, that is, the phase noise floor. This noise floor can be calculated by referring the power spectral density of the internal voltage noise, , which is usually referred at the detector output (as seen in Figure 4), to the system input. From (27), the relation between the power spectral density of the voltage noise at the detector output and the power spectral density of the phase fluctuations, , can be obtained. Then, with the help of (22), the phase noise floor can be calculated as where is the power spectral density of the input-referred phase noise floor, .

As shown in (28), the phase noise floor at the measurement system input can be improved in three ways: (i) increasing the gain of the system by using a larger delay line, with higher ; (ii) increasing the input power and the sensitivity of the power detectors; and (iii) directly reducing the internal noise. The following subsection will be devoted to discussing the third option.

##### 3.1. Internal Noise Reduction

As it is shown in (23), there are two important sources of internal noise. On the one hand, the noise generated in the data acquisition system depends on the number of quantification bits, the internal electronics, and the input voltage range. These characteristics are determined for each digital acquisition system. On the other hand, the noise generated by the power detection and baseband circuitry can be minimized if a deep analysis and a careful design and implementation of these blocks are performed. The schematic of the power detection and baseband circuitry used in this work is depicted in Figure 7.

The resistances and are the resistances of the amplifier and they set the amplifier gain, . They are related as . For low frequencies, the matching network at the power detector input, with a shorted stub, works as a short circuit. Therefore, applying the operational amplifiers topologies noise analysis [17], the voltage noise at the output is where is the Boltzmann constant, is the temperature (in Kelvin), is the video resistance of the diode, and and are the equivalent input current noise and the equivalent input voltage noise of the amplifier, respectively. If the approximation is fulfilled, (29) simplifies to

From (30), it can be concluded that there are two important noise sources in this circuit (): the diode detector and the operational amplifier. The election of a diode with a lower video resistance () and the operational amplifiers with low equivalent current and voltage noise is one of the keys to minimize internal noise and, therefore, to reduce the noise floor of the measurement setup.

#### 4. Simulation and Measurement Results

Simulations with the design environment Advanced Design System (ADS) from Agilent Technologies have been done. For these simulations, each functional block has been modeled depending on their nature. The power divider and the six-port have been characterized by full-vectorial 3D electromagnetic simulation using Ansys HFSS. In doing so, their parameters have been calculated. Then, they have been used to model both devices as multiport parameters defined components into ADS circuit schematic. A 15 ns delay line (LabFlex290 coaxial cable) has been modeled from the manufacturer parameters. These three blocks are passive and linear so they are fully defined from their parameters. For the rest of the blocks, diodes and amplifiers, ADS models based on the manufacturer specifications have been used. The diode model is based on HSCH-9161 zero bias diode from Avago Tech. and the amplifier has been modeled from the specifications of LMH6622 amplifier from Texas Instrument.

A prototype has been fabricated in order to assess the simulated results, which is shown in Figure 8. The prototype comprises a Wilkinson power splitter with 3.4 dB insertion loss at 5 GHz, a 3.779 m LabFlex290 delay line from EMC-RFLabs with a delay of 15 ns and 1.5 dB insertion loss at 5 GHz, an ultrawideband six-port circuit with state-of-the-art characteristics in amplitude and phase imbalance [14], four power detector blocks that use a zero bias diode HSCH-9161 from Avago Technologies, and an operational amplifier LMH6622 from Texas Instruments. The gain of the amplifiers is set to 20 dB (10 V/V). The output of each power detector is sampled and digitalized with a data acquisition board NI-6115 from National Instruments.

After the digitalization, the signals from the power detectors are processed with MatLab in order to obtain the phase noise spectrum. A simplified diagram of the mathematical process is shown in Figure 9. To take advantage of the full dynamic range of the data acquisition board, two measurements are done, with DC and AC coupling, in order to get the DC component and the time-varying part of the (16) and (17) signals. From these two measurements, the complete signals are built and the subsequent complex variable is obtained from (18). From the phase of , the power spectral density is calculated with (20) and the phase noise power spectral density at the input is obtained applying the FFT, as described in (22). The applied FFT uses a rectangular window of 200 ms and a frequency rate of 10 MHz. These parameters are enough for the measurement of phase noise between 100 Hz and 1 MHz from the RF carrier. As the six-port junction imbalances are very low [14], no calibration is required and the signals are directly obtained by subtracting the corresponding power readings, as seen in (16) and (17).

In this section, four measurements are shown. They are compared with simulation, theoretical calculations, and measurements using commercial instrumentation. The measurements are (i) the static response of the system, (ii) the internal noise, (iii) the dynamic response to a RF input signal with a frequency modulation, and (iv) phase noise measurement of a commercial RF oscillator.

##### 4.1. Static Response

To characterize the static response, a slow carrier frequency sweep, between 5 and 5.2 GHz, is done. This sweep is repeated for input powers between −10 dBm and −4 dBm. The results of the diagrams are shown in Figure 10, for simulation and measurements. In both graphics and for a concrete input power, an amplitude difference between and components that flatten the diagram can be seen. Besides, the module of the circle is not constant for different cycles. Both effects are because the system has not been calibrated.

**(a)**

**(b)**

##### 4.2. Internal Noise

The mathematical estimation of the internal noise has been separated in two parts, one for the power detector circuitry and one for the data acquisition system. The noise of the data acquisition board depends on the minimum voltage step, , and this minimum step depends at the same time on the quantification bits and the voltage range of the board. The data acquisition board NI-6115 uses 12 bits of quantification and its minimum input voltage range is ±200 mV, so

The specified noise spectral density in the sampling board datasheet is , where is the sampling frequency. For NI-6115, . Therefore, the board noise level is −147.9 dBV/Hz.

The power detector circuitry noise is calculated from (30). For this block, = 680 Ω and = 33 Ω resistances have been used in the amplifier in order to achieve a 20 dB gain. The diode HSCH-9161 video resistance is . The internal noise sources of the amplifier are and . From (30), the amplifier expected added noise is −136.5 dBV/Hz. Thus, noise limit is set by the power detector circuitry and can be estimated with (23) to a value of −136.5 dBV/Hz approximately.

As it can be seen in Figure 11, the measured total noise, at one data acquisition channel input, is 5 dB higher at 1 kHz (due to the effect of flicker noise) and 1 dB lower than the theoretically expected for frequency offsets higher than 100 kHz. The measured internal noise is approximately −137.5 dBV/Hz for the four channels.

##### 4.3. FM Demodulation

The dynamic characterization of the system is done by the detection of a frequency modulated RF signal. Using the N5181A vector signal generator from Agilent Technologies, a carrier of 5 GHz and −6 dBm is modulated with a frequency modulation, , of 30 kHz and frequency deviation, , of 100 kHz. This modulated RF signal is introduced in the six-port delay line discriminator. Since the main function of the system is the frequency fluctuations detection, a phase shift signal, , proportional to the baseband signal is expected. This phase shift signal is shown in Figure 12. It has a time period of 33 *μ*s, so its frequency is 30 kHz as expected and its maximum phase deviation, , is ±9.7 mrad. The proportionality factor of the system is

The calculation of the frequency deviation is as follows: where is the maximum phase deviation of . Therefore, for a value of 9.7 mrad, the measured frequency deviation is 102.9 kHz, very close to the original 100 kHz.

##### 4.4. Phase Noise Measurement

For the phase noise measurement, a voltage controlled oscillator (VCO) from HITTITE (HMC587LC4B) has been chosen as the device under test. This VCO features 4 dBm output power and it has an output carrier frequency range from 5 GHz to 10 GHz. Its phase noise is −95 dBc/Hz and −118 dBc/Hz for offset frequencies of 100 kHz and 1 MHz, respectively.

Two different phase noise measurements have been performed for two different VCO output power levels. Hence, 6 dB and 10 dB attenuators have been used at the VCO output to change the VCO output power, obtaining −2 dBm and −6.6 dBm, respectively. The VCO output frequency has been set at 5 GHz. For both setups, the VCO phase noise has been measured with the prototype shown in Figure 8 and with the signal source analyzer E5052B from Agilent technologies. The E5052B measurements have been performed using the internal PLL of the E5052B (that the E5052B uses in order to phase noise measure); an averaging factor of 20 and an IF gain of 50 dB have been used. For the measurements with the prototype, no PLL has been used (free running oscillator) and an averaging factor of 20 has been used.

The results of the measurements for both configurations are depicted in Figures 13 and 14. In each figure, the input-referred phase noise floor of the measurement system (voltage noise at the diode outputs referred to the input with (28)) is also plotted for comparison. The measurements results with the six-port delay line discriminator architecture (without calibration) are in very good agreement with the measurement results of commercial systems. Moreover, the proposed method is able to measure phase noise from sources with low carrier power without preamplification as well. Although the measurement system noise floor appears to be very high, as it can be seen in (28), it can be reduced using a longer delay line as it was done in [11].

#### 5. Conclusions

In this work, a new phase noise measurement architecture has been proposed. This architecture makes use of a six-port frequency discriminator. This circuit is based on a delay line discriminator configuration, thus not needing a reference oscillator. Most of the elements of the proposed architecture are passive (only the baseband operational amplifiers need DC supply), highly reducing its power consumption. The use of a six-port in spite of RF mixers allows phase noise measurement of low power oscillators without RF amplification.

A detailed study of the architecture has been performed including the system noise sources. Besides, the gain and sensitivity of the system have been calculated. In order to assess the proposed architecture, a prototype of the proposed six-port-based delay-line frequency discriminator has been fabricated. Several tests have been carried out to evaluate its performance. The phase noise of a commercial RF VCO has been measured and compared with the results obtained with a commercial phase noise measurement system. Even though no calibration algorithms have been used, the measured results of the proposed architecture are in good agreement with the commercial system ones, which show the validity of the proposed architecture for the characterization of RF VCOs.

#### Conflict of Interests

There is no conflict of interests regarding the publication of this paper.

#### Acknowledgment

This work has been funded by the Consejería de Economía, Innovación, Ciencia y Empleo of the Junta de Andalucía (Project P09-TIC-5268).

#### References

- B. Razavi, “Challenges in the design of frequency synthesizers for wireless applications,” in
*Proceedings of the IEEE Custom Integrated Circuits Conference*, pp. 395–402, May 1997. View at: Google Scholar - Agilent Technologies, “Product Note 11729C-2—Phase noise characterization of microwave oscillator,” 2007, http://www.agilent.com/. View at: Google Scholar
- S. Wu and Y. Bar-Ness, “OFDM systems in the presence of phase noise: consequences and solutions,”
*IEEE Transactions on Communications*, vol. 52, no. 11, pp. 1988–1996, 2004. View at: Publisher Site | Google Scholar - D. Petrovic, W. Rave, and G. Fettweis, “Effects of phase noise on OFDM systems with and without PLL: characterization and compensation,”
*IEEE Transactions on Communications*, vol. 55, no. 8, pp. 1607–1616, 2007. View at: Publisher Site | Google Scholar - F. M. Gardner,
*Phaselock Techniques*, Wiley Interscience, 3rd edition, 2005. - D. Banerjee,
*PLL Performance, Simulation and Design Handbook*, National Semiconductors, 4th edition, 2006. - E. Ayranci, K. Christensen, and P. Andreani, “Enhancement of VCO linearity and phase noise by implementing frequency locked loop,” in
*Proceedings of the International Conference on Computer as a Tool (EUROCON '07)*, pp. 2593–2599, Warsaw, Poland, September 2007. View at: Publisher Site | Google Scholar - A. Esswein, G. Dehm-Andone, R. Weigel, A. Aleksieieva, and M. Vossiek, “A low phase-noise SiGe colpitts VCO with wide tuning range for UWB applications,” in
*Proceedings of the 13th European Microwave Integrated Circuits Conference (EuMIC '10)*, pp. 452–455, October 2010. View at: Google Scholar - D. Owen, “Good practice guide to phase noise measurement,” National Physical Laboratory, Middlesex, UK, Tech. Note, 68, 2004, http://www.npl.co.uk. View at: Google Scholar
- K. P. Faith, F. Labaar, A. L. Lance, and F. G. Mendoza, “Frequency switching speed measurements using a delay line discriminator,”
*IEEE Transactions on Instrumentation and Measurement*, vol. 37, no. 4, pp. 620–625, 1988. View at: Publisher Site | Google Scholar - H. Gheidi and A. Banai, “Phase-noise measurement of microwave oscillators using phase-shifterless delay-line discriminator,”
*IEEE Transactions on Microwave Theory and Techniques*, vol. 58, no. 2, pp. 468–477, 2010. View at: Publisher Site | Google Scholar - F. Wiedmann, B. Huyart, E. Bergeault, and L. Jallet, “A new robust method for six-port reflectometer calibration,”
*IEEE Transactions on Instrumentation and Measurement*, vol. 48, no. 5, pp. 927–931, 1999. View at: Publisher Site | Google Scholar - S. O. Tatu, E. Moldovan, S. Affes, B. Boukari, K. Wu, and R. G. Bosisio, “Six-port interferometric technique for accurate W-B and phase-noise measurements,”
*IEEE Transactions on Microwave Theory and Techniques*, vol. 56, no. 6, pp. 1372–1379, 2008. View at: Publisher Site | Google Scholar - A. Moscoso-Martir, I. Molina-Fernandez, and J. M. Avila-Ruiz, “Ultra-wideband multilayer six-port junction with common phase error cancellation,” in
*Proceedings of the 42nd European Microwave Conference (EuMC '12)*, pp. 17–20, October 2012. View at: Google Scholar - G. F. Engen, “The six-port reflectometer: an alternative network analyzer,”
*IEEE Transactions on Microwave Theory and Techniques*, vol. 25, no. 12, pp. 1075–1080, 1977. View at: Google Scholar - A. Moscoso-Martir, I. Molina-Fernandez, and A. Ortega-Moñux, “Signal constellation distortion and BER degradation due to hardware impairments in six-port receivers with analog I/Q generation,”
*Progress in Electromagnetics Research*, vol. 121, pp. 225–247, 2011. View at: Google Scholar - Texas Instruments, “Application report SLVA043B—Noise analysis in operational amplifier circuits,” 2007, http://www.ti.com/. View at: Google Scholar

#### Copyright

Copyright © 2014 J. M. Ávila-Ruiz et al. This is an open access article distributed under the Creative Commons Attribution License, which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.