Mixed FEM for quantum nanostructured solar cells
Introduction
Quantum dots (QDs) are tiny nanocrystals made of a semiconducting material, which are buried into a piezoelectric matrix. They have wide range applications in microelectronic, optoelectronic devices and solar cell technology [1], [2]. There are also plenty applications in bioengineering and biomedical studies [3]. The QD solar cells are called as third-generation solar cells. In last years their properties are intensively studied to achieve a high-efficiency solar energy conversion [4]. Recently, Zheng et al. [5] have published a progress towards quantum dot solar cells with enhanced optical absorption in a comprehensive review.
A reliable analysis of elastic and electric fields are crucial to the design and fabrication of such structures [6], [7], [8], [9]. These fields are created in this system due to the lattice mismatch between the QD and matrix [10]. The intrinsic strains have remarkable effects on electronic and optical properties of quantum nanostructure cells. Strains change the interatomic distances and consequently energy levels and bonding of electrons. Then, knowledge of induced strains coupled with electric fields is essential for both understanding and the design of such photovoltaic components [11]. The QDs can work at extreme thermal conditions, therefore, Patil and Melnik [12] applied the classical thermo-piezoelectricity theory for QDs under stationary boundary conditions. Recently, Duc et al. [13] have developed the first analytical approach to investigate the nonlinear dynamic response and vibration of imperfect rectangular nanocomposite multilayer organic solar cell subjected to mechanical loads using the classical plate theory. However, classical continuum mechanics neglects the interaction of material microstructure and the results from it are size-independent. In a nanocomposite structure the size effect should be considered.
In many research works authors deal with uniform initial misfit strains. However, consideration of nonuniform misfit strains would be of great importance since the lattice mismatch occurs only at the interface of the QD and matrix. Inside of the QD it should be considered a functional variation with vanishing misfit at the center. Recently, Bishay et al. [14] have analysed the grading of the material composition as well as the lattice mismatch strain between the QDs and the host matrix. High-density arrays of quantum dots are extensively investigated as components of high performance solar cells. At high densities of QDs the interaction between them should be considered. However, in literature there are missing such analyses. Recently, Rashidinejad and Naderi [15] have studied electro-elastic fields in quantum nanostructure solar cells with consideration of the inter-nanostructure couplings and geometrical effects. One can observe a strong influence of the inter-nanostructure couplings significantly on the induced electro-mechanical fields. Authors [14] applied the classical theory of electro-elasticity for this problem.
The size of quantum dots is only several nanometers. Thus, the material length scale is comparable with the dimensions of the QDs and classical electro-elasticity should not be applied to analyze QD systems. In the classical theory the interaction with the material structure is neglected and results are size-independent. The first attempt to apply the gradient theory to a QD system is given by authors [16]. The gradient theory of thermo-piezoelectricity there is based on the simplified Aifantis theory [17], [18], [19] with one scaling parameter.
This paper presents, for the first time, the application of the gradient theory of piezoelectricity for 3D analyses of QDs with functionally graded lattice mismatch between the QD and matrix. The influence of distances between QDs is investigated too. For this purpose the FEM is developed and applied. In gradient theory it is needed to use C1-continuous elements to guarantee the continuity of the problem variables and their derivatives at the element boundaries. Since it is not an easy task to develop C1 continuous elements in 3D problems, the mixed FEM formulation is developed here. The C0 continuous interpolation is independently applied for displacement and displacement gradients. The independent interpolation requires to satisfy the kinematic constraints between strains and displacements. A collocation at some cleverly chosen internal points in elements is applied [20]. The present approach reduces the number of DOFs with respect to the Lagrangian mixed FEM formulation [21].
Numerical results for high-density arrays of quantum dots with functionally graded lattice mismatch between the QD (InAs) and matrix (GaAs) are presented, and the influence of the size effect parameter is discussed.
Section snippets
The gradient theory for electro-mechanical fields
A periodic distribution of QDs (InAs) in the matrix (GaAs) is considered and a representative volume element (RVE) for the quantum nanostructure solar cell is illustrated in Fig. 1.
The misfit strains induce electro-mechanical fields in this nanostructure solar cell. This eigenstrain is created by the lattice mismatch, which are occured only at the interface of the QD and matrix. Inside of the QD functional variations should be considered with vanishing misfit at the center. A quasistatic
The mixed finite element method (FEM)
The principle of virtual work is applied to derive the FEM equations for a boundary value problem in gradient theory of piezoelectricity. It means that variation of the electro-elastic enthalpy has to be equal to the work of the external generalized forces on generalized displacementswhere overbar is used for prescribed boundary values.
The order of derivatives in governing equations in gradient theory of piezoelectricity is higher
Numerical results
The representative volume element (RVE) in Fig. 1 is assumed to have a cubic geometry with side length of b = 40 nm, and a cubic QD with side length of is embedded in it. The top surface of the InAs QD is 4 nm from the top surface of the matrix. Material properties of the QD correspond to InAs [30]. The substrate is made of GaAs, whose properties are
The four side faces of the matrix are fixed along their normal
Conclusions
The 3D mixed finite element model has been developed for QDs with functionally graded lattice mismatch between the QD and the surrounding matrix. The periodic QD array is replaced by a representative volume element (RVE), which is analyzed by the gradient theory of piezoelectricity. Since governing equations in the gradient theory contain higher order derivatives than in conventional approaches, a continuous approximation of strains is required. The mixed FEM with the C0 continuous
Acknowledgment
The authors acknowledge the support of the Slovak Science and Technology Assistance Agency registered under number APVV-18-0004, and VEGA-2/0046/16.
References (31)
Quantum dot solar cells
Physica E
(2002)- et al.
Nonlinear dynamic response and vibration of nanocomposite multilayer organic solar cell
Compos Struct
(2018) - et al.
On the structure of the mode III crack-tip in gradient elasticity
Scripta Metall Mater
(1992) - et al.
Gradient elasticity in statics and dynamics: An overview of formulations, length scale identification procedures, finite element implementations and new results
Int J Solids Struct
(2011) - et al.
Stress concentrations in fractured compact bone simulated with a special class of anisotropic gradient elasticity
Int J Solids Struct
(2010) - et al.
Fracture mechanics analysis of size-dependent piezoelectric solids
Int J Solids Struct
(2017) - et al.
Quantum dot solar cells
Prog Photovoltaics Res Appl
(2002) - et al.
Potential clinical applications of quantum dots
Int J Nanomed
(2008) - et al.
Quantum dot solar cells
(2014) - et al.
Recent progress towards quantum dot solar cells with enhanced optical absorption
Nanoscale Res Lett
(2016)
InAs/GaAs pyramidal quantum dots: strain distribution, optical phonons, and electronic structure
Phys Rev B
Three-dimensional strain field calculations in multiple InN/AlN wurtzite quantum dots
J Appl Phys
Elastic and piezoelectric fields around a quantum dot: fully coupled or semicoupled model
J Appl Phys
Elastic and piezoelectric fields in substrates GaAs (001) and GaAs (111) due to a buried quantum dot
J Appl Phys
Quantum Dot Heterostructures
Cited by (5)
XFEM for multiphysics analysis of edge dislocations with nonuniform misfit strain: A novel enrichment implementation
2023, Computer Methods in Applied Mechanics and EngineeringMulti-physical modelling, design optimization and manufacturing of a composite dielectric solar absorber
2022, Composites Part C: Open AccessCitation Excerpt :Recently, a lot of interest has been directed at composite metal–insulator based absorbers, as they offer high solar selectivity, depending on the constituent materials for metal and insulator. A reliable analysis of the complex coupled physical fields are crucial to the design and fabrication of such composite structures [1–4]. Selective solar absorbers are being categorized into different types, in which each type of structure plays a different role in maximizing absorption and minimizing emittance for best efficiency [11–14].
The MLPG Method in Multiphysics and Scale Dependent Problems
2021, Mechanisms and Machine ScienceFinite Element Numeric Simulation of Organic Solar Cells with Gold Thin Film
2020, International Journal on Advanced Science, Engineering and Information Technology