abstract: While it is understood by women that tech fluency might act as a powerful career accelerator or even a new career direction towards software engineering, this awakening often comes after graduation from a different field, when it is difficult for the women to make the shift towards tech and computing. In this paper, we report on our experience with running a successful education non-profit called Czechitas, which shows that women in their 20s and 30s are (maybe surprisingly) highly interested in learning tech, they just need a helping hand and tailored assistance, encouragement and guidance.
author:
title: Women Want to Learn Tech:
  Lessons from the Czechitas Education Project

Shell *et al.*: Bare Demo of IEEEtran.cls for IEEE Journals

# Introduction

Across the world, as few as 10–25% of tech professionals are women, with many women dropping their interest in tech for preventable reasons. This is a major issue for both the current and future jobs market, as Europe alone faces a shortage of hundreds of thousands skilled tech engineers to date . Moreover, the concerns are related also to the impact a tech background can have on a women's career, and the economic potential that accompanies it .

Many studies have been conducted to date to reveal that in early education, girls and boys are equally touched by technology, computing and software engineering. But as adolescents (age 13 to 18), girls become much less interested in these topics, with limited recovery during the next 10 years that are crucial in terms of career choice .

While these *vulnerable years* and all factors that influence whether girls keep their interest in tech or not are very well studied and understood , barely any attention is paid to the period in life of young women some years after graduation when they reconsider their career choice (we call it *reconsideration years*, typically between 25 to 30 years of age). In our experience, this time might be linked to the moment when some of these women plan starting a family and consciously rethink what shall be the career they want to come back to after maternity leave.

Experience shows that during these *reconsideration years*, the interest in tech and software engineering among young woman raises again due to prospective opportunities that tech offers . If however a helping hand is not in reach to the young women, the interest goes in vain.

In this paper, we report on our experience with helping young women during these reconsideration years to find their way into tech, encourage and educate them. All that within the Czechitas non-profit education project, started specifically for this purpose in the Czech Republic, where the ratio of women in tech is the second worst in EU (see Figure ). Within Czechitas, we assist these young women in their learning path, we guide them towards specific technologies, and advise them in their career decisions. We connect them with professionals in the respective fields, and create welcoming community environment for them to learn from each other. Besides sharing this experience, this paper explains what kind of helping hand is needed during the reconsideration years so that the interest of these women in tech turns effective.

# The Story of Czechitas

Czechitas is a Czech non-profit organization that has emerged in 2014 from a simple idea to bring tech closer to girls, and girls closer to tech. Over the time, this idea has attracted a strong community of tech professionals, companies and volunteers, and gave rise to a portfolio of female-tailored courses in various areas of tech and software engineering, such as programming, web development, mobile app development, data science, testing, digital marketing and graphic design.

Thanks to the success of our education activities, consisting of hundreds of events a year (each receiving multiple times more registrations than its capacity), we have become recognized as the leading platform in the Czech Republic actively addressing gender diversity in STEM. We have also been recognized abroad with numerous European as well as overseas awards, e.g. becoming the first organization from the Central Europe awarded the prestigious Google.org grant in 2015, and in 2018 becoming the first organization worldwide being awarded the Google.org grant for the second time .

We have influenced over 10,000 women who graduated from our courses to either choose or change their career path to tech or use their new tech skills to advance their careers. Besides students, our community includes lecturers, professionals, partner companies, volunteers, the core team and social media audience (community of 15,000+). Even though the core team includes only 30 employees, we have the support of more than 350 volunteers actively involved in our activities (tech professionals, 40% female).

# Key Factors when Reconsidering Tech

Tech by itself is an appealing career choice, with high salary, job security, working flexibility, space for creativity and career growth. What is interesting then are the factors that negatively influence the decision to pursue tech career among women.

During the adolescence of girls, social encouragement and tech exposure in school, together with self perception and career perception are reported to be the key controllable indicators for whether or not girls decide to pursue a Computer Science or tech path . While numerous other studies exist that examine the *vulnerable years* and confirm these findings , little is known about the factors that play the key role during the *reconsidearation years*, hindering the decision in favour of tech.

In 2017 we have surveyed 302 participants of Czechitas events together with representatives of 36 companies to understand these factors . Our findings (see Figure ) very well match the conclusions of Harvey Mudd's President, Maria Klawe, who compiled their experience in three points, saying "Number one is they think it's not interesting. Number two, they think they wouldn't be good at it. Number three, they think they will be working with a number of people that they just wouldn't feel comfortable or happy working alongside ." Let us elaborate on these three points and add our observations, resulting also in a fourth point, which is that young women during reconsideration years critically miss guidance in learning.

## They think it's not interesting

Alarming effect of low exposure of girls to tech was described by Google in 2014, showing that the unfamiliarity with tech alone causes inclination to negative connotations with tech as a domain . The study found that no matter if the girls were exposed to tech via compulsory or elective classes, the exposure itself was enough for them to characterize tech in much more positive terms. While girls unfamiliar with tech were using words like "boring, hard, difficult", girls previously exposed to tech used words like "future, fun, interesting, exciting" (see Figure ).

One of the first corrective tasks we need to address is then to turn around the perception of tech that young women have due to their low exposure to tech in the past. Without that, young women have difficulty to perceive tech as a career that fulfills both the creative passion (inventing, problem solving, exploration, etc.) and the intangible, social passions (helping people, conservation, medical breakthroughs, etc.) that make a profession personally rewarding .

## They think they wouldn't be good at it

The key factor in our study among women participants of Czechitas courses is the perception of their own proficiency in tech (see Figure ) . Interestingly, the perception is rooted more strongly in their confidence than objective reasons, which is commonly observed as the Confidence Gap between men and women reported in various fields .

Similar conclusions have been drawn in mathematics, where although standardized tests indicate that girls score just as well as boys in math, and that in school grades, girls often even outperform boys , girls tend to have higher levels of math anxiety and lower levels of confidence in their math skills. In other words, even when girls show similar performance levels to boys, they are often less sure of themselves, which applies to math as well as tech .

## Who they will be working with

In our experience, women we teach are highly concerned about the work environment they shall once belong to, being afraid that they might end up working with people that they would not feel comfortable working alongside . And indeed, some studies report that 27% of women cited discomfort with their work environment as a factor why they left their IT job. Furthermore, workplace policies not suited to women also play a role . Our study moreover shows that women perceive work in tech jobs as absent of communication among people (only 1% considered soft skills useful in tech career in our study ), which they do not feel comfortable with.

## Missing guidance

In one form or another, young women during their reconsideration years feel lost when navigating tech education on their own. Helping hand or guidance during that time is at the end what decides whether they engage in tech education or not. Sometimes, they only need a direction, such as somebody who helps them to decide what programming language to start with. There is no universal answer to this question, which is why they often feel confused when they try to find the answer online. Moreover, they often ask how much time they shall be ready to invest before seeing results (e.g. in terms of an application they would like to develop), what learning path takes them to a specific job interview, which of the thousands of online courses they shall start studying. Then during learning, they often get stuck on trivial mistakes, which they find impossible to debug without help of another person. And they miss community of like minded female students.

# Our Solution

Research by Accenture and Girls Who Code concludes that improving universal access of women and girls to tech and computing will not by itself address the gender gap. They instead suggest that only by tailoring courses to girls' specific needs can we boost their commitment to computing . In this section, we discuss what such tailoring means to us.

## Community

We build a community of girls and women who learn to code, as we see that community feeling is one of the key factors in retaining these tech students. All our events have very informal and friendly atmosphere, with team building activities and food sharing, to enforce the community feeling.

## Safe environment and encouragement

Besides a lecturer, each of our coding events (for 30 participants) features around 5 more teacher mentors (we call them couches) who provide individualized help, support and encouragement to the students. Couches often sit by the students, helping them with whatever difficulty they have, which creates friendly atmosphere in which the students feel comfortable asking questions and reporting troubles they are facing (i.e. safe environment for trial and error). It moreover facilitates varying pace of each participant, thanks to which everybody can study on the speed that best fits their needs. This is in line with a study that shows that over a half of the girls learning tech would like to receive more active encouragement from teachers .

## Practical experience and hands-on exercises

All our events feature many fun and creative hands-on exercises. Studies show that the more practical experience a girl receives during her education––inside or outside the classroom––the higher interest in STEM she builds . Creativity in the classroom is also key. Girls who like STEM enthuse about being able to choose their own projects or go on field trips where STEM subjects are brought to life.

## Real-life applications

Our experience supported by studies also shows that women become more interested in tech once they are able to conceive how tech can be applied to real-life situations and how relevant it might be to their future . This is why we integrate real-life applications and projects in our activities.

## Role models

Having visible female role models sparks women interest in tech careers and helps them to picture themselves pursuing these fields. This is why we organize meet-ups with inspirational female engineers and organize Czechitas Thesis Award, which showcases examples of excellent graduation theses in tech from Czech universities, drawing attention to talented girls who are passionate about tech.

## Career guidance

We offer women a helping hand on their career path, from the initial decision about the field of study to the job interview. We organize annual Job Fair, which is a vibrant event where our community of students, mentors, lecturers and volunteers can meet our partner companies. The vision of this event is to get these two groups together for one day and provide them with specialized content, including career and IT workshops for our community and gender or educational workshops for our partner companies.

## Partnering with companies

It has been recognized that young women are more likely to pursue tech careers when they are confident that men and women will be treated equally working in these disciplines . We have observed that women feel much more confident when filing their job application to a Czechitas partner company, because the fact that the company partners with us gives them the reassurance that they are welcome at job interview. We also educate our partner companies about gender bias within the hiring process and what makes them more/less attractive employer to women.

# Format of Events

## One-day workshops

The most popular format of our events are full-day workshops, taking place typically on Saturdays or Sundays. These focus on introduction to programming, web and mobile apps development, embedded systems, graphic design, data science and digital marketing. Most popular programming languages are Javascript, Python, Java, C#, PHP, Kotlin, and Swift. The topics also include various software engineering essentials, such as testing, agile development, project management, security engineering, or version control.

The aim of one-day workshops is to spark interest in the attendees about coding and other areas of software engineering, and to give them necessary encouragement and guidance to study further.

## Long-term courses

We yearly organize over 20 spring and autumn courses, focusing mainly on various programming languages and technologies (most popular are Java, C#, Javascript, PHP, HTML/CSS, Android, UX design, testing). Each of these courses comprises of 10-14 evening sessions of 2-3 hours each. The aim of long-term courses is to teach women a specific programming language or tech skill, and equip them with confidence for starting their own software development.

## Re-qualification academies

For the women who made the decision to change the career, we offer an intense re-qualification program, which involves three months of lectures, hands-on exercises, and mentoring, called Digital Academy. Throughout the Academy, students participate in internships with companies and work on their final graduation projects. We offer them career services focused on professional transit to tech.

We have developed this unique type of re-qualification program (with successful 70% placement of the graduates on junior job positions in partner companies) for two fields: (1) data science and (2) backend web development in Java. The *Digital Academy: Data* includes teaching blocks on data mining, data processing and visualization tools, statistics, digital marketing, SQL and databases, data science and discovery, and Python as the programming language. The *Digital Academy: Java* focuses on key skills and technologies needed to develop web application with the focus on its backend in Java.

## Hackathons

To promote creativity and confidence among participants of our courses, we organize hackathons, where the women can test their knowledge during intense work on a full frame team project under the supervision of experienced mentors. The hackathons are organized both as stand-alone events and as part of the academies and other intense education programs.

## Coding clubs and Internships

Coding clubs and internships give girls the opportunity to work on a long-term team coding project under the supervision of an experienced professional. In case of the coding club, the team meets 2-4 times per month to discuss their progress on the joint project. In case of the internship, the woman becomes part of a team in a company that provides support for her to learn in practice.

## Czechtias New Generation

Over the years, we have become recognized as a platform with the ability to guide newcomers to tech. Thanks to that, we have been invited to engage in tech education of children and youth in the Czech Republic. Since then, we have not only been active in advancing tech education among kids (engaging female role models in teaching), but also in equipping their school teachers with the necessary skills to make tech education accessible and fun.

# Conclusion

This paper reports on our successful education project, called Czechitas, which supports women when restarting their career and changing direction towards tech and computing. The project enjoys large popularity, and hence is now growing into multiple new cities within the Czech Republic (from 5 cities where we have our teams now, growing to 3 more cities in the next two years), while enlarging its portfolio with new topics and forms of learning (building an online learning platform now), as well as support for educators (training of tech teaching skills, organizational platform for knowledge sharing and methodical support). Moreover, we are now putting effort in improving the way Czech schools approach tech education, to prevent the low female engagement in tech, which we are correcting now.


abstract: Cause of a rapid increase in technological development, increasing identity theft, consumer fraud, the threat to personal data is also increasing every day. Methods developed earlier to ensure personal the information from the thefts was not effective and safe. Biometrics were introduced when it was needed technology for more efficient security of personal information. Old-fashioned traditional approaches like Personal identification number( PIN), passwords, keys, login ID can be forgotten, stolen or lost. In biometric authentication system, user may not remember any passwords or carry any keys. As people they recognize each other by the physical appearance and behavioral characteristics that biometric systems use physical characteristics, such as fingerprints, facial recognition, voice recognition, in order to distinguish between the actual user and scammer. In order to increase safety in 2005, biometric identification methods were developed government and business sectors, but today it has reached almost all private sectors as Banking, Finance, home security and protection, healthcare, business security and security etc. Since biometric samples and templates of a biometric system having one biometric character to detect and the user can be replaced and duplicated, the new idea of merging multiple biometric identification technologies has so-called multimodal biometric recognition systems have been introduced that use two or more biometric data characteristics of the individual that can be identified as a real user or not.
author:
title: Analyze and Development System with Multiple Biometric Identification

# Introduction

In today's, security is the main problem and everybody will be looking for reliable ways to ensure safety. Old ways of securing our homes and things used key and lock systems, which is quite inefficient keys can be duplicated. Later, the next level of protection was provided by modern technologies with PIN and password. Modern systems are quite effective, but they can also be hacked and gain access. To achieve with high accuracy and reliability, a biometric identification technology has been introduced that uses physical or behavioral characteristics to verify the individual. The main reason for choosing biometric identification The technology above the traditional identification method is that traditional identification techniques are vulnerable security, unreliable and uncertain. Most developed biometric systems use one biometric technology for biomass identification and verification. These unimodal systems are unable to meet some of the challenges such as enrollment large population, high performance etc. Later advanced biometric identification technology with better Performance and a high degree of accuracy, called a multimodal biometric recognition system, are introduced. Multimodal systems are capable of using two or more physical or behavioral characteristics of a person write, verify and provide a decision. Multimodal biometric systems also combine two or more identification technologies, such as systems, fingerprint systems and face recognition, with more sensors and provides decisions by combining results from each subsystem. When using multimodal biometric systems, if one sensor or the system fails to carry out identification, other systems can be used to make decisions. Multimodal systems can be used in various fields, for example in homes for the biometric door locks to avoid theft, hospitals to record details of patient in software companies and educational institutions, to ensure participation in the banks for safe and secure payment processing, biometric alarm systems in case of data security in different sectors. Multimodal systems are more convenient and overcome all limitations of unimodal systems with increased security, high accuracy and liveliness.

# General Statement of The Multimodal Authentication Problem

The task of multimodal authentication consists of several steps:. Get results from individual systems.. To normalize these results.. Perform the merger of normalized components.

# Related works

In recent years, a lot of work has been done on different types of biometric identification techniques and their work, advantages and functioning of multimodal biometrics compared to unimodal and traditional approaches. Biometrics refers to the automatic authentication of a person based on his physiological and behavioral property. The use of biometrics as a reliable method of verification occurs in almost all EU countries government and private sector. Some limitations of unimodal biometric systems may be reduced by the use of multimodal biometric systems that integrate information at different levels to improve performance. In this the paper proposed a new multimodal biometric recognition system with a fusion of fingerprints and irises recognition to achieve higher accuracy. The main goal of this work is to show that the merger of unrelated, independent identification techniques achieve better accuracy than any unimodal biometric systems. Output scores from two different systems merge into a single score at the score level using three standardizations methods and four Fusion approaches. They also used standardization steps to transform unrelated, different scores to the common scale of the two techniques. In this article, they told about the types of mergers level in biometric systems. Sensor level, function level, score level are types of Fusion levels. First step on the score level is the normalization of the score, when the scores are transformed into a common format. Step two the Fusion level has two categories, such as classification and combination. Classification approach classifies decisions to a cheater or real. Combined approach combines score multiple decisions into one comparative score. The accuracy achieved by this work is more than 60 percentage. By this article, I understand that a merger at the rank level can be a fairly new merger access. Fusion is carried out at the level of rank, when the output of each biometric system was a subset of possible matches sorted by declining self-esteem. The intention of the merger at rank level was to combine the output individual biometric subsystems to achieve greater accuracy. They used facial recognition, ear signature recognition and recognition systems and processed ones that use Eigen and Fisher face projection projection for registration and identification of the face, projection of own ears for registration and identification of ears, and Self-Signature projection for registration and signature identification. So use the assignment function based on in the ranks of the face, ear and signature there was a merger at the level of rank, and the final exit from them was obtained. From their by experiment, they came to the conclusion that the fishing face works more efficiently than the technique of its own face and of the three biometric face recognition systems, it worked or worked better than its own ear and your own signature. In this article, they considered the extraction and evaluation of fingerprints, as well as the extraction of the veins of the fingers and conversion to the knowledge of a person. They also combined two biometric recognition systems using Enhanced partial discrete Fourier transformation (EP-DFT), since the DfT Matrix has its own characteristics, such as orthogonality. Like The DfT Matrix shortens the time for its multiplication by The Matrix, so that it becomes a partial discrete Fourier transformation compact and computationally efficient. They considered three different merger options where they have considered a set of fingerprint functions based on markants and a set of fingerprint functions based on an image using an element level Fusion Strategy. Their analysis concluded that EP-DFT significantly increases system security compared to EP-DFT original P-DFT. They also concluded that it is necessary to design a good, effective irreversible transformation functions for cancellable multi-biometric systems.

## EIGENFACES METHOD

The basis of the eigenfaces method is the Principal Component Analysis (PCA). Eigenfaces and PCA have been used by Sirovich and Kirby to represent the face images efficiently. They have started with a group of original face images, and calculated the best vector system for image compression. Then Turk and Pentland applied the Eigenfaces to face recognition problem. The Principal Component Analysis is a method of projection to a subspace and is broadly used in pattern recognition. An objective of PCA is the replacement of correlated vectors of large dimensions with the uncorrelated vectors of smaller dimensions. Another objective is to calculate a basis for the data set. Main advantages of the PCA are its low sensitivity to noise, the reduction of the requirements of the memory and the capacity, and the increase in the efficiency due to the operation in a space of smaller dimensions.

## VIOLA-JONES

This method is highly effective for search objects in images and video sequences in real time. This detector has an extremely low probability of false detection of the face. The method works well and detects facial features even when observing an object at a slight angle, up to about 30°. Recognition accuracy using this method can reach values over 90 percentages, which is a very good result. When the angle of inclination is more than 30°, the probability of detecting a face drops sharply. The indicated feature of the method does not allow in a standard implementation to detect a person's face, turned at an arbitrary angle, which greatly complicates or makes it impossible to use the algorithm in modern production systems, taking into account their growing needs.

# PROPOSED APPROACH

Individual or unimodal biometric systems are not enough for a person to effectively recognize a person, so I suggest multi biometric system for authentication to achieve better safety and reliability. I focus in particular on the use and degree of accuracy achieved by multibiometric recognition systems in domestic security and various area. I will explore more about biometric recognition systems in different areas.

## METHODOLOGY

The following steps are procedures for identifying an individual.

1.  **Enrollment of the biometric samples of a candidate** Enrollment is when a candidate first uses a biometric system, sample a biometric character as a fingerprint or face recognition is captured and stored in a database for later comparison.

2.  **Obtain/Acquire Live Samples from Candidate** First, get a biometric sample from the user using sensors.

3.  **Extract feature** Process and extract prominent features from the captured sample.

4.  **Comparison of samples** Process and extract prominent features from the captured sample....

5.  **Display decision** When the input samples coincide with the registered samples, it accepts the candidate or refuses.

# FLOW DIAGRAM

The below is the flow of working procedure of multi biometric recognition systems.

-   Initially, the sample biometric function is captured from the user when the user wants access. Biometric features can be fingerprints, voice recognition, face recognition.

-   The captured biometric data is pre-processed, which includes the removal of unwanted data, noise and a prominent feature is obtained on which a person is to be recognized

-   The extracted function is compared to registered samples that are stored in the database both samples have similarities.

-   The logged sample and the captured sample are compared at different angles using some defined algorithms and corresponding accuracy are displayed.

-   Based on the degree of accuracy, the decision is made as a real user or fraudster

-   If one biometric system provides 60 percent accuracy, which is not sufficient for any decision other a system with a different biometric character is used in the same procedure.

-   The final results obtained from all subsystems are scaled and transformed into a normal format; and the decision appears on the screen.

Nowadays, biometrics have become a new trend in public and private workspace, where data or any other personal security is a major concern. Most nations move away from the traditional way to secure houses of technology road. Instead of the lock and key system, more biometric door locks are used, which is a combined system of two or more biometric identification techniques for verification and verification. The biometric door will have a latch with fingerprint recognition or palm print recognition in order to identify the homeowner. And were introduced videos, which offer an easy way to find out who's at the door, without get close to the door. Bells with video can be connected to smartphones via Wi-Fi and the user can get alert when someone approaches the door. Video Bell records video when the bell is pressed and it offers two-way audio communication that allows you to communicate with the visitor from anywhere through your phone. Other features include facial recognition technology, which identifies visitors by name, motion sensing technology that knows the difference between humans, cars and animals.  Biometrics is also used in many organizations and educational institutions to record employees participation. The Multi biometric attendance system provides employees with flexibility in recording their presence; and working time at the right time interval. The main advantage of using multimodal biometric systems is fast identification and authentication. If one of the modalities fails to determine and submit a decision, the user can choose different modality or based on the convenience of the user, he can choose a modality to gain access. Biometrics are it is also used in financial sectors such as banks, ATMs, insurance companies for Secured Transactions. Financial sectors are in terms of transactions very sensitive, because spoofing can lose a huge amount of money. To ensure secured transactions, financial sectors adopt biometric identification techniques identify account holders and also for KYC. Scanning the veins and scanning the iris are popular today where the user does it it does not require contact with the sensor and is more reliable and faster than other biometric technologies. Healthcare the Centers also use biometrics to identify the right patient, to provide a medical history and patient treatment plan.

# Conclusion

In this digital era, most existing biometric systems use information from a single biological function verification and identification. The captured sample is due to the vulnerability of the biometric sensor to noisy or bad data biometric functions could become fake and it will be possible to reject the Registered User incorrectly and falsely accept the fake user. In the case of a fingerprint recognition system, scanners may not recognize dirty or injured fingers also underdeveloped fingerprint combinations in the case of young children and faded fingerprint combinations in the case of the elderly. Where, as in the case of a facial recognition system, sensors may not be able to detect identical twins, whose appearance will be almost similar. Also, in the case of voice recognition systems, a person can imitate the voice of another person and gain access to their personal information. In all of the above cases due to an erroneous decision from the scanning devices gives unauthorized access and leads to fraudulent activity. In order to overcome these problems, a biometric system with a merger of two or more biological characteristics and sensors called multi-biometric systems are used. The main advantage of using multimodal biometric recognition there is an increased accuracy of verification. Thus, the use of multi-biometric recognition systems helps us in identifying and more precisely, in comparison with the use of a single biometric recognition system. I learned about different biometric recognition systems in order to use the necessary multimodal method biometrics for smart homes to have safe and reliable access to our homes. I can do this multimodal biometrics are more reliable and accurate compared to unimodal biometric recognition systems.


author: Eduardo Avila; Márcio Dorn; Clarice Sampaio Alho; Alessandro Kahmann
bibliography: main.bib
title: Hemogram Data as a Tool for Decision-making in COVID-19 Management: Applications to Resource Scarcity Scenarios

# Introduction

Since its first detection and description , COVID-19 expansion has brought worldwide concerns to governmental agents, public and private institutions, and health care specialists. Declared as a pandemic, this disease has deeply impacted many aspects of life in affected communities. Relative lack of knowledge about the disease particularities has led to significant efforts devoted to alleviating its effects .

Alternatives to mitigate the disease spread include social distancing . Such a course of action has shown some success in limiting contagion rates . However, isolation policies manifest drawbacks as economic impact, with significant effects on macroeconomic indicators and unemployment rates . To address this, governments worldwide have proposed guidelines to manage the essential workforce, considered pivotal for maintaining strategic services and provide an appropriate response to the pandemics expansion .

Widespread reports of threats to critical national infrastructure have been presented, with significant impact associated with medical attention . Significant pressure is being faced by emergency response workers, with some countries on the brink of collapse of their national health systems . The main concern associated with COVID-19 is the lack of extensive testing capacity. Shortage of diagnostic material and other medical supplies pose as a major restraining factor in pandemics control .

The most common COVID-19 symptoms are similar to other viral infectious diseases, making the prompt clinical diagnostic impractical . Official guidelines emphasize the use of quantitative real-time PCR (qRT-PCR) assays for detection of viral RNA in diagnosis as the primary reference standard . In Brazil, test results are hardly available within at least a week, forcing physicians and health care providers to take strategic decisions regarding patient care without quality information.

Previous reports have described alterations in laboratory findings in COVID-19 patients. Hematological effects include leukopenia, lymphocytopenia and thrombocytopenia, while biochemical results show variation on alanine and aspartate aminotransferases, creatine kinase and D-dimer levels, among other parameters . Some efforts have been applied to evaluate clinical and epidemiological aspects of this disease using computational methods, such as diagnosis, prognosis, symptoms severity, mortality, and response to different treatments. A useful review of some of these methods is presented by Wynants and collaborators .

The main objective of this article is to provide insights to healthcare decision-makers facing scarcity situations, as a shortage of test capacity or limitations in the essential workforce. A useful method of doing so is using hemogram test results. This clinical exam is widely available, inexpensive, and fast, applying automation to maximize throughput. To do so, we have analyzed hemogram data from Brazilian symptomatic patients with available test results for COVID-19. We propose a framework using a *Naïve-Bayes* model for machine learning, where test conditions can be adjusted to respond to actual lack of resources problems. Finally, four distinct scarcity scenarios examples are presented, including handling of the essential workforce and shortage of testing and treatment resources.

# Material and Methods

## Data Collection

5644 patients admitted to the emergency department of *Hospital Israelita Albert Einstein* (HIAE - São Paulo, Brazil) presenting COVID-19-like symptoms were tested via qRT-PCR. A total number of 599 patients (10.61%) presented positive results for COVID-19. The full dataset contains patients anonymized ID, age, qRT-PCR results, data on clinical evolution, and a total of 105 clinical tests. Not all data was available for all patients, therefore the number of missing information is significant. All variables were normalized to maintain anonymity and remove scale effects. No missing data imputation was performed during model generation to avoid *bias*. Considering the significant ammount of missing data, only 510 patients presented values for all 15 parameters evaluated in hemogram results (comprising the following cell counts or hematological measures: hematocrit, hemoglobin, platelets, mean platelet volume, red blood cells, lymphocytes, leukocytes, basophils, eosinophils, monocytes, neutrophils, mean corpuscular volume (MCV), mean corpuscular hemoglobin (MCH), mean corpuscular hemoglobin concentration (MCHC), and red blood cell distribution width (RDW). Data for the above parameters were used in model construction, along with qRT-PCR COVID-19 test results. The full dataset is available in <https://www.kaggle.com/einsteindata4u/covid19>.

## Machine Learning Analysis - *Naïve Bayes* Classifier

Machine learning (ML) is a field of study in computer science and statistics dedicated to the execution of computational tasks through algorithms that do not require explicit instructions but instead rely on learning patterns from data samples to automate inferences . These algorithms can infer input-output relationships without explicitly assuming a pre-determined model . There are two learning paradigms: supervised and unsupervised. Supervised learning is a process in which the predictive models are constructed through a set of observations, each of those associated with a known outcome (*label*). In opposition, in unsupervised learning, one does not have access to the labels, it can be viewed as the task of "spontaneously" finding patterns and structures in the input data.

Our objective with this study is to predict in advance the results of the qRT-PCR test with machine learning models using data from hemogram tests performed on symptomatic patients. The main process can be divided into four steps: (1) *pre-processing of the data* (2) *selection of an appropriate classification algorithm*, (3) *model development and validation*, i.e., the process of using the selected characteristics to separate the two groups of subjects (positive for COVID-19 vs. negative for COVID-19 in qRT-PCR test), and (4) test generated model with additional data. Steps are detailed as follows:

*Data Pre-processing:* Samples presenting a missing value in any of the 15 evaluated features were removed. A total of 510 patients (73 positives for COVID-19 and 437 negatives) presented complete data and were considered for the model construction.

*Classification Algorithm:* In this work, we use the *Naïve Bayes* (NB) classifier, which is a probabilistic machine learning model used for classification tasks. The main reasons for choosing this classifier are due to their low computational cost and clear interpretation. In medicine, the first computer-learn attempts in decision support were based mainly on the Bayes theorem, in order to aggregate data information to physicians' previous knowledge . The *Naïve Bayes* (NB) method combines the previous probability of an event (also called **prior* probability*, or simply **prior**) with additional evidence (as, for example, a set of clinical data from a patient) to calculate a combined, conditional probability that includes the **prior** probability given the extra information. The result is the *posterior probability* of an outcome, or simply *posterior*. This classifier is called "naïve" because it considers that each exam result (variables) is independent of each other. Once this situation is not realistic in medicine, the model should not be interpreted . Besides this drawback, it can outperform more robust alternatives in classification tasks, and once it reflects the uncertainty involved in the diagnosis, Bayesian approaches are more suitable than deterministic techniques .

*Model Development and Validation:* A classifier is an estimator with a predict method that takes an input array (test) and makes predictions for each sample in it. In supervised learning estimators (our case), this method returns the predicted labels or values computed from the estimated model (positive or negative for COVID-19). *Cross-validation* is a model evaluation method that allows one to evaluate an estimator on a given dataset reliably. It consists of iteratively fitting the estimator on a fraction of the data, called training set, and testing it on the left-out unseen data, called test set. Several strategies exist to partition the data. In this work, we used the *Leave-one-out* (LOO) cross-validation model, as in Chang et al. . The number of data points was split N times (number samples). The method was trained on all the data except for one point, and a prediction was made for that point. The proposed approach was implemented in *Python v.3* (<https://www.python.org>) code using *Scikit-Learn v. 0.22.2* as a backend.

*Model Test:* In order to evaluate the adequacy and generalization power of the proposed model, a set of 92 samples (10 positives for COVID-19 and 82 negatives) was extracted from the patient database. Those samples were not initially employed in model delineation, considering they present a single missing value among all 15 employed hemogram parameters. Missing data for this training set was imputed using the average value of the missing parameter within the resulting group (positive or negative). The test set was submitted to the previously generated model in order to evaluate classification performance.

# Results

## Descriptive Analysis

For data description, probability density function (PDF) of all 15 hemogram parameters were estimated through the original sample by kernel density estimator. Some hemogram parameters present notable differences between the distributions of positive and negative results, mainly regarding its modal value (distribution peak value) and variance (distribution width). Differences are summarized in Table . Regarding basophiles, eosinophils, leukocytes and platelets counts, qRT-PCR positive group distribution shows lower modal value and lower variance. On the other hand, monocyte count displays opposite behavior, once lower modal value and variance are observed for the qRT-PCR positive group. Lower variance may depict a condition pattern, therefore it is expected that negative cases present higher variance once it may contain a higher variety of conditions (reasons for symptom presence). The remaining nine hemogram parameters did not show a notable difference between negative and positive groups. PDF analysis results are presented in Supplementary Material Figure S1.

| **Parameter** | **Modal value** | **Variance** |
|:---|:--:|:--:|
| Basophiles | Reduced in positive cases | Reduced in positive cases |
| Eosinophiles | Reduced in positive cases | Reduced in positive cases |
| Leukocytes | Reduced in positive cases | Reduced in positive cases |
| Monocytes | Augmented in positive cases | Augmented in positive cases |
| Platelets | Reduced in positive cases | Reduced in positive cases |
| Parameters not shown displayed no difference between groups |  |  |

Descriptive analysis of hemogram parameters used in present study.

## *Naïve Bayes* Model Results

A NB classifier based on training set hemogram data was developed. Under the model, the complete range of *prior* probabilities (from 0.0001 to 0.9999 by 0.0001 increments) was scrutinized, and *posterior* probability of each class was computed for different *prior* conditions. A *posterior* probability value of 0.5 was defined as the classification threshold in one of the positive or negative predicted groups. Resulting model showed a good predictive power of the qRT-PCR test result based on hemogram data. Figure shows the accuracy, sensitivity, and specificity curves derived from the model for different *prior* probabilities of each class (positive or negative for COVID-19). Reported *prior* probabilities refer to positive COVID-19 condition.

When setting the **prior** probability to the maximum defined value (0.9999), the NB classifier correctly diagnosed all PCR positive cases. On the other hand, such configuration improperly predicted 77.3% of negative PCR results as positive. Regarding the lower possible **prior** probability setting, it does not classify a single observation as positive. This result can be explained by the unbalanced number of observations for each class, tending to over classify samples as the class with more observation, i.e. negative results. Such characteristics can also be noticed in the general accuracy, since the decrease in the *prior* ponce the classifier tends to diagnose all observations as belonging to the dominant class (negative) and consequently raising the total of correctly classified samples. The break-even point is met when *prior* probability is set to 0.2933. Under this condition, all metrics are approximately 76.6%.

Regarding the model sensitivity, the rate of positive samples correctly classified is over 85% within 0.999 to 0.5276 range, with small decrease of it when the *prior* probability of positive result is diminished within this range. When *prior* is set to under 0.0606, the number of positive predicted samples decrease rapidly, yielding lower sensitivity. As for specificity, it presents linear growth as tested *priors* decrease. Ultimately, the accuracy results profile are similar to specificity, due to the negative patients dominance.

As mentioned above, *prior* probability choice has a critical relevance in proposed model use. It is clear that, when extreme values of positive probability are applied (close to 0 or 1), specific classes (positive or negative qRT-PCR test results predictions) are favoured, increasing its ability of correct detection. As an example, when a value of 0.9999 is set for **prior** probability of positive result is set, an increase in misclassification in negative class results is observed. At the same time, it is possible to properly identify samples where hemogram evidence strongly indicates a negative result, according to the model. This is based on the fact that evidence used in the model construction (in present case, hemogram data) must strongly support the reduction of *posterior* probability of disease to values under 0.5, therefore leading to a negative result. This logic can be applied to fine tune the *prior* probability used in the model, in order to improve correct classification of positive or negative groups prediction. Examples of how to use this feature is provided in the "Discussion" section. Test samples (n=92, including 10 qRT-PCR positives) were used to test the proposed model. Figure presents results obtained from the model application to test dataset.

# Discussion

Laboratory findings can provide vital information for pandemics surveillance and management . Hemogram data have been previously proposed as useful parameters in diagnosis and management of viral pandemics . In the present work, an analysis concerning hemogram data from symptomatic patients suspected of COVID-19 infection was executed. A machine learning model based on *Naïve Bayes* method is proposed in order to predict actual qRT-PCR from such patients. The presented model can be applied to different situations, aiming to assist medical practitioners and management staff in key decisions regarding this pandemic. Figure summarizes model construction and application. Predictions are not intended to be used as a diagnostic method since this technique was designed to anticipate qRT-PCR results only. As such, it is highly dependable on factors affecting qRT-PCR efficiency, and its prediction capacity is dependent on the sensitivity, accuracy, and specificity of the original laboratory exam .

Descriptive analysis of hemogram clinical findings shows differences in blood cell counts and other hematological parameters among COVID-19 positive and negative patient results. Differences are conspicuous among three measures (leukocytes, monocytes and platelets) and more discrete to additional two (basophiles and eosinophiles). It is possible that differences are also present across the complete data spectrum, even though they are not clearly visualized with PDF data. These results are in accordance to previous reports of changes in laboratory findings in COVID-19 infected patients, where conditions as leukopenia, lymphocytopenia and thrombopenia were reported . It is important to highlight that data analysis is not sufficient to characterize clinical hematological alterations in evaluated patients (when compared to demographic hematologic parameters data), once data was normalized for the evaluated sample set only. However, even within this particular quota of population (individuals presenting COVID-19-like symptoms), differences were found between individuals presenting negative or positive qRT-PCR COVID test results. The proposed NB-ML model can be helpful in accessing different levels of information from hemogram results, through inferring non-evident patterns and parameter relationships from this data.

Bayesian techniques are based on the choice of a *prior* probability of an event (in present case, positive result for qRT-PCR test). The method considers actual evidence (hemogram data) to result in a *posterior* probability of the outcome (prediction of a positive result). By changing the selected *prior* probability, we can derive an uncertainty analysis of the model to understand its distribution. Uncertainty can be then applied to adequately adapt the classifier to a particular ongoing context. This option allows the evaluation of different decision-making scenarios concerning diverse aspects of pandemics management. During a crisis situation, measures should be taken seeking to maximize benefits and achieve a fair resource allocation . To illustrate the model flexibility and how it can be used to help on this matter, a general framework of application is proposed, followed by a simulation of four scenarios where resource scarcity is assumed.

### Application Framework

The proposed NB model can be applied in two distinct situations. When clinical data is available for a particular patient, it is highly recommended that medical staff determine the *prior* probability on a case-by-case basis. When no clinical or medical data is available, or when decisions regarding resource management involving multiple symptomatic patients are necessary, the model can be used in multiple individuals simultaneously, aiming to identify those with higher probabilities of presenting positive qRT-PCR results.

### Individual Assessment

Individual risk management and personal evaluation is essential for COVID-19 response . Individuals presenting COVID-19 symptoms are medically evaluated where no COVID-19 test is available for appropriate diagnosis confirmation. Medical practitioners can determine a probability of disease based on anamnesis, symptoms, clinical exams, laboratory findings and other available data. This probability of infection, as determined by the physician or medical team, can be considered as the *prior* probability. Using hemogram data as input, and informing the *prior* probability of COVID-19 based on medical findings, the model will consider hemogram data to inform a *posterior* probability, which can be higher or lower than the original, and based on the hemogram alterations caused by the virus infection. It is important that hemogram data would not be included in original medical assessment and *prior* determination, in order to avoid bias and reduce model overfit.

### Multiple Patients Evaluation

It can be used in situations where decisions are necessary for resource management including multiple individuals. Choice of a target group (positive or negative qRT-PCR result prediction) should be defined. The model can be applied to multiple individuals simultaneously, with the choice of *prior* probability carefully adjusted to result in a specific number of predicted individuals from the target group, according to the desired outcome. This method increases the correct selection of candidates belonging to the target group, when compared to random selection. When additional clinical data is available, or become available later, patients selected during bulk evaluation should be reassessed individually as proposed in the general framework, in order to reduce misclassifications.

### Applications to Scarcity Scenarios:

Examples of proposed model use are presented for some specific scarcity scenarios in Table . As can be seen, the model sensitivity can be adjusted by selecting *prior* probability employed, according to desired outcome or interest group. *prior* selection should be carefully decided, based on current context or situation proposed, and must consider the classification group where higher accuracy is intended.

```latex
\begin{sidewaystable}\caption{\normalsize Strategies for {NB-ML} model applications and symptomatic patient selection in scarcity conditions. Hemogram test results are available for all symptomatic patients. Scenarios proposed for situations where test results are not available (no testing or waiting {qRT-PCR} test results). Prediction results were appraised in a binary form, with positive or negative classification based on \textit{posterior} probability threshold of 0.5. Results are presented in reference to random patient selection.}\resizebox{\columnwidth}{!}{

\label{tb2}
\begin{tabular}{cccccccc}
\hline\hline
\multirow{4}{*}{\textbf{Condition}}                                                        & \multirow{4}{*}{\textbf{Context example}}                                                                                & \multirow{4}{*}{\textbf{Objective}}                                                                                & \multirow{4}{*}{\textbf{Strategy}}                                                                                & \multirow{4}{*}{\textbf{Action}}                                                                                & \multirow{4}{*}{\textbf{\begin{tabular}[c]{@{}c@{}}Starting / \\ fixed \\ \textit{prior}\end{tabular}}} & \multirow{4}{*}{\textbf{\begin{tabular}[c]{@{}c@{}}Results in training \\ set (positive \\ misclassified \\ among cleared)\end{tabular}}} & \multirow{4}{*}{\textbf{\begin{tabular}[c]{@{}c@{}}Results in test set\\ (positive \\ misclassified \\ among cleared)\end{tabular}}} \\
                                                                                &                                                                                &                                                                                &                                                                                &                                                                                &                                                                                &                                                                                &                                                                                \\
                                                                                &                                                                                &                                                                                &                                                                                &                                                                                &                                                                                &                                                                                &                                                                                \\
                                                                                &                                                                                &                                                                                &                                                                                &                                                                                &                                                                                &                                                                                &                                                                                \\ \hline\hline
\multirow{4}{*}{\begin{tabular}[c]{@{}c@{}}Testing \\ shortage\end{tabular}}               & \multirow{4}{*}{\begin{tabular}[c]{@{}c@{}}Testing capacity \\ is limited to \\ a third of \\ candidates only\end{tabular}}               & \multirow{4}{*}{\begin{tabular}[c]{@{}c@{}}Maximize \\ number \\ of infected \\ patients tested\end{tabular}}                      & \multirow{4}{*}{\begin{tabular}[c]{@{}c@{}}Prioritize \\ TP \\ identification\end{tabular}}                                        & \multirow{4}{*}{\begin{tabular}[c]{@{}c@{}}Fine-tune \textit{prior} \\ until positive \\ reach testing \\ capacity\end{tabular}}                               & \multirow{4}{*}{0.5}                                                                           & \multirow{4}{*}{\begin{tabular}[c]{@{}c@{}}130\% increase \\ in actual infected \\ patients tested\\ (\textit{prior}=0.3482)\end{tabular}}         & \multirow{4}{*}{\begin{tabular}[c]{@{}c@{}}100\% increase in \\ actual infected \\ patients tested \\ (\textit{prior}=0.9607)\end{tabular}}   \\
                                                                                &                                                                                &                                                                                &                                                                                &                                                                                &                                                                                &                                                                                &                                                                                \\
                                                                                &                                                                                &                                                                                &                                                                                &                                                                                &                                                                                &                                                                                &                                                                                \\
                                                                                &                                                                                &                                                                                &                                                                                &                                                                                &                                                                                &                                                                                &                                                                                \\ \hline
\multirow{4}{*}{\begin{tabular}[c]{@{}c@{}}Lack of \\ essential \\ workforce\end{tabular}} & \multirow{4}{*}{\begin{tabular}[c]{@{}c@{}}Professionals with\\ high risk of nocosomial \\ or  work-related \\ transmission\end{tabular}} & \multirow{4}{*}{\begin{tabular}[c]{@{}c@{}}Keep symptomatic, \\ non-infected \\ essential \\ workers in duty\end{tabular}}          & \multirow{4}{*}{\begin{tabular}[c]{@{}c@{}}Search for evident \\ non-infected \\ workers (TN \\ identification)\end{tabular}}      & \multirow{4}{*}{\begin{tabular}[c]{@{}c@{}}All workers are\\  considered as \\ infected, unless \\ model says otherwise\end{tabular}}                 & \multirow{4}{*}{0.9999}                                                                        & \multirow{4}{*}{\begin{tabular}[c]{@{}c@{}}19.4\% of total \\ workforce \\ cleared (0\%)\end{tabular}}                                    & \multirow{4}{*}{\begin{tabular}[c]{@{}c@{}}50\% of total \\ workforce \\ cleared (6.5\%)\end{tabular}}                               \\
                                                                                &                                                                                &                                                                                &                                                                                &                                                                                &                                                                                &                                                                                &                                                                                \\
                                                                                &                                                                                &                                                                                &                                                                                &                                                                                &                                                                                &                                                                                &                                                                                \\
                                                                                &                                                                                &                                                                                &                                                                                &                                                                                &                                                                                &                                                                                &                                                                                \\ \hline
\multirow{4}{*}{\begin{tabular}[c]{@{}c@{}}Lack of \\ essential \\ workforce\end{tabular}} & \multirow{4}{*}{\begin{tabular}[c]{@{}c@{}}Professionals with \\ medium to \\ low risk of \\ transmission\end{tabular}}                   & \multirow{4}{*}{\begin{tabular}[c]{@{}c@{}}Keep symptomatic, \\ non-infected \\ essential\\  workers in duty\end{tabular}}          & \multirow{4}{*}{\begin{tabular}[c]{@{}c@{}}Find ideal balance \\ to simultaneously, \\ maximize both \\ TN and TP\end{tabular}}     & \multirow{4}{*}{\begin{tabular}[c]{@{}c@{}}Use intersection of \\ sensitivity and \\ specificity curves \\ from training set\end{tabular}}            & \multirow{4}{*}{0.2933}                                                                        & \multirow{4}{*}{\begin{tabular}[c]{@{}c@{}}69.0\% of total \\ workforce \\ cleared (5\%)\end{tabular}}                                    & \multirow{4}{*}{\begin{tabular}[c]{@{}c@{}}81.5\% of total \\ workforce \\ cleared (6.6\%)\end{tabular}}                             \\
                                                                                &                                                                                &                                                                                &                                                                                &                                                                                &                                                                                &                                                                                &                                                                                \\
                                                                                &                                                                                &                                                                                &                                                                                &                                                                                &                                                                                &                                                                                &                                                                                \\
                                                                                &                                                                                &                                                                                &                                                                                &                                                                                &                                                                                &                                                                                &                                                                                \\ \hline
\begin{tabular}[c]{@{}c@{}}Limited \\ medical \\ access\end{tabular}                       & \begin{tabular}[c]{@{}c@{}}Medical assistance \\ limited to \\ 20\% of symptomatic \\ individuals only\end{tabular}                       & \begin{tabular}[c]{@{}c@{}}Avoid contagion \\ exposure of non-infected \\ patients in  ER during \\ medical assistance\end{tabular} & \begin{tabular}[c]{@{}c@{}}Eliminate non-infected \\ from candidates for \\ medical assistance \\ (TN identification)\end{tabular} & \begin{tabular}[c]{@{}c@{}}Fine-tune \textit{prior} to select \\ most likely negative results. \\ Select remanining set \\ for medical assistance\end{tabular} & 0.5                                                                                & \begin{tabular}[c]{@{}c@{}}35.6\% decrease in \\ non-infected \\ patients exposure \\ (\textit{prior}=0.0954)\end{tabular}                         & \begin{tabular}[c]{@{}c@{}}18.8\% decrease in \\ non-infected \\ patients exposure \\ (\textit{prior}=0.4652)\end{tabular}                    \\ \hline\hline
\multicolumn{8}{l}{TP: True Positive; TN: True Negative}                                                                                \\ \hline\hline
\end{tabular}}
\end{sidewaystable}
```

High accuracy in qRT-PCR result prediction is achieved based on hemogram information only. Further analysis performed on the original data (not shown) suggest that additional clinical results can improve prediction efficiency. This conclusion is in accordance with previous findings suggesting biochemical and immunological abnormalities, in addition to hematologic alterations, can be caused by COVID-19 disease . In this context, the relevance of data employed to generate ML models is emphasized. The use of large and comprehensive datasets, containing as much information as possible regarding clinical and laboratory findings, symptoms, disease evolution, and other relevant aspects, is crucial in devising useful and adequate models. The development of nationwide or regional databases based on local data is essential, in order to capture epidemiological idiosyncrasies associated with such populations . Also, natural differences in hemogram results from distinct demographic groups (as seen in reference values according to age, sex, or other physiological factors) can aggregate noise to the model, which can be reduced when large database are employed in model construction, and results can be devised for each ethnographic strata.

Despite having high overall accuracy, performance metrics obtained with proposed model show unequal ability to predict positive or negative results. This situation is caused by a significant imbalance in number of samples belonging to each of this qRT-PCR result groups in original data. The use of balanced data in machine learning model design is important to assure high prediction quality . The option of maintaining original data in model construction was adopted, since it better represents actual COVID-19 prevalence among symptomatic patients, and therefore seems to represent a more realistic situation. Additional simulations applying a balanced model (data not shown) using positive group oversampling (to compensate its insufficiency in original data) have devised alternative models with superior prediction power. Alternative balanced model results are presented in Supplementary Material Figure S2.Therefore, additional positive samples will be added to the data and used in future model versions.

As a perspective, collection of hemogram results from asymptomatic patients (in addition to symptomatic individuals) can be used to evaluate the utility of this approach on the detection of asymptomatic infections, in order to provide alternatives in diagnostics, especially in a context of testing deficiency. A web-based application was developed by the authors, in which hemogram data can be introduced for a single individual, along with prior probability of infection, based on data used to generate the present model. The online tool is available at <http://sbcb.inf.ufrgs.br/covid>. Future implementation will allow the upload of multiple patients simultaneously, and construction or testing of user data-derived models. This service will allow easy access and practical application of the proposed model.

Supplementary Material


abstract: We studied the electronic structure of a Shastry-Sutherland lattice system, HoB$_4$ employing high resolution photoemission spectroscopy and *ab initio* band structure calculations. The surface and bulk borons exhibit subtle differences, and loss of boron compared to the stoichiometric bulk. However, the surface and bulk conduction bands near Fermi level are found to be similar. Evolution of the electronic structure with temperature is found to be similar to that observed in a typical charge-disordered system. A sharp dip is observed at the Fermi level in the low temperature spectra revealing signature of antiferromagnetic gap. Asymmetric spectral weight transfer with temperature manifests particle-hole asymmetry that may be related to the exotic properties of these systems.
author: Deepnarayan Biswas$^1$; Nishaina Sahadev$^1$; Ganesh Adhikary$^1$; Geetha Balakrishnan$^2$; Kalobaran Maiti$^1$
date: 2024-09-30
title: Spectral evolution in a Shastry-Sutherland lattice, HoB$_4$

# Introduction

Physics of disorder is one of the most challenging issues in material science as it is hard to avoid disorder in real systems arising from defects, imperfections *etc.* While it is hard to probe these natural effects, some lattice structures possess intrinsic disorder on the basis of electronic interactions leading to exotic properties providing pathways to probe disorder induced effects on the electronic properties. Shastry-Sutherland lattice is one such realization, where the strong diagonal exchange interaction, $J_2$ forms spin-dimers. Competition between $J_2$ and inter-dimer interactions, $J_1$ leads to interesting phase diagram involving spin-liquid phase, long-range ordered phase *etc.* Plethora of exotic phases could be derived via tuning some control parameters. E.g. J. Liu *et al.* showed that charge carrier doping (electron or hole) in such systems leads to strong asymmetry in the electronic properties - while hole doping renormalizes spin-spin interactions leading to superconductivity, electron doped systems can be metallic but non-superconducting. Varieties of fields have emerged based on these spin-disordered parent systems. Experimental realization of these effects are not so common due to paucity of materials exhibiting such structures. Most of the studies so far are centered around the compound, SrCu$_2$(BO$_3$)$_2$, where Cu$^{2+}$ has spin-half state and forms the Shastry-Sutherland lattice exhibiting interesting phenomena.

Recently, it was observed that rare-earth tetraborides possess structure and magnetic interactions akin to Shastry-Sutherland lattice as shown in Fig. 1. The diagonal magnetic interaction, $J_2$ between the rare earth moments shown in the figure are strong forming spin-dimers. These materials exhibit diversified electronic and magnetic properties depending on interaction, $J_1$ among these dimers relative to $J_2$. For example, CeB$_4$ and YbB$_4$ do not exhibit long range order, while PrB$_4$ orders ferromagnetically at low temperatures and most of the other tetraborides order antiferromagnetically. Some of them exhibit multiple magnetic phase transitions at low temperature. For example, TbB$_4$ exhibits non-collinear magnetic structure with two antiferromagnetic phase transitions. Non-collinear magnetism with quadrupole strain and multiple phase transition is observed in DyB$_4$.

HoB$_4$ is one among the complex tetraborides exhibiting multiple magnetic transitions and non-collinear magnetism. The crystal structure of HoB$_4$ is tetragonal with $a$ = 7.085 Å and $c$ = 4.004 Å. Depending upon the symmetry of the boron sites, there are three non-equivalent borons: B1 (apical) and B3 (equatorial) form octahedra, and B2 forms boron dimers connecting the boron octahedra. The interaction among Ho moments are mediated by the borons. The spin-spin interactions in the $ab$-plane as shown in Fig. 1(b) correspond to Shastry-Sutherland type lattice. HoB$_4$ exhibits two antiferromagnetic transitions - at T$_{N1}$ = 7.1 K, it exhibits an incommensurate antiferromagnetic ordering. The magnetic structure becomes commensurate at T$_{N2}$ = 5.7 K. Resistivity and transport measurements in the presence of magnetic field exhibit anomalous phase diagram involving varieties of metamagnetic phase. Evidently, HoB$_4$ is complex exhibiting plethora of interesting phases and one of the ideal systems to study disorder induced effects in the electronic structure. Here, we studied the electronic structure of HoB$_4$ employing high resolution photoemission spectroscopy and *ab initio* band structure calculations. Interesting surface-bulk features are observed in the experimental spectra. The evolution of the electronic structure appears to be similar to that observed for charge-disordered systems.

# Experiment

Single crystalline sample of HoB$_4$ was grown using mirror furnace and characterized by $x$-ray diffraction, Laue pattern, transport and magnetic measurements revealing high quality of the sample. Photoemission measurements were carried out using R4000 WAL electron analyzer from Gammadata Scienta and monochromatic photon sources, Al K$\alpha$ (1486.6 eV), He I (21.2 eV) & He II (40.8 eV) radiations. The resolution of the instrument was set to 450 meV, 5 meV and 10 meV for Al K$\alpha$, He I and He II measurements respectively. An open cycle He cryostat was used to get the temperature variation from room temperature to 10 K. The base pressure of the chamber was better than $5\times10^{-11}$ torr during measurements. The sample was extremely hard with no cleavage plane. We fractured the sample by putting a top-post to expose clean surface. We also verified the sample surface by scraping using small grain diamond file. The angle integrated spectral features obtained from both type of surface preparations were identical. Reproducibility of the spectra were ensured after each trial of surface cleaning procedures.

The energy band structure was calculated using full potential linearized augmented plane wave method within the local density approximation (LDA) using WIEN2k software. The energy convergence was achieved using 1000 $k$-points within the first Brillouin zone. The tetragonal crystal structure of HoB$_4$ was considered for the calculations. Convergence was achieved for the paramagnetic ground state considering both energy and charge convergence criteria. In order to introduce electron-electron Coulomb repulsion among Ho 4$f$ electrons, $U$ into the calculation, we have employed LDA+$U$ method considering an effective electron interaction strength, $U_{eff}$ ($= U - J$; $J$ = Hund's exchange integral was set to zero following Anisimov et al.) for various values of $U_{eff}$ upto about 7 eV.

# Results and Discussions

Probing depth of the photoemission spectroscopic technique can be varied by changing the electron emission angle. Following this procedure, we studied the core level spectra using Al $K\alpha$ photon energy, where the mean escape depth, $\lambda$ for valence electrons is close to $\sim$ 20 Å. $\lambda$ can be made half of its normal emission value if the photoelectron emission angle is made 60$^o$ with respect to the surface normal. The Ho 4$d$ core level spectra shown in Fig. 2(a) exhibit identical lineshape for normal emission and 60$^o$ angled emission geometries indicating similar surface and bulk electronic structures for the electronic states associated to Ho. B 1$s$ spectra, however, exhibit modifications with the change in $\lambda$. There is an overall decrease in intensity of the B 1$s$ peak by about 10% in the 60$^o$ angled spectrum indicating loss of boron on the surface compared to bulk. B 1$s$ spectra exhibit two distinct features around 187 and 187.7 eV binding energies - the relative intensity of these features changes with the change in surface sensitivity suggesting the surface borons to be somewhat different from the bulk borons.

The valence band spectra shown in Fig. 2(c) are identical at both the emission geometries. The near Fermi level, $\epsilon_F$ part is rescaled and shown in Fig. 2(c) exhibiting unaltered lineshape with the change in surface sensitivity. This is curious as the intensities near $\epsilon_F$ are largely constituted by B 2$p$ states and surface borons exhibit subtle differences. These results suggest that the influence of disorder on the conduction electrons dominates over the translational symmetry breaking occurring at the surface.

Photoemission being a quantum mechanical process, the spectral features represent the eigenstates of the excited state Hamiltonian. If the electron correlation is finite in a system, these final states, which has a hole due to photoemission, differ significantly from the initial eigenstates and multiple features appear in the photoemission spectra associated to different level of screening of the photoholes, multiplet splitting, etc. The well screened feature is usually termed as the main peak and other features are called satellites. Conduction electrons consisting of the hybridized states of Ho (5$d$`<!-- -->`{=html}6$s$) states and B 2$p$ states are expected to be weakly correlated as both possess large degree of itineracy due to large orbital extensions. It is to note here that electron correlation among Ho 4$f$ electronic states is strong. However, these states primarily appear far away from $\epsilon_F$ making Ho (5$d$`<!-- -->`{=html}6$s$) states and B 2$p$ contributions dominant at $\epsilon_F$. Thus, satellites are not expected in this system. In contrast, the core level spectra of Ho exhibit significant final state effects with distinct satellite features. This is found in all the core level spectra as depicted in Figs. 3(a) and 3(b) for the 4$d$ and 5$p$ spectra, respectively. Presence of the satellites indicates that electron-electron interaction cannot be neglected in this system as also observed in various other 4$d$ and 5$d$ systems. Change in temperature down to 10 K does not influence the spectral lineshape.

B 1$s$ spectra, on the other hand, are significantly influenced by the temperature - a gradual shift towards lower binding energy is observed in Fig. 3(c) with the decrease in temperature, which cannot be represented by a rigid shift of the experimental spectrum. Three distinct features, A, B and C are required to simulate the experimental spectra at any temperature as shown in Fig. 3(e). The linewidth ($\sim$ 0.74 eV), energy separation ($\sim$ 0.7 eV) and relative intensity of the features B and C remain almost the same across all the temperatures studied. The peak position of the feature, A exhibit different temperature dependence as shown in Fig. 3(f). The intensity of B relative to C increases in the 60$^o$ angled emission spectrum compared to the normal emission one suggesting dominant surface contribution in the feature B. Higher binding energy of the surface boron signal compared to the bulk signal can be associated to surface reconstruction, defects, surface inhomogeneity that may lead to an effective reduction of negative boron valency.

There are three inequivalent borons in the lattice as shown in Fig. 1. Since the number of neighbors and bond-lengths are significantly different for octahedral borons (B1 and B3) and boron dimers (B2), the Madelung potential is expected to be different for these non-equivalent sites. Interestingly, the ratio of the integrated area of the peaks A and C is about 1 : 2.8, which is similar to the number ratio of B2 and B1+B3 (1 : 3). Thus, the peak A can be attributed to B2 and C to the photoemission signals from B1 and B3. B2 dimers bridge the boron octahedra possessing large delocalized character compared to the octahedral borons leading to a larger width of the feature A. The binding energy of the B 1$s$ peak increases with the increase in temperature - the increase is much less for dimer borons relative to the octahedral ones. Such binding energy change indicates enhancement of negativity and/or ionicity of the corresponding borons with the decrease in temperature. Thus, 2$p$ electrons of octahedral borons will gain significant local character at lower temperatures that could be related to the long range order in the ground state.

The temperature evolution of the valence band spectra is shown in Fig. 3(d). The spectral intensities in the vicinity of $\epsilon_F$ appears very similar within the energy resolution of Al $K\alpha$ spectroscopy. The features beyond 5 eV binding energy exhibit a small but gradual shift towards lower binding energies with the decrease in temperature (only two temperatures are shown for clarity in the figure). The energy shift is close to about 160 meV between 300 K and 10 K.

Valence band spectra collected using Al $K\alpha$ and He II photon energies are shown in Fig. 4. There are three distinct features around 5.8, 6.6 and 9.1 eV binding energies. The intensity of these features enhances significantly at Al $K\alpha$ photon energy. The comparison of the sensitivity of the peak intensities with the incident photon energy and the photoemission cross section suggests large Ho 4$f$ photoemission contribution in these peaks. The intensities near $\epsilon_F$ are primarily due to the B 2$p$ contributions.

In order to clarify these assertions further, the energy band structure of HoB$_4$ has been calculated within the local spin density approximations in full potential linearized augmented plane wave method. Different effective on site Coulomb correlation strength, $U_{eff}$ for the Ho 4$f$ electrons has been considered to reproduce the experimental spectra. The calculated results are shown in Fig. 4(b). It is found that the electron correlation strength of about 7 eV produces the Ho 4$f_\uparrow$ and 4$f_\downarrow$ partial density of states at the peak positions observed in the experimental spectrum.

The partial density of states (PDOS) corresponding to B 2$p$ and Ho 5$d$ states from the LDA calculations are shown in Fig. 5. B 2$s$ states dominate between 5 and 10 eV binding energies. Thus, the features around 5.8 and 9.1 eV in the experimental spectrum, can be attributed to Ho 4$f_\downarrow$ and 4$f_\uparrow$ contributions, respectively. The intensities around 6.6 eV is due to B 2$s$ photoemission signal. Boron 2$p$ states are dominated in the energy range -5 eV to $\epsilon_F$. B1 and B3 2$p$ PDOS have relatively larger contributions at higher binding energies possessing energy distributions similar to that of Ho 5$d$ PDOS indicating stronger hybridization between the octahedral borons with Ho relative to the coupling of Ho with dimer borons, B2. Interestingly, the B1 contribution at higher binding energies is the largest suggesting strongest coupling with Ho moments.

The changes close to $\epsilon_F$ has been investigated with high energy resolution using He I radiations. Since the electron escape depth at this energy is similar to the 60$^o$ angled emission in Al K$\alpha$ spectra, the identical spectral shape for normal emission and 60$^o$ angled emission shown in Fig. 2(c) established that the He I spectra would be representative of the bulk electronic structure. The raw data are shown in Fig. 6(a) exhibiting spectral evolution with significant spectral weight transfer across $\epsilon_F$ with temperature as expected in a system of Fermions. Interestingly, the difference spectra ($I(T) - I(300 K)$) exhibit larger change in spectral weight above $\epsilon_F$ relative to that below $\epsilon_F$ indicating a signature of particle-hole asymmetry. Intensity at the Fermi level obtained by the symmetrization of the experimental spectra ($I(\epsilon - \epsilon_F) + I(\epsilon_F - \epsilon)$) exhibits gradual decrease with the decrease in temperature. An additional sharp dip is observed in the 10 K data reflecting signature of the long range order in the ground state. Spectral density of states (SDOS) obtained by the division with the Fermi-Dirac distribution function exhibit similar lineshape with $|\epsilon - \epsilon_F|^{0.5}$ dependence at all the temperatures.

Disorder due to the charge degrees of freedom has extensively been studied theoretically and experimentally. Altshuler and Aronov showed that charge disorder induced effect in a correlated system exhibit $|\epsilon - \epsilon_F|^{0.5}$ dependence of the density of states near Fermi level. Such disorder effect has been observed in varieties of systems such as Ge$_{1-x}$Au$_x$, Ni(S,Se)$_2$, oxides, *etc*. Magnetic fluctuations often leads to a $|\epsilon - \epsilon_F|^{1.5}$ dependence of the spectral function or $|\epsilon - \epsilon_F|^{0.25}$ in the vicinity of quantum critical behavior. Evidently, the lineshape of the density of states close to the Fermi level is sensitive to the electronic interaction parameters as well as disorder induced effects. In the present case, although the disorder is with respect to the spin-spin interactions, the spectral function exhibit a 0.5 exponent, similar to that predicted in Altshuler-Aronov theory. HoB$_4$ exhibits multiple long range antiferromganetic order with non-linear magnetization direction. A change in the temperature and magnetic field influences the transport and magnetism in this system significantly suggesting intimate relation between these two phenomena. In addition, geometrical quadrupolar frustration has been observed in this compound indicating a strong spin-lattice interaction that might be the reason for $|\epsilon - \epsilon_F|^{0.5}$ dependence of the spectral lineshape. It would be interesting to see how the spectral function evolves in the proximity of spin liquid phase.

# Conclusion

In summary, the electronic structure of a non-collinear Shastry-Sutherland lattice system, HoB$_4$ is studied using high resolution photoemission spectroscopy. The results show subtle differences and loss of boron at the surface compared to the bulk. Ho core level spectra exhibit signature of satellites suggesting role of electron correlation among the conduction electrons. The valence band spectra corresponding to the surface and bulk electronic structures are found to be similar. The spectral intensity near the Fermi level possess dominant 2$p$ orbital character and exhibit disorder induced depletion of density of states. An additional sharp dip at 10 K is observed that can be attributed to the long range order in the ground state. Spectral weight transfer with temperature is found to be asymmetric with respect to $\epsilon_F$ revealing signature of particle-hole asymmetry. Spectral lineshape in the vicinity of the Fermi level in this spin-disordered system exhibit 0.5 exponent, which is similar to that found in charge disordered system.

# Acknowledgements

The authors N. S. and K. M. acknowledge financial support from the Dept. of Science and Technology, Govt. of India under the Swarnajayanti fellowship programme. G. B. acknowledges financial support from EPSRC, UK (EP/I007210/1).