Iteratively Estimating Error for the Least Absolute Shrinkage and Selection Operator

Information
Project Year: 2022
Student(s):   Shivesh Mehrotra | Yale University CT  
Mentor(s): Pierre Bellec - Statistics

We studied how to estimate out-of-sample and generalization error for Least Absolute Shrinkage and Selection Operator (LASSO) at each iteration of the Iterative Shrinkage-Thresholding Algorithm (ISTA). Our goal was to derive novel estimators for the aforementioned quantities as well as study their derivative structures. Due to the iterative nature of the problem, there was a natural connection to multitask learning. Thus the approach developed is analogous to an approach used to estimate these quantities in multitask learning. The results were confirmed with initial simulations.

L-functions and the Langlands Program

Information
Project Year: 2022
Student(s):   Max Lind | Princeton University NJ  
Mentor(s): Alex Kontorovich - Mathematics
Our goal was to introduce L-functions and explain their relevance to the Langlands conjectures. Along the way, we defined Godemond-Jacquet L-functions, which generalize all other known L-functions.

Markoff Surfaces and Strong Approximation

Information
Project Year: 2022
Student(s):   Enver Aman | Rutgers University-New Brunswick NJ  
Mentor(s): Alex Kontorovich - Mathematics
The Markoff equation is the equation x2 + y2 + z2 = 3xyz: we are interested in the solutions of this equation over Z or Fp, for a prime p. The Strong Approximation Conjecture states that, for each prime p and nonzero solution (x1; x2; x3) to the Markoff equation modulo p, there is a corresponding solution (x1'; x2'; x3') to the same equation over Z such that xj=xj' (mod p). We worked on the properties of these solutions.

Motifs in billion-edge graphs

Information
Project Year: 2022
Student(s):   David Wang | Hofstra University NY  
Mentor(s): James Abello - Computer Science
We worked on a method which allows for human visualization and exploration of large graphs (i.e. billion edge graphs). Our method relies heavily upon the core decomposition algorithm developed by Abello and Nakhimovich. Using this method, we create interactive Graph Cities which makes human interpretation of exponentially larger graphs possible. The Graph Cities interface also makes it easier to catalogue patterns and subgraphs of the parent supergraph. We also studied compelling avenues in which this work could continue.

Non-manipulable tournament rules

Information
Project Year: 2022
Student(s):   David Miksanik | Charles University (Prague, Czech Republic)   ,   Jan Soukup | Charles University (Prague, Czech Republic)  
Mentor(s): Ariel Schvartzman - DIMACS
We consider a manipulability of tournament rules on $n$ teams, where the rules select (possibly randomly) a single winner based on the result of all $\binom{n}{2}$ matches. A tournament rule is said to be $k$-SNM-$\alpha$, if no $k$ teams can increase their joint probability of winning by fixing the $\binom{k}{2}$ matches between them; hence, $k$-SNM-$\alpha$ is a measure of manipulability of a tournament rule. Previous works show that among all Condorcet-consistent (i.e., undefeated teams always win with probability 1) monotonous (i.e., no team can increase its probability to win by throwing matches) rules the best we can hope for are $k$-SNM-$\frac{k-1}{2k-1}$ rules and show several examples of rules achieving this bound for $k=2$. In the case of $k=3$, there exists only one rule achieving non-trivial non-manipulability, specifically $3$-SNM-$\frac{31}{60}$. Our main result is an existence of a new rule that is Condorcet-consistent, monotonous, $2$-SNM-$\frac{1}{3}$, and $3$-SNM-$\frac{1}{2}$. Additionally, the analysis of this rule is tight and uses a new technique with the possibility of finding similar rules achieving even stronger non-manipulability. Our second result shows that rules for tournaments on $n$ teams satisfying a stronger version of Condorcet-consistency can be extended to rules for any number of teams with only slightly worse manipulability. Finally, for every $d \geq 3$, we generalize a~random single binary elimination bracket rule to a~random single $d$-ary elimination bracket (RS$d$EB) rule (the complete binary tree is replaced by the complete $d$-ary tree). We show that, for every $k \geq 3$, the rule RS$k$EB is Condorcet-consistent, monotone, and $k$-SNM-$\alpha$ for $\alpha < 1$ depending on $k$.

On Inapproximability of Steiner Tree Computation in Hamming and Rectilinear metrics

Information
Project Year: 2022
Student(s):   Henry Fleischmann | University of Michigan-Ann Arbor MI  
Mentor(s): Karthik Srikanta - Computer Science
We consider the Steiner tree problem in Hamming space. While this problem is known to be APX-hard, no hardness of approximation factor is known. We show that the Hamming Steiner tree problem is NP-hard to approximate within a factor of 104. In showing this, we derive a structural correspondence between vertex covers of 4-regular graphs and optimal Hamming Steiner trees of particular point configurations in Hamming space. As a corollary of our results, we show that the Rectilinear Steiner tree problem is NP-hard to approximate within a factor of 104. We also show analagous results for discrete variants of the Steiner tree problem. That is, the Hamming and Rectilinear Discrete Steiner tree problems are NP-hard to approximate within a factor of 104.

Proving Descartes Theorem in Lean

Information
Project Year: 2022
Student(s):   Archana Mohandas | Massachusetts Institute of Technology MA  
Mentor(s): Alex Kontorovich - Mathematics
Developing a proof of Soddy-Gossett's Theorem that uses inversive geometry is a useful step towards advancing the Lean library and allowing more complex theorems to be proven in Lean. In this project, we produced a concise proof of Soddy-Gossett's Theorem using inversive geometry and began the implementation of the proof in Lean. Thus far, we have defined important concepts that are necessary in the proof such as the definitions of the co-radius, inversive coordinates, the inner product, and the inner product matrix. However, the more complex preliminary lemmas and the main statement of Soddy-Gossett's theorem are still in the process of being proved. A future goal after the conclusion of this project would be to complete the proof of Soddy-Gossett's Theorem in Lean and eventually extend the ideas of this proof to other theorems whose proofs can be simplified using inversive geometry.

Rate-1 Non-malleable Codes for Polysize Tampering

Information
Project Year: 2022
Student(s):   Guillermo Gamboa | Charles University (Prague, Czech Republic)  
Mentor(s): Marshall Ball - Department of Computer Science, Courant Institute of Mathematical Sciences, NYU
In this project, we attempted to construct a rate-1 compiler for a non-malleable code with respect to tampering from functions that can be computed by circuits of polynomial size. We extend the ideas presented by Ball, Dachman-Soled and Loss who presented such a code with rate 0. Ourc ompiler does not preserve the non-malleability properties of the code and we present an attack that shows this. Finally, we present some possible improvements in this direction.

Realizing a Fake Projective Plane as a Degree 25 Surface in P5

Information
Project Year: 2022
Student(s):   Zachary Lihn | Columbia University in the City of New York NY  
Mentor(s): Lev Borisov - Mathematics
Fake projective planes are smooth complex surfaces of general type with Betti numbers equal to that of the usual fake projective plane. Recent explicit constructions of fake projective planes embed them via their bicanonical embedding in P9. In this paper, we study Keum's fake projective plane (a = 7, p = 2, {7},D327) and construct an embedding of the fake projective plane in P5. We also simplify the 84 cubic equations defining the fake projective plane in P9.

Semi-external approach for fixed point decomposition

Information
Project Year: 2022
Student(s):   Jan Bronec | Charles University (Prague, Czech Republic)  
Mentor(s): James Abello - Computer Science
Decomposition of a graph into sub-graphs that are invariant with respect to degree peeling, i.e. fixed points, is a useful method of community detection in networks and visualisation of massive graphs. However, with the analysis of massive graph comes the problem of memory size limitations. Currently, present algorithms for fixed point decomposition cannot be used in the scenario where the whole graph doesn't fit into memory,since those approaches would result in too large I/O time. We present two methods for the computation of these decompositions and we also propose a different decomposition that is not based on fixed points but can be computed more efficiently.

Sequence-Dependent Geometric Modeling of Fluctuations of DNA

Information
Project Year: 2022
Student(s):   Travis Pence | University of South Carolina-Columbia SC  
Mentor(s): Wilma Olson - Chemistry and Chemical Biology
DNA minicircles of base pair length 336 have recently been used to inject genes into mice. The step-parameter model describes a DNA strand's configuration in space and is packaged into the software emDNA, which currently finds the minimal energy state for an input sequence. This project expands upon emDNA and adds the emDNA-move tool, allowing for the thermal fluctuation of constrained linear or circular sequences of DNA. This tool was developed in C++ and allows for the customization of the degree of random movement and fluctuation length. This project also compares the deformability of each unique tetramer and the coupling of step parameters.

Topological Analysis of Connectome Data

Information
Project Year: 2022
Student(s):   Iris Horng | University of Pennsylvania PA  
Mentor(s): Jie Gao - Computer Science
Topological Data Analysis is a promising approach for analyzing brain connectome data due to its ability to abstract underlying geometric structures in complex networks. The research presented in this study investigates a proper filtration that can be used to generate meaningful barcodes in order to represent the functional brain network of humans. Filtration using curvature was also used to explore the neural network of C. Elegans. Analysis of barcodes established a pattern among brain networks. Based on these visualizations, it was determined that differences in connectivity can be observed among various connectomes.

Tumor evolutionary model incorporating spatial heterogeneity

Information
Project Year: 2022
Student(s):   Sycamore Herlihy | Worcester Polytechnic Institute MA  
Mentor(s): Subhajyoti De - Rutgers Cancer Institute of New Jersey
Intratumoral heterogeneity is a feature in the majority of human cancers. The capability to observe the spatiotemporal clonal dynamics in the tumor of a patient is severely limited, raising the need for a mathematical model to further the study of ITH. We present STEDISH, a stochastic mathematical model for tumor growth incorporating the concept of spatiotemporal heterogeneity. STEDISH allows examining the patterns of spatial ITH under different models of tumor evolution and other constraints.

Two-way Communication Complexity of Weighted Maximum Cut

Information
Project Year: 2022
Student(s):   Liubov Samborska | Yale University CT  
Mentor(s): Sepehr Assadi - Computer Science
In this project, we consider the well-known Weighted Max-Cut problem in the context of two-party communication. We establish a two-way randomized communication complexity lower bound of Ω(n^2) for Weighted Max-Cut on graphs with n vertices, which differs from the trivial upper bound by only a logarithmic factor. To establish the desired lower bound, we construct a chain of communication reductions: Disjointness to 3-Coloring, 3-Coloring to NAE 3-SAT, and finally NAE 3-SAT to Weighted Max-Cut. The reductions we present are non-arbitrary, as the size of constructed instances and the amount of communication used in the reduction impact the strength of the resulting lower bound. We thus apply reductions in a white-box manner, establishing the communication complexity lower bound of Ω(n^2) for Weighted Max-Cut and additional lower bounds for the intermediate problems of 3-Coloring and NAE 3-SAT in the communication reduction chain.

Understanding non-monotonicity for two independent items

Information
Project Year: 2022
Student(s):   Jachym Mierva | Charles University (Prague, Czech Republic)   ,   David Sychrovsky | Charles University (Prague, Czech Republic)  
Mentor(s): Ariel Schvartzman - DIMACS
Sale of independent items is an often studied problem. Myerson's lemma provides an efficient way to realize an optimal trade of a single item. However, in a multi-item setting, the results can be non-intuitive. It has been shown before that in some cases, buyers having a higher valuation for sold items can lead to lower revenue of the seller. This non-monotonicity is not yet understood. In this paper, we provide conditions which guarantee monotonicity and an algorithm to generate non-monotone examples. We show that such examples are rare, even though the revenue can change significantly in some cases.

Visibility graphs of polygons

Information
Project Year: 2022
Student(s):   Gaurav Kucheriya | Charles University (Prague, Czech Republic)  
Mentor(s): James Abello - Computer Science
The visibility graph recognition problem asks to determine if for a given graph G there is a polygon P having G as its visibility graph. It is not known to be in NP. Here we show that persistent graphs having a block of consecutive blocking vertices can be realized as a polygon.We conjecture that this technique can be extended to realize a larger subclass of persistent graphs.

A note on the involutive concordance invariants of certain (1,1)-knots

Information
Project Year: 2021
Student(s):   Anna Antal | Rutgers University-New Brunswick NJ   ,   Sarah Pritchard | Georgia Institute of Technology-Main Campus GA  
Mentor(s): Kristen Hendricks - Mathematics
In this project we worked with involutive Heegaard Floer homology, developed by Hendricks and Manolescu as a refinement to Heegaard Floer homology. Involutive Heegaard Floer homology introduced two new knot concordance invariants. These invariants are interesting, because unlike other concordance invariants like τ, ε, and ν, they do not necessarily vanish on knots of finite concordance order, like the figure-eight knot. We computed these two new invariants for all non-thin and non-L-space 10 and 11-crossing (1,1)-knots for which they were not yet known. We also studied (1,1)-diagrams, and the relationships between them and the resulting CFK^∞ chain complexes.

Bayesian inference of the origin of an infection process over a network

Information
Project Year: 2021
Student(s):   Hwai-Liang Tung | Brown University RI  
Mentor(s): Min Xu - Statistics
Many processes such as the diffusion of information or fake news over a social media network or the transmission of disease may be modeled with an infection process. We investigate the problem of finding the patient zero using noisy observations of a set of infected nodes. Given a set of detected infected nodes and the background graph, we present a MCMC algorithm to calculate the probability a node is the patient zero. We also prove that for an infinite lattice Z^d and a set of n infected nodes, for any ε-level, there exists a credible set C_ε such that there are O(√ n log^(2d)n) nodes in C_ε.

Co-evolution of Opinions and Signed Network dynamics in Real World Scenarios

Information
Project Year: 2021
Student(s):   Soham Palande | Rutgers University-New Brunswick NJ  
Mentor(s): Jie Gao - Computer Science
With the pervasion of social media in modern society, evolution of opinions and network dynamics in social networks have gained special interest from industry, academia, and governments. The study of opinion and network dynamics has several applications- social sciences, behavioral economics, game theory- with consequential implications. Adversarial attacks through fake news and troll bots are increasingly common and there is a need for a rigorous understanding of their impact on individuals and the social networks. We study opinion and network dynamics and characterize social networks at their limit states extending recent work in co- evolution of signed networks and opinions to better represent the real world. We build upon intuition of the real world, in which individuals hold multiple opinions on a range of issues for example- healthcare, gun control and seek to better represent complex relationships/social ties between individuals in a network. We formulate these intuitions mathematically and incorporate them in a co-evolution model which combines signed network dynamics based on structural theory and opinion dynamics. We conduct a comprehensive study through numerous simulations and validate the model on the benchmark Zachary's Karate Club dataset to achieve 100% accuracy. We mathematically address divergence of social ties to infinity, which in the real world would translate to limitless polarization, by introducing normalization which keeps the relative magnitude of the social ties intact but prevents divergence.

Combinatorial Options Markets

Information
Project Year: 2021
Student(s):   Jacob Gorenburg | Haverford College PA  
Mentor(s): David Pennock - DIMACS
This paper is an extension of the work done by Xintong Wang, David Pennock, and others on combinatorial options markets. Their original paper focused on the design and time complexity of an exchange that accepts call and put options on multiple stocks in a single bundle. These options give traders greater flexibility and precision in their trading. The original work focused on determining the optimal match in single-instance auctions. We extend this in two ways: providing a system for transferring the surplus from these trades back to the traders and designing a continuous combinatorial auction. One of the key principles of modern exchanges is that their profits come from a small fee on every trade, not from any surplus that occurs during the course of trading. If one trader is willing to buy a stock at a given strike price for $10 and another is happy to sell for $7, the $3 surplus is given to whichever trader placed the second bid. In a normal market, the surplus is always a flat amount since identical options are being bought and sold. However, in a combinatorial market, the surplus potentially also contains a variable surplus that depends on the relative successes of the various options. In this paper, we propose a method of transferring the surplus in a combinatorial market from the exchange to a trader by providing them with a bundle of options. In addition to dealing with the surplus, this paper extends the idea of a combinatorial market to the continuous-style auctions found in real-world exchanges. Wang's work focuses on a single-instance auction where all offers are collected and a single optimal match is returned. In practice, most markets are running continuously and decide on trades when new bids are placed. This style of auction allows for faster trading and benefits quick traders over those with better offers. In this paper we introduce a mechanism for a continuous combinatorial market and make some progress on an algorithm for running one.

Comparing Physics-Informed Loss Functions for Porosity Prediction in Laser Metal Deposition

Information
Project Year: 2021
Student(s):   Erin McGowan | Rutgers University-New Brunswick NJ  
Mentor(s): Weihong 'Grace' Guo - Industrial and Systems Engineering
Laser metal deposition (LMD) is a type of additive manufacturing (AM) during which metal components are created using a laser beam that fuses metal powder by melting it as it is deposited. Porosity, or small cavities that form in this printed structure, is generally considered to be one of the most destructive defects that can occur in metal AM. Currently, porosity can be measured after printing with CT scans. While this is useful for understanding the nature of pore formation and its characteristics, purely physics-driven models lack real-time prediction ability. Meanwhile, a purely deep learning approach to porosity prediction leaves valuable physics knowledge behind. Here we create a hybrid model that takes advantage of both empirical and simulated LMD data to show how various physics-informed loss functions impact the accuracy, precision, and recall of a baseline deep learning model for porosity prediction. In particular, we find that some versions of the physics-informed model are able to improve upon the precision of the baseline deep learning-only model (albeit at the expense of overall accuracy). This work lends itself to a wide array of applications as LMD is frequently used in the automotive, aerospace, energy, petrochemicals, and medical industries.

Data Analysis on COVID-19 patterns

Information
Project Year: 2021
Student(s):   Yun Huen Cheng | Mount Holyoke College MA  
Mentor(s): Lazaros Gallos - DIMACS
The COVID-19 pandemic has a global impact due to its high infectivity and due to human travel behaviors. Focusing on geographic patterns to analyze the evolution of the pandemic, we can determine which areas are at risk of an outbreak. By analyzing the spatial correlations of new active cases in the US at the county level we can show the extent of these correlations at different times. Our results show that the epidemic was largely not uniform but formed clusters. We found that in the first few months long-range spreading was mainly located in cities rather than rural areas. In addition, there exists a percolation transition in November 2020, and a smaller transition in January 2021, both corresponding to peaks of the epidemic.

Data-Driven Dynamics

Information
Project Year: 2021
Student(s):   Anna Cusenza | University of California-Los Angeles CA  
Mentor(s): Konstantin Mischaikow - Mathematics
Our goal is to interpret the dynamics of multi-parameter nonlinear dynamical systems. Such systems can have very complicated behavior. Because of this and the fact that in applications data collected may not accurately represent the generating system, it is desired to gather from the data robust patterns which are consistent with respect to perturbations in parameter values. A method to analyze data generated by such systems, as introduced in previous work, is to discretize the dynamics via a decomposition of the compact phase space into square grids. With a grid decomposition of the phase space, a directed graph which approximates the behavior of the system for a fixed set of parameter values is created. Each node in the directed graph represents a cubical cell in the phase space, and the directed edges denote where each cell maps to under the system. With graph algorithms we can find regions which contain recurrent dynamics, and using Conley index theory we can further identify special properties of the system. In this project we use the method of analysis described above, with the change being that we decompose the phase space using Voronoi cells rather than cubical cells. A Voronoi cell complex is generated by a finite set $P = \{x_1,x_2,\dots, x_n\}$, where each $x_i$ is chosen in the phase space. Each Voronoi cell $V_i$ is defined by $V_i = \{x|d(x,x_i)\leq d(x,x_j) \text{ for all } j\neq i,\ 1\leq j\leq n\}$. Our goal is to determine whether using Voronoi cells to decompose the phase space yields comparable results to the results generated by the cubical complex; in visual quality and accuracy. The goal is to apply these methods to real data collected from systems generated by weather patterns and robotics.

Data-Guided Modeling of Cancer and Intratumor Heterogeneity

Information
Project Year: 2021
Student(s):   Sivasomasundari Arunarasu | Emory University GA  
Mentor(s): Subhajyoti De - Rutgers Cancer Institute of New Jersey
Although all cells in a tumor come from one original mutated cell, there is genetic, nongenetic, and microenvironmental variation among these cells, referred to as intratumor heterogeneity. Nongenetic heterogeneity, which includes epigenetic changes, is closely related to the cancer stem cell (CSC) model, which states that certain tumor cells are characteristic of regular stem cells, and are likely to contribute to tumor survival and regrowth after treatment. Being able to identify phenotypic features that can classify tumor cells into a stemlike and differentiated state can help to better understand the CSC model, and the characteristics of cancer stem cells. Using the data of cells grown from the T24 cell line, PCA (principal component analysis) and UMAP analysis in R was used to group tumor cells into two different classes, which are thought to correspond with more holoclone (stemlike) and more differentiated cells. This classification can be used to analyze the colony makeup of other tumors, and might help determine the likelihood that a tumor will regenerate after treatment, based on the proportion or density of holoclone and meroclone-like cells.

Deep Genomic Analysis of Tumor Specimen

Information
Project Year: 2021
Student(s):   Yixuan (Connie) Zhang | Macalester College MN  
Mentor(s): Hossein Khiabanian - Rutgers Cancer Institute of New Jersey
Understanding complexity, dynamics, and stochastic patterns in genomic data - concepts native to physics and mathematics - is critical for elucidating how disease states originate and evolve. In this project, we focus on the application of the Tunable Biclustering Algorithm (TuBA) to examine genetic and clinical data of cancer patients, aiming to identify genetic markers that can provide us insight into how cancers evolve. Our ultimate goal is to develop novel statistical platforms for fast translation of genomic data into clinical practice.

Explorations in Tensor Phase Retrieval

Information
Project Year: 2021
Student(s):   Cheng-Hao Fu | University of Michigan-Ann Arbor MI  
Mentor(s): Anand Sarwate - Electrical and Computer Engineering
This paper outlines a project that started in the summer of 2021 with the goal of coming up with a better method than existing ones to solve the Low Rank Phase Retrieval Problem, primarily in terms of sample complexity, but also potentially in runtime. It remains a work in progress.

Graph Cities Experimentation

Information
Project Year: 2021
Student(s):   Jiajun Chen | Stony Brook University NY   ,   Moriah Walker | Rutgers University-New Brunswick NJ  
Mentor(s): James Abello - Computer Science
The applications Graph Cities and Graph Strata are able to decompose graphs of over 1 billion edges into smaller subgraphs that can fit on a screen and be interactively explored without losing the context of the original graph. Observing these local subgraphs reveals a variety of patterns. This report serves to find a way to describe these patterns as well as unpack meaning from these patterns by generating data stories and detecting semantically important vertices through the vertex statistics, activity and diversity.

Health Disparities & COVID-19

Information
Project Year: 2021
Student(s):   Fiona Shafer | Rutgers University-New Brunswick NJ  
Mentor(s): Christie Nelson - Masters of Business and Science in Analytics and CCICADA
As the COVID-19 pandemic progresses across the nation and the world, data suggests the virus disproportionately impacts marginalized communities. Specifically, disparities exist in infection rates and disease outcomes by race and ethnicity in the United States. However, it is unknown to what degree these differences have impacted specific communities in the past year. During the course of this research project, I intend to research various avenues related to COVID-19 including but not limited to long term care facilities and education systems.

Identifying Bifurcations for Combinatorial Dynamical Systems

Information
Project Year: 2021
Student(s):   Yousef Sayes | New Jersey Institute of Technology NJ  
Mentor(s): Konstantin Mischaikow - Mathematics
The identification of higher order bifurcations in combinatorial dynamical systems allows for deep insights into the underlying structures of these systems. The modeling of such phenomena, specifically biological ones, in a combinatorial fashion has been explored; however, a systematized process to which higher order bifurcations may be found through is still nonexistent. The use of multiple paths between given adjacent parameter nodes and applying slow-fast dynamics to represent transitions between these nodes point us toward the discovery of such bifurcations. A process to identify bifurcations therefore allows us far greater insight into these systems, especially as their complexity increases.

Lower Bounds for Deterministic Graph Coloring in the Streaming Model

Information
Project Year: 2021
Student(s):   Glenn Sun | University of California-Los Angeles CA  
Mentor(s): Sepehr Assadi - Computer Science
In the graph semi-streaming model of computation, the input graph arrives as a stream and the goal is to compute solutions in O(n polylog n) space. With randomization, Assadi, Chen, and Khanna recently showed that the coloring problem for graphs of maximum degree Δ admits a solution in this model using an optimal Δ + 1 colors. However, little is known about this problem in deterministic settings. In this project, we showed that when Δ = Ω(n^{2/3}), no deterministic streaming algorithm can output proper colorings using O(n^{1-ε}) colors, for all ε > 0. In particular, any general algorithm for deterministic graph coloring must use at least Ω(Δ^{3/2 - ε}) colors.

Machine Learning for SAT-solving heuristics

Information
Project Year: 2021
Student(s):   Emily Chin | Harvey Mudd College CA  
Mentor(s): Periklis Papakonstantinou - Management Science and Information Systems
3-SAT is a known NP-Hard problem in theoretical computer science. In other words, it takes exponential time to deterministically solve this problem. However, there is a fairly effective random algorithm that finds the solution quickly most of the time, if there is one. This "Random Walk Algorithm" finds the solution by acting randomly at each timestep, but the results may not be completely random. This leads us to explore what features about an instance of the problem are important and how we can use these features to predict if a solution can be found at the end of a certain number of iterations of the algorithm. We used various machine learning techniques to make predictions using features that we found to be significant.

Morse Flow Trees and Chain Complexes

Information
Project Year: 2021
Student(s):   Sangjun Ko | Rutgers University-New Brunswick NJ  
Mentor(s): Christopher Woodward - Mathematics
Morse flow trees are a tool to compute Legendrian contact homology, and Ekholm has showed that these flow trees have a correspondence with certain types of rigid pseudoholomorphic disks in T*M if M is a manifold. There are still some open questions regarding the space of these flow trees, some of which we investigated this summer. Our aim here will be to describe some of the background material required to understand the definition of Morse flow trees. First we discuss the elements of Morse Theory and define the Morse complex; then we discuss a little bit of symplectic topology so that we define what a Lagrangian submanifold is; finally, we give the definition of a flow tree and introduce a problem related to it.

On Spanning Tree Counts for Bipartite Graphs

Information
Project Year: 2021
Student(s):   Aaditya Raghavan | Georgia Institute of Technology-Main Campus GA  
Mentor(s): Bhargav Narayanan - Mathematics
Ehrenborg conjectured that, for a simple bipartite graph G, the number of spanning trees of G, is bounded above. The bound is known to hold with equality for a class of bipartite graphs known as Ferrers graphs, and not necessarily with equality in other, specific classes of bipartite graphs. In this paper, we introduce another inequality motivated by the Cartesian product of bipartite graphs - in particular, an inequality that relates the eigenvalues of the Laplacians of the multiplicand graphs, with products of the vertex degrees and the part sizes of the factor graphs. Furthermore, we introduce an inductive framework through which the conjecture can be approached by means of a determinantal inequality, along with some preliminary results that can be proved using such a technique.

On the compactification of bounded edge flow trees

Information
Project Year: 2021
Student(s):   Kenneth Blakey | Brown University RI  
Mentor(s): Christopher Woodward - Mathematics
Let L⊂J^1(M) be a closed Legendrian submanifold of the 1-jet space of a closed Riemannian manifold. We conjecture that the moduli space of flow trees with at most n edges has a natural compactification given by the moduli space of broken flow trees with at most n edges. In order to solve this one would have to compactify the various moduli spaces of gradient trajectories of local function differences determined by L. We establish a compactification of the moduli space of Morse trajectories - gradient trajectories connecting critical points of a fixed local function difference - similar to the compactification in standard Morse theory.

Predicting Dissolution Rates of Volcanic Glass Using Graph Neural Networks

Information
Project Year: 2021
Student(s):   Emily Thompson | Southwestern University TX  
Mentor(s): Shashanka Ubaru - IBM Watson Research Center
Graph neural networks (GNNs) in recent years have gained popularity due to the increasing need for machine learning models to accommodate graph structured data. Interdisciplinary fields such as material informatics use graph neural networks to conduct research on materials that are otherwise difficult to analyse. One of these materials, volcanic glass, is used in storing nuclear waste due to their durability in extreme conditions. This durability, measured by its dissolution rate, is difficult to determine in a traditional lab environment due to the extensive time and resources needed to collect and analyse samples. I seek to implement a GNN model that performs regression in order to predict the dissolution rate of ten different types of volcanic glass. In addition to implementing the GNN, I explore varying methods of optimizing the performance of the model. Results demonstrate that unknown glass materials and their accompanying dissolution rate can accurately be determined the GNN in a short period of time. The long term goal for this research is to be able to discover new glass materials based on their dissolution rates.

Programmable Routers: RMT and P4

Information
Project Year: 2021
Student(s):   Tavis Johnson | Iona College NY   ,   Delta Lyczak | Lafayette College PA  
Mentor(s): Srinivas Narayana - Computer Science
Internet routers deal with significant amounts of internet traffic in everyday use. This increasing amount of traffic over the internet affects the performance, and speeds, at which these routers operate. The goal of our research project is to understand how routers manage processing packets at such high speeds while remaining flexible and reconfigurable. Through this project, we have come to understanding layers of networking-abstraction, key networking concepts, pipelined systems, along with the history of routing software, hardware, and information on completed P4 programming exercises. With this knowledge, we hope to be able to design improvements for the performance of routers and the management of campus network resources.

Simple reductions to circuit minimization

Information
Project Year: 2021
Student(s):   Noah Singer | Harvard University MA   ,   Vishal Ramesh | Charles University (Prague, Czech Republic)   ,   Sasha Sami | Charles University (Prague, Czech Republic)  
Mentor(s): Eric Allender - Computer Science

The complexity of the minimum circuit size problem (MCSP) - and its many variants, including the minimum Kolgomorov time-complexity problem (MKTP) - is linked intricately to countless other questions in theoretical computer science. For instance, NP-hardness of MCSP or MKTP is known to imply ZPP ≠ EXP, and even if such reductions exist, they cannot be nearly as simple as the standard NP-completeness reductions for other problems. In this project, we study whether various variants of circuit minimization can be complete for NP, or smaller classes, under simple types of reductions. Specifically, we investigate the following questions. Can recent hardness results for MKTP be extended to MCSP? How important is adaptivity in NP-hardness reductions for circuit minimization? How robust is the recent definition of non-interactive statistical zero-knowledge with logspace-bounded verifiers and simulators?

Tetramer-Dependent Elastic Energy Model for DNA Conformation

Information
Project Year: 2021
Student(s):   Zoe Wefers | McGill University  
Mentor(s): Wilma Olson - Chemistry and Chemical Biology
DNA folding plays an important role in gene regulation and genome packaging. The ability for a piece of DNA to bend is largely impacted by the elastic energy stored between base pairs in a segment of DNA. There are well-established techniques to mathematically model the elastic energy of a collection of DNA base pairs. It has long been known that the intrinsic elasticity parameters of DNA are sequence-dependent, but past work only models DNA folding using elasticity parameters that are specific to dimers, a sequence of two base pairs. Recent studies have suggested the intrinsic elasticity parameters of DNA are actually unique to tetramers, a sequence of four base pairs. We present software that can simulate DNA folding by optimizing elastic energy using tetrameric steps. With this new software, we optimized the configuration of DNA minicircles and confirmed that increasing the scope of sequence dependence from dimers to tetramers has an effect on DNA conformation.

The Quantum Mechanics of an Electron Surrounded by Two Photons

Information
Project Year: 2021
Student(s):   Hallie Byatt | Georgia Institute of Technology-Main Campus GA  
Mentor(s): Shadi Tahvildar-Zadeh - Mathematics
I formulated a relativistic Lorentz-covariant system of wave equations for a three-body, three-time one-dimensional system consisting of two photons and one electron. These particles interact only upon contact, where a boundary condition is applied to preserve the conservation of probability current. The solutions for various configurations of the initial-boundary value problem were found, and hopefully will pave the way to finding the solutions for a system of one electron surrounded by N photons.

Toward a fair and strategyproof tournament rule for selfish agents

Information
Project Year: 2021
Student(s):   Eric Xue | Yale University CT  
Mentor(s): Ariel Schvartzman - DIMACS

A tournament on n agents is a binary relation that describes the outcomes of the $\binom{n}{2}$ matches played between two distinct agents. Tournaments play an important role not only in sporting events but in any setting where any two agents are comparable and some "most qualified" agent should be selected, such as elections. The winner of a tournament is determined by a tournament rule that maps tournaments to probability distributions over the agents. We want these rules to be fair (i.e., choose a high quality agent) and robust to strategic manipulation. Prior work has shown that under minimally fair rules, manipulations between two agents can be prevented when utility is nontransferable but not when utility is completely transferable. Motivated by the fact that neither utility model reflects reality, we consider pairwise manipulations under partially transferable utility. These situations arise, for example, when agents care about winning themselves but are still willing to collude if joint gains offset perceived individual losses. We show that previously studied tournament rules require high levels of selfishness in order to prevent pairwise manipulations and that for fixed levels of selfishness, there are always two agents who stand to gain a constant positive amount of probability under these rules by manipulating in large enough tournaments. We also show that for stronger notions of fairness, non-manipulable tournament rules are closely related to tournament rules that witness decreasing gains from manipulation as the number of agents increases. A major open question that remains is whether there is a tournament rule that is fair and non-manipulable for a fixed level of selfishness that meets the lower bound we prove. We conjecture that indeed such a rule exists.

Unique Solution to the House Assignment Problem for N = 4

Information
Project Year: 2021
Student(s):   David Ryzak | Charles University (Prague, Czech Republic)   ,   David Sychrovsky | Charles University (Prague, Czech Republic)  
Mentor(s): Ron Holzman - Mathematics, Princeton University
We study a problem from mechanism design where we try to find a mechanism which assigns N houses to N players based on their preferences with desirable properties. These aim to unsure fairness and efficiency while limiting exploitation by the players and give incentives to tell the truth. We aim to show that the Random serial dictatorship (RSD) is the only mechanism with these prescribed properties. It has been shown before that RSD in unique for N=3. We provide a computer aided proof of uniqueness for N=4 as well as rigorous for selected preference profiles for arbitrary N.

Δ-Vertex Coloring in Graph Streams: Towards a Streaming Brooks' Theorem

Information
Project Year: 2020
Student(s):   Pankaj Kumar | Charles University (Prague, Czech Republic)   ,   Parth Mittal | Charles University (Prague, Czech Republic)  
Mentor(s): Sepehr Assadi - Computer Science
The celebrated Brooks' theorem in graph theory states that every connected graph which is not a clique nor an odd-cycle can be Δ colored; here Δis the maximum degree of the graph and we use n to denote the number of vertices. We investigate the possibility of obtaining a Δ-coloring via a space-efficient streaming algorithm. That is, an algorithm that makes one pass over the edges of the input graph while using a limited memory - much smaller than the input size which can be as large as Θ(n2) - and at the end of the stream, outputs a coloring of the input graph. Our main result is a randomized algorithm for this problem that uses O(n7/4) space.

A mathematical model of pancreatic cancer development and the immune response

Information
Project Year: 2020
Student(s):   Chloe Shiff | Brandeis University MA  
Mentor(s): Subhajyoti De - Rutgers Cancer Institute of New Jersey
Although tumors arise as a result of the rapid proliferation of cancerous cells, these cells comprise only a fraction of the tumor. Many types of cells, each of which possesses a specific role in the growth of solid tumors, make up the tumor microenvironment (TME). In recent years, immunotherapy, a treatment that utilizes the significant presence of immune cells within the TME, has proven effective for liquid tumors such as leukemia and lymphoma. However, in solid tumors, the complex relationships within the TME can make the response to immunotherapy challenging to predict and the treatment difficult to use effectively. Thus, we created a mathematical model of pancreatic adenocarcinoma growth which includes cancer cells, T-cells, and macrophages, and used it to study the dynamics of the cellular interactions within the TME over time and throughout treatment. The model behavior indicates that the ability of CAR T-cells to kill cancer cells is more critical than their persistence. However, when used in combination with chemotherapy, the tumor can be eliminated using less effective T-cells than would be needed if using immunotherapy alone. These results indicate that CAR T-cell research should focus on engineering T-cells with high efficacy in killing cancer cells rather than long persistence in the body. Doing so could allow for effective immunotherapy use in solid tumors, either as a lone treatment or used in combination with chemotherapy.

A Pandemic Model: Looking at Policies and Historic COVID Data to Promote a Resilient Future

Information
Project Year: 2020
Student(s):   Coleman Bligh | Colgate University NY   ,   Elisavet Gallou | Rutgers University-New Brunswick NJ   ,   James Robinson III | Albany State University GA   ,   Yuqing Zhou | Rutgers University-New Brunswick NJ  
Mentor(s): Ava Majlesi - Miller Center for Community Protection & Resilience/Center for Critical Intelligence Studies

In the recent light of the SARS-CoV-2 pandemic, the United States' response has proven to be incredibly poor with respect to policy decisions made on the local and federal levels. With deaths totaling at 142,755, this is more than the number of battle deaths in the Vietnam War, Korean War, and World War I combined. One of the easiest ways to prevent such an outbreak from reoccurring has shown to be a stronger implementation of policies on all governmental levels. Three categories which require more heavy implementation include Test, Track, and Treat (TTT), Distance, Isolate, Mask (DIM), and Prepare & Organize, Lead & Inform, Coordinate & Yaager results (POLICY). The American response was not sufficient enough to prevent mass infections and fatalities, and analysis of state and federal policies with respect to COVID-19 containment outstandingly point to a need for more drastic policy decisions.

AES and Expansion Properties

Information
Project Year: 2020
Student(s):   Eli Goldin | Columbia University in the City of New York NY  
Mentor(s): Periklis Papakonstantinou - Management Science and Information Systems
AES is one of the best known and most widely used block ciphers to date. However, the primary argument towards the security of AES is simply the fact that it has not been broken in practice. We present a new way of theoretically analyzing the security of AES and related block ciphers by looking at the expansion of a related Cayley graph. More specifically, we provide a generalization of a single round of AES, which we consider as a generating set for the symmetric group. This then gives a Cayley graph for which one could theoretically explicitly calculate the spectral gap. We explicitly perform this calculation for small security parameters, and we compare these results to related expansion results for the ideal cipher. We hope that this analysis, when continued, will either give some indication of the security of AES or will lead to an explicit attack.

Analyzing gene regulatory networks by comparing the dynamics obtained via DSGRN (Dynamic Signatures Generated by Regulatory Networks) and RACIPE (Random Circuit Perturbation)

Information
Project Year: 2020
Student(s):   Aaron Scheiner | Rutgers University-New Brunswick NJ  
Mentor(s): Konstantin Mischaikow - Mathematics
Understanding gene regulatory networks is key to targeting diseases with drugs that don't include dangerous risks. In this research project, we studied the analyses generated by Huang et al. using random circuit perturbation (RACIPE). After gaining a comprehensive understanding of Huang's paper, we employed RACIPE to reproduce these results. We then worked to produce analogous results in DSGRN (Dynamic Signatures Generated by Regulatory Networks). Moreover, we found the DSGRN parameter nodes corresponding to the RACIPE data. We added weights to the parameter nodes based on what percentage of the RACIPE models occurred in each parameter node. Our preliminary low dimensional results lead to enticing conjectures. When including weighting, the DSGRN results are very similar to RACIPE for high Hill coefficients. Thus, weighting in DSGRN with sampling from RACIPE yields comparable results to RACIPE with far quicker computations. Hence, with additional work we hope to show that DSGRN can achieve comparable results to RACIPE for a fraction of the computational cost of RACIPE.

Antipodal paths in 2-colorings of hypercubes

Information
Project Year: 2020
Student(s):   Tomas Hons | Charles University (Prague, Czech Republic)   ,   Marian Poljak | Charles University (Prague, Czech Republic)   ,   Tung Anh Vu | Charles University (Prague, Czech Republic)  
Mentor(s): Ron Holzman - Mathematics, Princeton University

Feder and Subi conjectured that for any 2-coloring of edges of the hypercube Qn, there always exists a pair of antipodal vertices connected by a shortest path which changes color at most once. Leader and Long proved that there always exists a path between antipodal vertices with at most n/2 changes and Dvorak improved this bound to (3/8+o(1))n. We give some partial results which may lead to further improvements of Dvorak's upper bound.

Authenticated Disk Key-Value Stores

Information
Project Year: 2020
Student(s):   Akif Patel | Rutgers University-New Brunswick NJ  
Mentor(s): Martin Farach-Colton - Computer Science
Security enclaves are sections of chips that run trusted and verified code that cannot be accessed by other parts of the computer. If code running on an enclave would like to use a key-value store that is stored on the disk, it will have to go through the untrusted main OS. We investigated how to verify integrity and maintain secrecy with regards to this key-value store that the enclave is using. We first showed how to do it in the case when our data structure is a tree. Then, we used this idea to construct a faster set membership proof on disk. Finally, we made this construction write optimized using ideas drawn from Bε-trees.

Classifying Integral Polyhedra

Information
Project Year: 2020
Student(s):   Brittany Gelb | Muhlenberg College PA  
Mentor(s): Alex Kontorovich - Mathematics
We systematically determine infinite circle packings that arise from 8 and 9-vertex polyhedra. We then classify packings where all radii are the reciprocals of integers.

Covid-19 Data Analysis and Risk Management Assessment

Information
Project Year: 2020
Student(s):   Rachael Tovar | Wheaton College MA  
Mentor(s): Christie Nelson - Masters of Business and Science in Analytics and CCICADA
In the United States alone, over 3.8 million people have contracted the Covid-19 virus. Due to the increased demand for testing, there have been shortages and issues for the tests themselves. Access to testing and the reporting of cases varies by region, so reported numbers are mainly consisting of people who exhibit symptoms of the virus. Taking into account those who are asymptomatic or do not have access to reliable tests, one can safely assume that the actual cases of this virus are actually much larger. Despite rates continuing to rise and testing becoming more scarce, the United States government pushes for the economy to reopen. This may lead to multiple consequences, to the cases of citizens surging to the negative social impact. Which makes research on methods to help protect individuals from infection a necessity. Continuing to see the changes and trends of this virus is crucial to both understanding causation of case spikes and to finding solutions to help flatten the curve as well as find effective ways to distribute resources. To understand the full scope of the Pandemic's impact on the United States, we considered a wide range of information regarding the virus. Focusing on the changing mobility patterns and state government mandates to understand the risk management of state governments as well as the social impact of Covid-19.