Texture and anisotropy of hot-pressed h-BN matrix composite ceramics with in situ formed YAG

Textured hexagonal boron nitride (h-BN) matrix composite ceramics were prepared by hot-pressing using different contents of 3Y2O3−5Al2O3 (molar ratio of 3:5) as the sintering additive. During hot-pressing, the liquid Y3Al5O12 (YAG) phase showing good wettability to h-BN grains was in situ formed through the reaction between Y2O3 and Al2O3, and a coherent relationship between h-BN and YAG was observed with [010]h−BN//[111]YAG\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$${\left[ {010} \right]_{{\rm{h}} - {\rm{BN}}}}{\rm{//}}{\left[ {111} \right]_{{\rm{YAG}}}}$$\end{document} and (002)h-BN//(321)YAG. In the YAG liquid phase environment formed during hot-pressing, plate-like h-BN grains were rotated under the uniaxial sintering pressure and preferentially oriented with their basal surfaces perpendicular to the sintering pressure direction, forming textured microstructures with the c-axis of h-BN grains oriented parallel to the sintering pressure direction, which give these composite ceramics anisotropy in their mechanical and thermal properties. The highest texture degree was found in the specimen with 30 wt% YAG, which also possesses the highest anisotropy degree in thermal conductivity. The aggregation of YAG phase was observed in the specimen with 40 wt% YAG, which resulted in the buckling of h-BN plates and significantly reduced the texture degree.


Introduction 
Hexagonal boron nitride (h-BN) has a layered structure similar to that of graphite, where B and N atoms are along the a-and b-axis, and the c-axis was calculated to be 4.1 and 537 W/(m·K), respectively [8], together with its insulating property [9], making h-BN matrix composites the promising thermal management devices in electronics [10][11][12][13][14][15].
The strong sp 2 covalent bonds combining B and N atoms result in a low atomic diffusion rate even at high sintering temperatures, and the plate-like shape of h-BN grains usually form the "house of cards" structure with high porosity, making h-BN ceramics hard to be densified during sintering [16][17][18][19]. Introducing sintering additives composed of two or more oxides which can form liquid phases at high sintering temperatures and fill these pores between plate-like h-BN grains, such as Y 2 O 3 -Al 2 O 3 , Y 2 O 3 -SiO 2 , Yb 2 O 3 -MgO, Yb-Si-Al(N)-O, and La-Al-Si-O [20][21][22][23][24][25][26], has been proven to be an effective method to improve the sinterability of h-BN ceramics.
Compared with plate-like h-BN powders with large grain sizes, fine h-BN powders with more defects are easier to be densified due to their high specific surface area and defect energy [27,28]. However, the final gain size of h-BN ceramics prepared from fine h-BN powders is limited even with sintering additives which can form liquid phases during sintering [29][30][31], and thus h-BN ceramics with highly textured microstructures cannot be obtained. Plate-like h-BN powders with larger grain sizes and sintering additives which can form liquid phases are necessary for the fabrication of textured h-BN matrix composite ceramics with significant anisotropy.
In our previous study, the textured h-BN matrix composite specimen with an index of orientation preference (IOP) of −530 was prepared by hot-pressing using 20 wt% 3Y 2 O 3 -5Al 2 O 3 with a molar ratio of 3 : 5 as the sintering additive, which was transformed into YAG phase showing good wettability to h-BN grains [22]. In this study, different contents (10,20,30, and 40 wt%) of 3Y 2 O 3 -5Al 2 O 3 were introduced to prepare h-BN matrix composite ceramics by hot-pressing to investigate the influence of YAG contents on the texture and anisotropy in the mechanical and thermal properties of the composite ceramics.

Materials and methods
The raw materials were h-BN powders with a platelike morphology (~15 μm in diameter and ~0.3 μm in thickness, purity > 99%), Al 2 O 3 powders (~0.2 μm, purity > 99.9%), and Y 2 O 3 powders (~2 μm, purity > 99.9%). 3Y 2 O 3 -5Al 2 O 3 (molar ratio of 3 : 5) were used as the sintering additive. Four specimens, whose mass ratios of h-BN to the sintering additive were set as 9 : 1, 8 : 2, 7 : 3, and 6 : 4, were prepared as follows: Raw powders were mixed using Al 2 O 3 balls in ethanol on a pot-mill machine with a rotation speed of 80 rpm for 24 h and then dried on a hot plate set at 70 ℃, and then the mixed powders were loaded into a Ф50 mm cylindrical graphite die and hot-pressed at 1900 ℃ under a uniaxial pressure of 30 MPa in N 2 atmosphere for 1 h. Hotpressed h-BN matrix composite ceramics with 10 wt%, 20 wt%, 30 wt%, and 40 wt% YAG (it will be proven in Section 3.1 that Y 2 O 3 and Al 2 O 3 fully reacted during hot-pressing, forming YAG with a cubic crystal structure) were marked as BN-1YAG, BN-2YAG, BN-3YAG, and BN-4YAG, respectively.
Bulk densities of the hot-pressed specimens were measured by Archimedes method. Phase compositions and crystallographic orientation were analyzed by the X-ray diffractometer (XRD; X′PERT, Panalytical, the Netherlands). Further analyses on phases and interface microstructure were realized by the transmission electron microscope (TEM; Tecnai G2 F30, FEI, USA). Properties were measured along three different loading directions by the electronic universal testing machine (Model 5569, Instron, USA). The flexural strength was tested on the 3 mm × 4 mm × 36 mm sample bars by three-point bending with a span of 30 mm and a crosshead speed of 0.5 mm/min; the elastic moduli were tested using strain gauges stuck on the tensile surfaces of sample bars used for strength testing; the fracture toughness was tested on 2 mm × 4 mm × 20 mm single-edge-notched beams with a span of 16 mm and a crosshead speed of 0.05 mm/min. Fracture morphologies were observed using the scanning electron microscope (SEM; Quanta 200 FEG, FEI, USA). Coefficients of thermal expansion (CTEs) were measured on 5 mm × 5 mm × 25 mm bars using a dilatometer (DTL 402C, Netzsch, Germany). Specific heat capacity and thermal diffusivity were measured on 12.7 mm × (2-3) mm discs by a laser flash apparatus (LFA 427, Netzsch, Germany) and the thermal conductivity was calculated from the specific heat capacity and thermal diffusivity.
top surface (TS) perpendicular to the sintering pressure direction and side surface (SS) parallel to the sintering pressure direction of h-BN matrix composite ceramics hot-pressed using different contents of 3Y 2 O 3 -5Al 2 O 3 additive. All specimens are composed of two phases: h-BN and YAG, indicating that Y 2 O 3 and Al 2 O 3 with the molar ratio of 3:5 have fully reacted. They all show a much higher h-BN(002) diffraction peak for TS than that for SS and a higher h-BN(100) diffraction peak for SS than that for TS, which can be seen in the illustrations showing the complete (002) diffraction peaks at the top-right side of these XRD patterns in Fig. 1, indicating that all specimens possess highly textured microstructures with the c-axis of h-BN grains preferentially oriented parallel to the sintering pressure direction, which were attributed to the rotation of plate-like h-BN grains under the uniaxial sintering pressure in the YAG liquid phase environment. Figure 2 gives the relative densities, IOP, and Lotgering orientation factors f 00l of hot-pressed h-BN matrix composite ceramics with different YAG contents. The relative density increases with the increase of the YAG content and reaches the highest for BN-4YAG (92.36%), as the liquid YAG phase formed during hot-pressing acted as the pore filler between plate-like h-BN grains.
IOP and f 00l are used to quantitively characterize the texture degree of h-BN matrix composite ceramics. IOP was calculated according to Eq. (1) [32]:  I  I  I  I  I  I  I  I   I  I  I  I  I  I where hkl I and hkl I are the intensities of (hkl) diffraction peaks for TS and SS, respectively. For a hot-pressed h-BN matrix composite specimen, IOP = ±1 indicates that it consists of randomly oriented h-BN grains; IOP < −1 and IOP > 1 indicate that its h-BN grains are preferentially oriented with the c-axis parallel and perpendicular to the uniaxial sintering pressure direction, respectively. In this study, a negative IOP with a larger absolute value indicates a higher texture degree.
The f 00l was calculated according to Eqs.   was calculated from the measured XRD patterns for TS of hot-pressed specimens, while P 0 was calculated from the standard JCPDS card (PDF#34-0421) of h-BN. P and f 00l are always lower than 1. A higher f 00l closer to 1 indicates more (00l) planes on TS. When increasing the YAG content in hot-pressed h-BN matrix composite ceramics from 10 to 40 wt%, the absolute values of IOP and f 00l both increase first and then decrease dramatically, and reach the highest for BN-3YAG (IOP = −1778; f 00l = 0.99), which means that BN-3YAG has the highest texture degree. Figure 3 shows the bright field images and selected area electron diffraction (SAED) patterns for SS of BN-3YAG and BN-4YAG. These plate-like h-BN grains are highly oriented with their (002) planes perpendicular to SS in BN-3YAG, which is confirmed by the SAED pattern of area A in Fig. 3(a). The SAED pattern of area B shows a secondary diffraction pattern caused by h-BN below the YAG area, and there are a row of diffraction spots deviating from the marked (002) direction, which was caused by the orientation difference of adjacent h-BN grains. The arrangement of h-BN grains in BN-4YAG is not as orderly as that in BN-3YAG, which is demonstrated by the existence of both (1 10) planes (area A in Fig. 3(b)) and (002) planes (area B in Fig. 3(b)) of h-BN grains. The aggregation of liquid YAG phase formed during hot-pressing (area C in Fig. 3(b)) caused the buckling of plate-like h-BN grains, which significantly reduced the texture degree of BN-4YAG. Figure 4 shows the high-resolution transmission electron microscopy (HRTEM) image, fast Fourier transform (FFT) pattens of several square regions in the HRTEM image, and corresponding inverse fast Fourier transformation (IFFT) images around the interface between h-BN and YAG in BN-3YAG. Areas 1-4 in Fig. 4(a) are h-BN, the interface between h-BN and YAG, YAG, and the h-BN area infiltrated by YAG, respectively. Stacking fault was observed in h-BN ( Fig. 4(b)), which has been proven to commonly exist in hot-pressed h-BN ceramics due to the "multi-area co-growth" of h-BN grains in our previous study [29]. The IFFT image of area 2 in Fig. 4(a) displays the coherent interface between h-BN and YAG (Fig. 4(c)) formed due to the very close interplanar spacing between the (002) plane of h-BN and the (100) plane of YAG, where the (002) and (100) planes of h-BN are parallel to the (321) and (145) planes of YAG, respectively. The IFFT image of area 3 in Fig. 4(a) displays the crystal structure of YAG observed from the [111] direction ( Fig. 4(d)), where the hexagons composed of Al and Y atoms can be easily distinguished, while the secondary diffraction pattern in the FFT pattern of area 3 confirms the existence of h-BN below the well-crystallized YAG area. The FFT pattern of area 4 indicates that this area contains both h-BN and YAG, and its IFFT image shows that the YAG phase infiltrates into the h-BN grain from multiple directions (Fig. 4(e)), further demonstrating the good wettability of YAG to h-BN grains. Figure 5 shows the surface morphologies for TS and SS of hot-pressed h-BN matrix composite ceramics with different YAG contents. Many continuous small particles (< 0.5 μm) exist at the edges of plate-like h-BN grains, especially for BN-1YAG and BN-2YAG. Figure 6 shows the EDS spectra of points 1 and 2 in Fig. 5(d), where point 2 has much lower B and N peaks and much higher Y, Al, and O peaks compared with point 1, indicating that these small particles outside h-BN grains are YAG.
For BN-3YAG with the most highly textured microstructure, YAG particles were mainly observed on SS which contains most (100) planes, i.e., the side surfaces of plate-like h-BN grains. Such a result corresponds to the phenomenon that the diffraction peak intensities of YAG in the XRD patterns for SS are much higher than those for TS of BN-3YAG (Fig. 1(c)). For BN-4YAG, almost no YAG particles were observed, while many plate-like h-BN grains were buckled, which is consistent with the TEM observation in Fig. 3(b), indicating that YAG phase could be hidden under these buckled h-BN plates. Figure 7(a) shows the schematic diagrams of the fracture modes and the fracture morphologies along different loading directions (D1, D2, and D3) of hotpressed h-BN matrix composite ceramics with different YAG contents. For D1, the long axis of the sample bar is perpendicular to h-BN layers (parallel to the c-axis orientation) and the loading direction is parallel to h-BN layers, resulting in a relatively flat fracture morphology. For D2, the long axis of the sample bar is parallel to h-BN layers and the loading direction is perpendicular to h-BN layers, resulting in a tortuous crack propagation path. For D3, the long axis of the sample bar and the loading direction are both parallel to h-BN layers, generating a tortuous initial crack and a fracture morphology similar to that along D2.
The fracture morphologies of BN-3YAG show the most significant texture (Figs. 7(h)-7(j)), which is in consistent with its highest absolute values of IOP and f 00l . For BN-1YAG, the side of (002) plane of h-BN was observed on the fracture morphology along D1, and the (002) plane of h-BN was observed on the fracture morphology along D2 and D3 (Figs. 7(b)-7(d)). For BN-4YAG, many buckling behaviors of h-BN plates were observed (Figs. 7(k)-7(m)), which was caused by the aggregation of YAG phase (Fig. 3(b)) and significantly reduced its texture degree. Figure 8 gives the mechanical properties of hot-pressed h-BN matrix composite ceramics with different YAG contents. The schematic diagram of samples with different orientations in the bulk specimen for testing anisotropy in mechanical properties was shown in Fig. 8(a). The work of fracture (γ WOF ) was calculated from the load-displacement curves during fracture toughness testing according to Eq. (4) [34]:

2 Anisotropic properties
 is the total fracture energy during the fracture toughness testing which is determined as the area under the load-displacement curves, and A is the cross-sectional area of the unnotched part of the singleedge-notched beams. As the sp 2 covalent bonds combining B and N atoms in the same atomic layer are much stronger than the vdW forces between different atomic layers, the mechanical properties along D1 of these hot-pressed h-BN matrix composite ceramics are much lower than those along D2 and D3. The elastic moduli, flexural strength, and fracture toughness along D3 of these specimens with different YAG contents are all a little higher than those along D2, which is attributed to the fact that the crack initiation along D2 only needs to break the sp 2 covalent bonds combining B and N atoms, while the crack initiation along D3 needs to break both the sp 2 covalent bonds combining B and N atoms and the vdW forces between different atomic layers.
If we considered the orientation of h-BN grains as the only factor influencing the mechanical properties, the increase of the texture degree would make it more difficult for the crack to initiate and propagate along D2 and D3 as more sp 2 covalent bonds need to be broken, and make it easier for the crack to initiate and propagate along D1, resulting in higher mechanical properties along D2 and D3 and lower mechanical properties along D1. However, the YAG phase and relative density also have significant influence on the mechanical properties of these composite specimens. As YAG has a cubic crystal structure, it shows isotropic properties along different directions. Polycrystalline YAG ceramics were proven to show an elastic modulus around 300 GPa [35] and a flexural strength of 260-420 MPa [36], which are higher than the highest values of these h-BN matrix composite ceramics in this study. Therefore, the existence of YAG phase plays a role in strengthening and toughening along all loading directions of these composite specimens. Besides, a higher relative density, i.e., a lower porosity, also has a positive effect on the mechanical properties.  As the YAG content increases from 10 to 30 wt%, the texture degree keeps increasing, which plays a role in strengthening and toughening along D2 and D3 and is harmful to the mechanical properties along D1, while the increased YAG content and relative density are beneficial to the mechanical properties along all loading directions. As a result, the elastic modulus, flexural strength, and fracture toughness along D1 show no significant change, while those along D2 and D3 keep increasing. As for BN-4YAG, although its texture degree drops dramatically compared with BN-3YAG, its higher relative density and YAG content give it the highest elastic modulus, flexural strength, and fracture toughness along all loading directions.
For h-BN matrix composite ceramics with the absolute value of IOP lower than 600 and f 00l lower than 0.98 (BN-1YAG, BN-2YAG, and BN-4YAG), their work of fracture along D2 is lower than that along D3, which is the same as the fracture toughness. However, for BN-3YAG with an IOP value of −1778 and a f 00l value as high as 0.99, its work of fracture along D2 (155.7± 5.1 J/m 2 ) is 66.7% higher than that along D3 (94.5± 4.5 J/m 2 ) and also 41.6% higher than the highest work of fracture (110.0±8.1 J/m 2 ) of other specimens despite its lower fracture toughness along D2 than that along D3, which owes to its high texture degree and more tortuous crack propagation path along D2. Figure 8(e) shows the load-displacement curves along D1, D2, and D3 of BN-3YAG during the fracture toughness testing. The highest load along D3 is higher than that along D2, resulting in the higher fracture toughness along D3 than that along D2. BN-3YAG behaved a non-catastrophic fracture mode along D2 due to the orientation of platelike h-BN grains, which resulted in step-by-step dropping of load with increasing the displacement during the tortuous crack propagation. Such a non-catastrophic fracture mode was also observed in the load-displacement curves along D2 of other specimens (Fig. 8(f)), where BN-3YAG behaved more significant step-by-step dropping of the load with the increase of the displacement, which is consistent with its high work of fracture and could improve its reliability in actual use as directional thermal conductive materials. Therefore, it could be speculated that a texture degree high enough (f 00l > 0.98) could significantly improve the work of fracture with the loading direction parallel to the c-axis orientation (D2) of h-BN matrix composite ceramics. Figure 9 shows the lateral fracture morphologies along D2 of hot-pressed h-BN matrix composite ceramics with different YAG contents, where crack deflection phenomena correspond to the step-by-step dropping of load in Fig. 8(f). Different from BN-3YAG, the texture degrees of other specimens are not that high (f 00l < 0.98) to induce crack propagation paths tortuous enough which can give them higher work of fracture along D2 than that along D3. Figure 10 shows the linear thermal expansion curves and CTE measured parallel and perpendicular to the c-axis orientation of hot-pressed h-BN matrix composite ceramics with different YAG contents. The schematic diagram of samples with different orientations in the bulk specimen for testing anisotropy in thermal expansion is shown in Fig. 10(b). The CTE of YAG was calculated to be 7 × 10 −6 K −1 at room temperature based on the lattice parameters obtained from XRD [37]. The CTE perpendicular to the c-axis of h-BN is −2.72 × 10 −6 K −1 at room temperature and approaches zero at high temperatures, while the CTE parallel to the c-axis of h-BN is 37.7 × 10 −6 K −1 at room temperature and tends to be invariable at high temperatures [38].
A higher texture degree tends to decrease the CTE perpendicular to the c-axis orientation and increase that parallel to the c-axis orientation, while a higher YAG content tends to increase the CTE perpendicular to the c-axis orientation and decrease that parallel to the c-axis orientation. A higher relative density tends to decrease the CTE along both directions. As a result, the CTE perpendicular to the c-axis orientation keeps increasing with the increase of YAG content due to the contribution of YAG, while that parallel to the c-axis orientation first decreases from 10 wt% YAG to 20 wt% YAG attributed to the increased YAG content and relative density and then increases from 20 wt% YAG to 30 wt% YAG due to the increased texture degree.
With the further increase of the YAG content to 40 wt%, the CTE parallel to the c-axis orientation decreases drastically due to the increased YAG content and relative density and the decreased texture degree.
www.springer.com/journal/40145 Figure 11 gives the specific heat, thermal diffusivity, and thermal conductivity (25-1200 ℃) parallel and perpendicular to the c-axis orientation of hot-pressed h-BN matrix composite ceramics with different YAG contents and their anisotropy degrees in thermal conductivity, which were determined by the ratio of the thermal conductivity perpendicular to the c-axis orientation to that parallel to the c-axis orientation. The schematic diagram of samples with different orientations in the bulk specimen for testing anisotropy in thermal conductivity is shown in Fig. 11(b). The specific heat at room temperature of YAG was measured to be around 0.6 J/(K·g) by the differential scanning calorimeter [39], which was lower than 0.8 J/(K·g) of h-BN [40]. As a result, the specific heat of these composite specimens decreases with the increase of YAG content. Besides, all specimens have higher specific heat at higher temperatures.  The thermal diffusivity and thermal conductivity perpendicular to the c-axis orientation are higher than those parallel to the c-axis orientation of these composite specimens due to the anisotropic thermal conductivity of h-BN [8]. A higher texture degree would result in higher thermal conductivity perpendicular to the c-axis orientation and lower thermal conductivity parallel to the c-axis orientation. Although the thermal conductivity of YAG was measured to be only −10 W/(m·K) at room temperature [39], BN-3YAG, which has the highest texture degree, still has the highest thermal diffusivity (82.4 mm 2 /s at 25 ℃ ) and thermal conductivity (163.2 W/(m·K) at 25 ℃) perpendicular to the c-axis orientation, which are 11.1 times those (7.4 mm 2 /s and 14.7 W/(m·K), respectively) parallel to the c-axis orientation at room temperature. However, the higher relative density of BN-3YAG gives it higher thermal conductivity parallel to the c-axis orientation than that of BN-2YAG (11.7 W/ (m·K) at 25 ℃), and the higher YAG content of BN-2YAG (102.9 W/(m·K) at 25 ℃) gives it lower thermal conductivity perpendicular to the c-axis orientation than that of BN-1YAG (103.4 W/(m·K) at 25 ℃), although the former has a higher texture degree than the latter.
The thermal diffusivity and thermal conductivity of these composite ceramics decrease and the decrease rate gets slower with the increase of the testing temperature (Figs. 11(a) and 11(b)), especially for those parallel to the c-axis orientation, which is attributed to the influence of temperature on the mean free path of phonons. As the thermal conduction in h-BN is achieved by phonon transport, the thermal conductivity of h-BN matrix composite ceramics can be evaluated by Eq. (5) [41,42]: where k is the thermal conductivity achieved by phonon transport, C is the specific heat which increases with the increase of the temperature, v is the velocity of phonons which is influenced by the grain orientation, and l is the mean free path of phonons which is influenced by both the temperature and grain size. With increasing the testing temperature from room temperature, the number of phonons increases dramatically, which causes more scattering phenomena between phonons and phonons, leading to the decrease of l. Further increase of temperature diminishes l to be closer to the lattice spacing, and thus the adverse effect of temperature on the thermal conductivity decreases; meanwhile, the specific heat increases, making thermal conductivity tends to be constant. For the thermal conduction parallel to the c-axis orientation, the grain size (the thickness of plate-like h-BN grains) is very small, and thus there exist many phonon scattering phenomena caused by grain boundaries, that is, many phonons get scattered by grain boundaries before encountering other phonons. Therefore, the decrease rate of the thermal conductivity parallel to the c-axis orientation is slower than that perpendicular to the c-axis orientation, leading to the decrease of the anisotropy degree in the thermal conductivity with the increase the testing temperature, especially for BN-3YAG having the highest texture degree (Fig. 11(c)).

Conclusions
Textured h-BN matrix composite ceramics with different contents of YAG formed in situ from the 3Y 2 O 3 -5Al 2 O 3 sintering additive were fabricated by hot-pressing. The YAG phase shows good wettability to h-BN grains, and a coherent relationship between h-BN and YAG was observed with [010] h-BN //[111] YAG and (002) h-BN // (321) YAG . The rotation of plate-like h-BN grains under the uniaxial pressure in the liquid YAG phase during hot-pressing led to textured microstructures with the c-axis of h-BN grains preferentially oriented parallel to the sintering pressure direction, giving these composite ceramics anisotropy in their mechanical and thermal properties. They show better mechanical properties along D2 and D3 than those along D1, higher coefficients of thermal expansion parallel to the c-axis orientation than those perpendicular to the c-axis orientation, and higher thermal conductivity perpendicular to the c-axis orientation than that parallel to the c-axis orientation.
With increasing the YAG content, the texture degree of the composite specimen first increases, reaching the highest with a f 00l of 0.99 for the specimen with 30 wt% YAG. Further increase of the YAG content (40 wt%) resulted in the aggregation of liquid YAG phase which caused the buckling of h-BN plates and significantly reduced the texture degree. The anisotropy degrees in mechanical and thermal properties of these composite specimens are not totally consistent with their texture degrees due to their different YAG contents and relative density.
Textured h-BN matrix composite ceramics behave a non-catastrophic fracture mode when the crack propagates through h-BN layers, among which the specimen with 30 wt% YAG has the highest work of fracture (155.7±5.1 J/m 2 ) due to its highest texture degree, which also has the highest anisotropy degree in thermal conductivity (14.7 and 163.2 W/(m·K) parallel and perpendicular to the c-axis orientation at 25 ℃, respectively), making it a promising candidate as thermal management devices in electronics.