New Ferroelectric Capacitance Model frmc from Ramtron Corporation in SmartSpice
Introduction
Ramtron International Corporation has developed a ferroelectric capacitance model with a new concept of double distributions of domain reversal voltages.
The model is now available in the public domain. It has been implemented in SmartSpice and can be invoked by setting LEVEL=6 in the capacitance model card. The previous model is also implemented in SmartSpice and can be invoked by setting LEVEL=5 in the capacitance model card.
Features
The ferroelectric capacitor is treated as a nonlinear capacitor, except the biases in the ferroelectric materials are reversed at reversal voltages with double distributions. The model has been tested under various combinations of input voltage waveforms. Compared to the previous model (LEVEL=5), this model is much more accurate in ferroelectric hysteresis loops and voltage pulse responses, its simulation speed is at least six times faster than the previous interpreter model, and the temperature dependence is included. Furthermore, the parameters in frmc can be easily extracted.
Ferroelectric Capacitor Device
Syntax form
Cxxx n1 n2 mname <IC=val> <L=val> <W=val> <IP=val> <A=val>
| Cxxx: | Capacitor element name. Must begin with a ?C?. |
| n1,n2: | Positive and negative terminal node names, respectively. |
| mname: | Model name (must be a ferroelectric capacitor model). |
| IC: | Initial voltage across the ferroelectric capacitor (V). Default is the applied voltage. |
| L: | Capacitor length (m). Default is 0.0. |
| W: | Capacitor width (m). Default is model parameter DEFW. |
| IP: | Initial polarization. Default is 0.8. |
| A: | Capacitor area (m2). Default is 1.0E-12. |
Output Device Variables
| Parameters | Description |
| cap (capval, ceff, cf) | Device capacitance |
| c (i, if) | Total current (DC, transient, AC) |
| p (power) | Power |
| vf | Device voltage |
| q | Charge |
| cq | Current due to charge |
Ferroelectric Capacitor Model
Syntax form
.MODEL mname C LEVEL=6 <parameter=value>
| C: | Specifies a capacitance model. |
| mname: | Model name. |
| parameter: | Any model parameter name. |
New Model Parameters
| Parameter | Description | Units |
Default |
| VMAX | Maximum voltage | V |
5.0 |
| PMAX | Maximon bias at Vmax | uC/cm2 |
29.85 |
| KPMAX | Temperature coefficient at Pmax | -3.47e-2 | |
| VSAT | Saturated model point | 2.0 | |
| VCR | Minimum voltage for domain switching | V |
0.5 |
| AS1 | curve fitting parameter | 2.5236e-1 | |
| BS1 | curve fitting parameter | 1.1090 | |
| CS1 | curve fitting parameter | 1.3287 | |
| AS2 | curve fitting parameter | c 5.1668e-2 | |
| BS2 | curve fitting parameter | 2.3961e-1 | |
| CS2 | curve fitting parameter | -1.4109e-1 | |
| DS0 | curve fitting parameter | 5.086e-1 | |
| AU1 | curve fitting parameter | 4.4077e-1 | |
| BU1 | curve fitting parameter | 5.1064e-1 | |
| CU1 | curve fitting parameter | 1.6479 | |
| AU2 | curve fitting parameter | 1.3386e-2 | |
| BU2 | curve fitting parameter | -1.9351e-2 | |
| CU2 | curve fitting parameter | -1.7674e2 | |
| DU0 | curve fitting parameter | 3.5208e-2 | |
| KAS1 | curve fitting parameter | -1.0107e-3 | |
| KBS1 | curve fitting parameter | 1.4579e-3 | |
| KCS1 | curve fitting parameter | -3.0895e-3 | |
| KAS2 | curve fitting parameter | 4.8937e-4 | |
| KBS2 | curve fitting parameter | 6.2892e-4 | |
| KCS2 | curve fitting parameter | 7.2631e-3 | |
| KDS0 | curve fitting parameter | 9.7153e-5 | |
| KAU1 | curve fitting parameter | -9.5819e-4 | |
| KBU1 | curve fitting parameter | 2.6618e-3 | |
| KCU1 | curve fitting parameter | -6.1257e-3 | |
| KAU2 | curve fitting parameter | 3.7970e-4 | |
| KBU2 | curve fitting parameter | 2.7597e-4 | |
| KCU2 | curve fitting parameter | 1.0525 | |
| KDU0 | curve fitting parameter | -3.4559e-2 |
Simulation Results
Before showing results obtained with the new ferroelectric capacitance model frmc (LEVEL=6), some simulation results are presented to show the limitations of the previous model (LEVEL=5). The previous model LEVEL=5 is successful in modeling the major hysteresis loops under continuous sinusoid and triangular input voltages but, the major hysteresis loops become unrealistic when the capacitor is driven by voltage pulses. To illustrate this the input waveform shown in Figure 1 was applied to a circuit containing a ferroelectric capacity shown in Figure 2. When the capacitor models were LEVEL=5 large jumps can be observed at three points (0V, +/- 5V) on the hysteresis loops shown in Figure 3, This means that the charge on the ferrocapacitance still increases even when the input voltage is set at 5V. This phenomenon is not observed in real ferro-electric capacitors. This discrepancy is due to the incorrect modeling for the domain reversal behavior.

Figure 1. Input pulse voltage waveform

Figure 2. Sawyer-Tower circuit at hystersis loop measurement

Figure 3. Hysteresis loops under input pulse voltage waveform
The ferroelectric capacitance model frmc has been modified to consider relatively large voltage oscillations which are just across the coercive field. The new model has passed tests on large circuit simulations with various voltage oscillations.
The new ferroelectric capacitance model (LEVEL=6) has corrected the previous problems with stepped waveforms. The same circuit shown in Figure 2 was applied to test the new model.
First a sinusoidal input waveform was applied which resulted in the output waveform shown in Figure 4. The hysteresis loops under the sinusoidal voltage at 1 MHz are correct. Next an input voltage step shown in Figure 5 was applied. Figure 6 shows the hysteresis loops obtained which are now clearly correct. Further, an unwanted peak was also present in the input pulse of Figure 5. Its influence can be observed on the hysteresis loop of Figure 6 which is again correct. Finally, to test the new model to the extreme the input voltage waveform shown in

Figure 4. Hysteresis loops under sinusoidal voltage form at 1MHz frequency.

Figure 5. Input pulse voltage waveform with widths of several ns

Figure 6. Hysteresis loops under pulse voltage waveform
Figure 7 was applied. The output waveform shown in Figure 8 shows the correct hysteresis behavior even under these extreme conditions.

Figure 7. Input voltage with unwanted oscillations

Figure 8. Hysteresis loops from unwanted oscillations
Conclusion
The previous ferroelectric capacitance model (LEVEL=5) unrealistically simulates ferroelectric subloops and voltage pulse responses. Thus, it is dangerous to apply for simulations at low voltage or voltage pulse responses.
The new ferroelectric capacitance model frmc (LEVEL=6), based on double distributions of domain reversal voltages, provides accurate simulation for ferroelectric hysteresis loops and sub-loops, transient responses to short voltage pulses (with widths in the nanosecond range), and temperature behavior of ferroelectric capacitors. The advantages can be summarized as:
- The new model is accurate and reliable for any applied voltage waveforms
and voltage levels. The applied voltages can be short pulses at the
nanosecond scale.
- The simulation speed is at least six times faster than the previous
interpreter model.
- The new model does not need gap capacitors and resistors. From the
user's point of view, it is a two-terminal device; thus, it is easy
to use.
- The parameters in frmc can be easily extracted.
- The temperature dependence has been included in the model.
References
- Ramtron International Corporation documents